Book a demo
  1. Home
  2. Articles
  3. Well engineering

Directional surveys · worked example · Llanos Basin, Colombia

Directional survey quality: magnetic corrections, sag and gross errors

DrillIQ drilling engineeringUpdated 1 October 202613 min read

Every position in a well file, every target hit and every anti-collision scan rests on the survey, and an MWD survey is only as good as its magnetic reference and its corrections. This article works through survey quality on a horizontal well near the equator: the corrections from magnetic to grid north, the field checks every station should pass, drill-string interference and sag, and the gross errors that no tolerance band catches.

Where survey error comes from

An MWD tool carries three accelerometers and three magnetometers. From the gravity vector it gets inclination and toolface; from the magnetic field, projected on the horizontal plane, it gets the azimuth from magnetic north. Each step adds error.

The reference. Magnetic north is turned to grid north with a declination from a model of the earth's field, which drifts year by year and does not see local magnetism in the crust.

The drill string. Steel above and below the non-magnetic collars magnetises along the string, so the axial magnetometer reads a field that is not the earth's. This axial bias is the largest single error in most MWD surveys.

The BHA. The sensor's collar bends under its own weight, so in a high-angle hole its axis points slightly below the hole's and it reads a lower inclination. This is sag.

People and files. A declination applied with the wrong sign, a mistyped inclination, a mixed-up north reference. Rare, and large when they happen.

The ISCWSA error models describe the first three as ellipses (see anti-collision: separation factor, ellipses and travelling cylinder). Survey QC checks that each real station behaves as the model assumes, and corrects what can be corrected.

From magnetic north to grid north

The tool measures an azimuth from magnetic north. Well plans and maps are in grid north. The conversion has two terms:

Agrid = Amag + D − γ

γ = atan( tan Δλ · sin φ )

D
magnetic declination, east positive, from a field model on the survey date
γ
grid convergence: the angle between grid north and true north at the well
Δλ, φ
longitude from the UTM zone's central meridian, and latitude

Near the equator the two terms behave very differently. Convergence scales with sin φ, so at 4.3° N it is small: +0.195° at the example well, in UTM zone 18. Declination there is large and moving. The World Magnetic Model 2025 gives −9.77° (west) on the example's survey date, changing by −0.18° a year. A reference computed when the well was planned, more than a year before it was drilled, is already a fifth of a degree out.

The field checks: gravity, total field and dip

Each MWD station measures the full gravity and magnetic vectors, so it can be checked against what the earth should give at that place and date. Three numbers are compared:

  • G total, the length of the gravity vector, against 1 g. Typical tolerance ±0.0025 g.
  • B total, the length of the magnetic vector, against the model's total field. Typical tolerance ±300 nT.
  • Dip, the angle of the magnetic vector below horizontal, against the model's dip. Typical tolerance ±0.45°.

A station outside a band is suspect: a sensor fault, steel nearby, or a reading taken while the string was moving. The trap is the reverse. A station inside every band can still carry a large azimuth error, because an axial bias changes the azimuth far more than it changes B total or dip in some directions. The azimuth error from an axial bias ΔBz is, to first order:

ΔA ≈ ΔBz · sin I · sin Am / ( B cos δ )

ΔBz
axial magnetic bias from the drill string (nT)
I, Am
inclination and magnetic azimuth of the station
B, δ
total field and dip; B cos δ is the horizontal field the compass works from

The error is largest in a horizontal hole running magnetic east or west, and smallest where the horizontal field is strong. That second factor is where the region comes in.

The example well

The example is a horizontal producer in a sandstone reservoir, set in typical conditions of the Llanos foreland of eastern Colombia: flat ground, a rotary table 260 m above sea level, the reservoir at about 2,340 m TVD and a 1,000 m lateral heading east-north-east. The top hole is gyro-surveyed; the 12¼″ and 8½″ sections are surveyed by MWD every stand.

