My understanding is that delta or N(d1) in the BSM is a decent approximation for probability of expiring in the money. N(d2) is what you really want. This guy gives a decent explanation:
"N is just the notation to say that we are calculating the probability under normal distribution.
D2 is the probability that the option will expire in the money i.e. spot above strike for a call. N(D2) gives the expected value (i.e. probability adjusted value) of having to pay out the strike price for a call.
D1 (delta) is a conditional probability. A gain for the call buyer occurs on two factors occurring at maturity. One, the spot has to be above strike price. (Direction). Two, the difference between spot and strike prices at maturity (Quantum). Imagine, a call at strike price $100. If the spot price of the stock is $101 or $150, the first condition is satisfied. The second condition is about whether the gain is $1 or $50. The term D1 combines these two into a conditional probability that if the spot at maturity is above, what will be its expected value in relation to current spot price."
Formulas:
N(d1) and
N(d2) are probability factors: N is the cumulative standard normal distribution function,
d2 = −log(X/S) − (r − 1 2σ2)τ σ √τ
d1 = d2 + σ √τ