Use eq 2.7, 1 - k^2 = Lsc / Loc. Since you have secondary referred leakage, use secondary referred magnetizing inductance (Loc) as well. Use coupled inductors, or Fig 3 or 4 (don't need to use separate leakage inductors if k ~= 1; then you put in the total as Lsc).
Is this just a two-winding transformer?
If you have more windings, then in general you need a triangular matrix* of k's for each set of windings, subject to certain limitations (the coupling between any two windings in a set of three, can't be less than the product of the other two coefficients).
*The matrix is symmetrical (k12 = k21, etc.) so we're only concerned about one triangle. At least, that's the usual case, as nonreciprocal transformers are very rare.
To illustrate that more clearly, suppose you create a transformer by wiring two independent transformers together. The coupling from primary 1 to secondary 1 (= primary 2) is k12, and the coupling from primary 2 to secondary 2 is k23. The coupling from primary 1 to secondary 2 is k13 = k12*k23, obviously enough**. You can't artificially have k13 less than this, but you can have more (say there were additional windings linking p1 to s2, but not to s1/p2).
**Unless it's not the straight product in k but a little finagling to get there, I'm not sure. In the k ~= 1 limit it should be (in which case k13 ~= 2 - k12 - k23 as well). Easy enough to prove what it actually is, anyway.
The full 2nd order transformer model also includes the DC resistance of the windings, the AC resistance of the core loss, the self-capacitance of each winding, and the isolation capacitance between windings. With this, you have a complete (if still approximate***) model of the efficiency, bandwidth and impedance of the transformer.
***This is the 2nd stage in an infinite series of circuits approximating a real component -- this is necessary because the actual EM fields within the transformer propagate at the speed of light, so can't be represented by mere lumped RLC elements. Except when they don't, which is relevant to skin effect in the wire and core, where the propagation velocity is very much slower than the speed of light. For these, the DCR and ACR components are modified (into RL networks). For the inter-turn and inter-winding fields, LC networks are used.
It's rarely worth modeling a transformer beyond 2nd order LC, plus say 3rd or 4th order losses. The reason is, you're mostly modeling the high frequency cutoff range, which is chock full of peaks and valleys, that are inconsistent between parts (depends on the exact number of turns per layer, etc.), and difficult or impossible to make use of in a practical circuit (the impedances, as well as the frequencies, are all over the place). So I just want to emphasize that, by taking it a step further with a couple resistors and capacitors, gets you to the all-around most useful wideband transformer model.
Tim