The example well and its magnetic reference
Locationabout 4.3° N, UTM zone 18; grid convergence +0.195°
Field (WMM2025, survey date)total field 29,627 nT, dip 25.25°, declination −9.77°; horizontal field 26,795 nT
Trajectoryvertical to 1,850 m, build at 3.5°/30 m to land at 90° on 065° grid at 2,621 m MD, 2,341 m TVD; lateral to TD at 3,621 m MD
Hole and surveys17½″ to 600 m (gyro); 12¼″ to the landing, 9⅝″ casing there (MWD run 1); 8½″ lateral (MWD run 2); 102 MWD stations
What the survey carriesaxial drill-string bias 180 nT (run 1) and 320 nT (run 2); sag 0.10° and 0.20° at 90°, typical values from a BHA sag model; sensor noise 40 nT and 0.0003 g per axis; one station 19 m below the shoe disturbed by the casing
Magnetic heading of the lateral075°: 15° from magnetic east

The raw readings were built with the module's own forward model from a known hole, so every error can be measured against the truth, then run through the same QC and corrections as a real survey.

A hand check of the lateral

  1. Azimuth error per station. In the lateral, 320 nT × sin 90° × sin 75° / 26,795 nT = 0.0115 rad, or 0.66°. The engine measures a mean of −0.66° across the lateral's stations.
  2. Sideways drift of the lateral. 1,000 m × 0.0115 = 11.5 m. The build adds about 2.5 m: 180 nT gives 0.0065 rad at 90°, weighted by sin²I, which averages 0.5 over the 771 m build. Together about 14.0 m to the left.
  3. Sag. 0.20° is 0.00349 rad, so 1,000 m of lateral read 0.20° low puts TD 3.49 m too deep. The build adds 771 m × 0.5 × 0.00175 = 0.67 m. Together 4.16 m.
  4. An old reference. The declination 14 months earlier was −9.56°, 0.21° from the survey date's value. Over the well's 1,491 m horizontal departure that is 1,491 × 0.00365 = 5.4 m.

The engine, running minimum curvature on every station, gives 13.96 m left, 4.13 m deep and 5.47 m. These are not ellipse sizes. They are systematic errors: every station shares them, so they add up along the hole instead of averaging out.

What the station checks show

Figure 1 Total-field check, station by station
  • As measured
  • After multi-station analysis
  • Tolerance ±300 nT
-400-300-200-10001002003001,0001,5002,0002,5003,0003,500Survey depth MD (m)B total − reference (nT)kick-offlanding, 9⅝″ shoe-347 nT at 2,640 m:casing at the shoe,station rejected

DrillIQ's magnetic QC (corrections.station_qc) against WMM2025 at the well on the survey date: 29,627 nT, dip 25.25°. Raw readings made with the engine's forward model from a known hole, with a drill-string bias, sensor noise and sag.

Show the data
MDAs measuredAfter MSA
63070-14
66028-59
72020-70
78028-60
81085-2
87012138
93014664
9609816
1,02018-68
1,08040-50
1,140921
1,17015062
1,23017287
1,290886
1,32015674
1,38073-12
1,440912
1,470-52-143
1,53010514
1,59030-58
1,6509-75
1,68040-42
1,74015470
1,800924
1,830-15-106
1,89012835
1,95091-2
1,98012230
2,04086-7
2,10010911
2,16035-66
2,19058-42
2,25076-20
2,310892
2,34060-23
2,40012850
2,46013054
2,49010432
2,5509835
2,6108233
2,640-347-36
2,67012841
2,73013243
2,79012848
2,82070-2
2,8807511
2,9409324
2,97011539
3,03072-17
3,09069-27
3,15012639
3,18058-24
3,24012957
3,30065-10
3,3302-78
3,39059-32
3,45084-9
3,48010716
3,54061-16
3,60012053
3,62168-4

Two of 102 stations fail. The first lateral station, 19 m below the 9⅝″ shoe, reads B total 347 nT low: casing steel close to the sensor. Its azimuth is 3.4° off. It is rejected and the next station carries the survey on.

Every other lateral station passes, between +3 and +143 nT. Yet each of them is 0.66° off in azimuth. Heading 075° magnetic, the tool's axis sits almost square to the earth's field, so a bias along it barely changes the field's length. A total-field check cannot see an axial bias in an east-west lateral.

