A method comprising: calculating an expected wellbore energy indicative of an amount of energy needed to deform a theoretical metallic cylinder into a shape traced out by a planned wellpath, wherein the metal cylinder has a length equal to a portion of the wellpath and the deforming comprises bending and/or twisting the theoretical metallic cylinder into a shape traced out by the planned wellpath, wherein the calculating is performed by a computer system and the expected wellbore energy is calculated using the equation: E exp = ∑ i = 1 n κ i 2 + τ i 2 Δ Di D n + Δ D n where E exp is the expected wellbore energy across a plurality of measurement locations n; i is an index representing measurement locations along the planned wellpath; κ i is a wellbore curvature at the measurement location i; τ i is the torsion at the measurement location i; Δ Di is a change in depth at the measurement location i; D n is a depth at survey point i; Δ D n is a change in depth across a plurality of measurement locations n; and where D n + ΔD n is a total length between depths; calculating an expected tortuosity using the equation: T = ∑ i = 1 n α n − 1 + ΔD × δ i D i − D i − 1 where T is the expected tortuosity of the planned wellpath; α is an inclination angle of the wellpath between two measurement locations, defined as angle between a tangent line along the planned wellpath and a vertical line; Δ D is a distance between a first measurement location within the planned wellpath and a second measurement location within the planned wellpath; δ i is a dogleg severity value; and D i is a depth along the planned wellpath at a survey location i; combining the expected wellbore energy with the expected tortuosity to create a first value indicative of probability of occurrence of a drilling event for drilling along the planned wellpath; and then drilling a wellbore along at least a first portion of the planned wellpath if the first value indicates the probability of the drilling event is less than a predetermined threshold probability. calculating an expected wellbore energy indicative of an amount of energy needed to deform a theoretical metallic cylinder into a shape traced out by a planned wellpath, wherein the metal cylinder has a length equal to a portion of the wellpath and the deforming comprises bending and/or twisting the theoretical metallic cylinder into a shape traced out by the planned wellpath, wherein the calculating is performed by a computer system and the expected wellbore energy is calculated using the equation: E exp = ∑ i = 1 n κ i 2 + τ i 2 Δ Di D n + Δ D n where E exp is the expected wellbore energy across a plurality of measurement locations n; i is an index representing measurement locations along the planned wellpath; κ i is a wellbore curvature at the measurement location i; τ i is the torsion at the measurement location i; Δ Di is a change in depth at the measurement location i; D n is a depth at survey point i; Δ D n is a change in depth across a plurality of measurement locations n; and where D n + ΔD n is a total length between depths; calculating an expected tortuosity using the equation: T = ∑ i = 1 n α n − 1 + ΔD × δ i D i − D i − 1 where T is the expected tortuosity of the planned wellpath; α is an inclination angle of the wellpath between two measurement locations, defined as angle between a tangent line along the planned wellpath and a vertical line; Δ D is a distance between a first measurement location within the planned wellpath and a second measurement location within the planned wellpath; δ i is a dogleg severity value; and D i is a depth along the planned wellpath at a survey location i; combining the expected wellbore energy with the expected tortuosity to create a first value indicative of probability of occurrence of a drilling event for drilling along the planned wellpath; and then drilling a wellbore along at least a first portion of the planned wellpath if the first value indicates the probability of the drilling event is less than a predetermined threshold probability.