Figure 2 Dip check, station by station
  • As measured
  • After multi-station analysis
  • Tolerance ±0.45°
-0.45-0.30-0.150.000.150.300.451,0001,5002,0002,5003,0003,500Survey depth MD (m)Dip − reference (°)kick-offlanding, 9⅝″ shoe+0.46° at 1,350 m:axial bias invertical hole

Same stations and engine as the total-field chart. In the vertical hole the axial bias shows in the dip; in the lateral it hardly shows at all.

Show the data
MDAs measuredAfter MSA
6300.22-0.13
6600.23-0.12
7200.30-0.04
7800.29-0.06
8100.35-0.00
8700.14-0.21
9300.380.03
9600.23-0.13
1,0200.17-0.18
1,0800.29-0.05
1,1400.410.06
1,1700.18-0.17
1,2300.20-0.15
1,2900.24-0.12
1,3200.370.02
1,3800.25-0.10
1,4400.34-0.01
1,4700.28-0.06
1,5300.380.03
1,5900.360.01
1,6500.32-0.03
1,6800.23-0.13
1,7400.30-0.05
1,8000.21-0.14
1,8300.420.07
1,8900.380.04
1,9500.32-0.01
1,9800.330.00
2,0400.320.00
2,1000.340.05
2,1600.290.04
2,1900.13-0.11
2,2500.340.14
2,3100.15-0.02
2,3400.160.01
2,4000.180.06
2,4600.02-0.06
2,4900.110.05
2,550-0.02-0.03
2,610-0.09-0.05
2,640-0.05-0.01
2,670-0.17-0.11
2,730-0.030.05
2,790-0.19-0.09
2,820-0.050.05
2,880-0.020.06
2,9400.010.06
2,970-0.08-0.04
3,030-0.11-0.06
3,090-0.10-0.03
3,150-0.010.08
3,180-0.19-0.09
3,240-0.040.04
3,300-0.020.05
3,330-0.050.01
3,390-0.060.00
3,450-0.11-0.03
3,480-0.10-0.01
3,5400.020.11
3,6000.040.12
3,621-0.17-0.13

The dip chart shows the other side. In the vertical hole, 630–1,850 m, the dip reads 0.29° high on average and one station fails at +0.46°. There the tool's axis is vertical and the earth's field only 25° below horizontal, so an axial bias tilts the field and moves the dip. Near the equator a bias in vertical hole shows in the dip; at high latitudes it shows in B total. In the lateral the dip moves less than 0.3° and every station passes.

After multi-station analysis, below, both checks tighten around zero: the run-1 RMS falls from 103 nT and 0.23° to 37 nT and 0.07°, and run 2 from 86 nT and 0.12° to 39 nT and 0.08°.

Why the latitude matters

The same 300 nT axial bias, in a horizontal hole on the same 065° grid heading, gives different errors in different fields. Each column below is DrillIQ's engine evaluating the station in WMM2025 at that place.

A 300 nT axial bias at 90° inclination, 065° grid
Llanos, 4.3° NNorth Oman, 22.4° N60° N
Total field, dip29,627 nT, 25.3°43,258 nT, 34.9°51,283 nT, 72.7°
Horizontal field26,795 nT35,464 nT15,244 nT
Magnetic azimuth075.0°063.1°062.8°
Azimuth error0.62°0.43°0.99°
B total and dip change, lateral+72 nT, −0.07°+112 nT, −0.10°+42 nT, −0.15°
B total and dip change, vertical hole+129 nT, +0.52°+172 nT, +0.32°+287 nT, +0.10°

Low dip helps: most of the field is horizontal, which is what the compass uses. At 60° N only 15,244 nT is horizontal and the same bias costs 0.99°. But the field over this part of South America is weak overall, so the Llanos has less horizontal field than North Oman. And with a 9.8° west declination, a lateral that looks north-east on the map runs closer to magnetic east than its grid heading suggests. Check every lateral's magnetic azimuth.

Correcting the survey

Multi-station analysis (MSA). All the stations of one run share the same sensor and the same drill string, so they share the same biases. MSA finds the bias vector that brings every station's B total and dip closest to the reference, by least squares weighted by the tolerances, then recomputes each azimuth from the corrected field. With five or more stations at toolfaces spread over at least 90°, the cross-axial biases can be solved too; otherwise only the axial one. On the example, run 1's 25 stations give an axial bias of 200.6 nT, run 2's 34 stations 336.0 nT, against 180 and 320 nT put in. The cross-axial terms come out under 10 nT. The lateral's mean azimuth error falls from −0.66° to +0.03°.

Single-station axial correction. The alternative distrusts the axial magnetometer at each station and solves its value from the reference field alone. It needs no run of stations, but it fails where the axial reading cannot be told from the field: near horizontal and near magnetic east-west. Its conditioning is √(cos²I + sin²I cos²Am), and DrillIQ leaves stations below 0.2 uncorrected and flags them. This lateral averages 0.26: usable, but only just. On a lateral turned 10° further east, MSA with spread toolfaces would be the only magnetic correction left.

Sag. The sensor reads the inclination low by the sag angle times sin I. The angle comes from a BHA deflection model for the actual collar sizes, stabilisers and weight on bit; DrillIQ applies it per run. On the example the lateral's inclination error falls from −0.197° to +0.003°.

In-field referencing (IFR). A local survey of the earth's field replaces the global model's values with measured ones. It removes the crustal part of the declination error and tightens the B total and dip tolerances. For small targets or close offsets, IFR with MSA is the usual combination.

Figure 3 What each uncorrected error does to the position
  • Drill-string bias not corrected: sideways (+ right)
  • Sag not corrected: TVD (+ deeper)
  • Reference from 14 months earlier: sideways
  • Everything corrected: 3D distance
-15-10-505102,0002,4002,8003,2003,600Survey depth MD (m)Difference from the true hole (m)landing−14.0 m left+4.1 m deep+6.4 m right0.9 m

Minimum curvature (DrillIQ's survey engine) on the example's stations, each line with one correction left out and the others applied, against the known hole. The shoe station is rejected in every case.

Show the data
MDAxial bias left in (sideways)Sag left in (TVD)Reference 14 months old (sideways)All corrected (3D)
1,8000.090.000.120.09
1,8300.090.000.120.09
1,8600.100.000.130.10
1,8900.110.000.150.11
1,9200.120.000.170.12
1,9500.120.000.190.13
1,9800.110.000.230.14
2,0100.100.010.270.15
2,0400.070.020.310.15
2,0700.020.030.330.13
2,100-0.040.040.370.12
2,130-0.110.050.410.10
2,160-0.210.070.440.07
2,190-0.300.090.490.05
2,220-0.410.110.540.03
2,250-0.550.140.600.01
2,280-0.630.170.710.04
2,310-0.730.200.840.08
2,340-0.860.240.950.11
2,370-1.000.281.070.13
2,400-1.150.311.180.14
2,430-1.300.361.320.18
2,460-1.490.401.410.17
2,490-1.660.451.540.20
2,520-1.800.511.710.26
2,550-1.990.551.840.28
2,580-2.180.601.970.30
2,610-2.370.652.110.33
2,621-2.470.682.160.34
2,670-3.030.842.350.36
2,730-3.731.052.600.38
2,760-4.101.162.700.37
2,790-4.451.262.820.38
2,820-4.781.372.960.42
2,850-5.111.473.100.45
2,880-5.441.573.250.48
2,910-5.791.673.380.50
2,940-6.161.773.480.50
2,970-6.501.873.620.52
3,000-6.811.973.780.58
3,030-7.112.063.950.64
3,060-7.442.164.090.67
3,090-7.802.264.200.67
3,120-8.152.364.320.68
3,150-8.472.474.460.71
3,180-8.832.574.570.72
3,210-9.192.684.680.72
3,240-9.492.804.860.78
3,270-9.802.915.020.84
3,300-10.163.015.140.84
3,330-10.543.115.220.82
3,360-10.903.215.340.83
3,390-11.253.315.460.84
3,420-11.583.435.600.87
3,450-11.883.555.760.92
3,480-12.233.655.880.93
3,510-12.623.755.970.90
3,540-12.963.856.090.92
3,570-13.333.956.190.91
3,600-13.724.056.290.89
3,621-13.964.136.380.91

The chart shows each correction left out in turn, with the others applied. The drill-string bias is the largest: 14.0 m to the left at TD, growing from the landing, invisible to the total-field check. Sag puts the lateral 4.1 m deeper than it is, which matters in a thin reservoir and for the TVD of the next well's target. A reference 14 months old puts it 6.4 m right: the 5.5 m of the declination change plus the 0.9 m the fully corrected survey still carries from sensor noise and the residual bias.

Gross errors: the checks no tolerance band makes

Gross errors happen after the measurement, so the field checks never see them.

The declination's sign. Applied as +9.77° instead of −9.77°, every azimuth turns 19.5°. The example's TD moves 508 m. B total and dip are unchanged. DrillIQ recomputes each station's azimuth from its raw readings and the run's own reference and flags any that differ from the recorded value by more than 0.5°, so a wrong sign fails every station at once: all 59 above 3° in the example.

A mistyped station. An inclination of 89.56° keyed as 98.56° at 3,000 m reads as 8.86°/30 m into the station and 9.07°/30 m out of it. The dogleg check against the well's 6°/30 m limit flags both intervals. A single spike in DLS, up and back down, is almost always a typing error rather than the hole.

A north reference mixed up. When the positions in a file were computed from true azimuths and the azimuths are stored in grid, each horizontal step is turned by the convergence. The survey QC compares the stored steps with its own minimum-curvature steps and reports a constant +0.195° turn from 1,860 m, and says it matches the grid convergence. Near the equator that turn is small. At 60° N, 3° from the central meridian, the same formula gives 2.6°.

Independent checks. Compare the MWD with a gyro run where one is run, and the lateral's landing TVD with the log's marker depths. Two independent surveys disagreeing by more than their ellipses allow is the strongest gross-error test there is.

A checklist for the rig team

  • Compute the magnetic reference for the survey date, not the planning date; the declination here moves 0.18° a year.
  • Check the declination's sign and the north reference of every survey file on its first station.
  • QC every station's G total, B total and dip, and reject stations near casing shoes or fish.
  • Know each lateral's magnetic azimuth; near east-west, passing B total says little about azimuth.
  • Run MSA per run once it has five or more stations at spread toolfaces, and watch the run's RMS fall.
  • Apply a sag correction from a BHA model in every high-angle run.
  • Treat a DLS spike that goes up and straight back down as a typing error until shown otherwise.
  • Compare MWD with a gyro or with log markers wherever the target is tight.

How DrillIQ handles it

DrillIQ's Directional & Surveys module takes each run's magnetic reference from WMM2025 on the run's date, an IFR set or a typed value, and converts between magnetic, true and grid north. Every station with raw readings gets the gravity, total-field and dip checks against its tool's tolerances, and its recomputed angles are compared with those recorded. Corrections are kept per run beside the recorded values: sag, multi-station analysis or the single-station axial correction with its conditioning. Positions, reports and anti-collision scans all use the corrected survey, and imported files are checked for depth order, doglegs, gaps, position jumps and north-reference rotations. Everything above was worked out with the same engines.

Book a demo See all capabilities

References

  1. Williamson, H.S. Accuracy Prediction for Directional Measurement While Drilling. SPE 67616, SPE Drilling & Completion, 2000.
  2. Russell, A.W. and Roesler, R.F. Reduction of Nonmagnetic Drill Collar Length Through Magnetic Azimuth Correction Technique. SPE/IADC 13476, 1985.
  3. Ekseth, R. et al. The Reliability Problem Related to Directional Survey Data. IADC/SPE 103734, 2006.
  4. Wolff, C.J.M. and de Wardt, J.P. Borehole Position Uncertainty: Analysis of Measuring Methods and Derivation of Systematic Error Model. SPE 9223, Journal of Petroleum Technology, 1981.
  5. Chulliat, A. et al. The US/UK World Magnetic Model for 2025–2030: Technical Report.

About this example. The well in this article is a representative scenario set in typical Llanos Basin, Colombia conditions and worked through with DrillIQ's engines. It is not a client's well, and its figures are not field results.