Skip to main content

Blog

From a 1999 Patent to Transformers: 24 Years of Trying to Automate Fracture Picking

The problem statement for automated borehole-image interpretation was written down and patented in 1999, and it has not changed since. Schlumberger's US 5,960,371 already named the same three feature classes we target today (fractures, beddings, elliptical vugs) and read them as sinusoids in the unrolled image, descending straight from Hough's 1962 patent. What changed in a quarter century is not the target but the method that finally hit field-usable accuracy: a well-scale supervised transformer, in 2023.

Tarry SinghNarendra Patwardhanby Tarry Singh, Narendra Patwardhan8 min read
EarthScan insight

We spent part of a roughly twenty-month engagement with a mid-sized Middle East carbonate operator building a transformer that picks fractures off borehole-image logs. Somewhere in the middle of that work we pulled an old patent out of the reference folder and read it properly for the first time, and it stopped us. Schlumberger's US 5,960,371, filed on 4 September 1997 and granted on 28 September 1999, is titled "Method of determining dips and azimuths of fractures from borehole images." Thirteen claims. It describes a Hough-transform method that reads planar features off an unrolled borehole image, and it names the same three feature classes we were regressing in 2023: fractures, beddings, and elliptical vugs. The patent even cites its own ancestor, Paul Hough's original 1962 patent US 3,069,654 on detecting lines and curves in a parameter space. The problem we thought we were solving had been written down, claimed, and filed a quarter century before we opened our first log.

That is the whole argument of this piece. The problem statement for automated borehole-image interpretation did not change in twenty-four years. The three feature classes did not change. The geometry that turns each of them into a sinusoid did not change. What changed, and only recently, is the one thing that matters commercially: a method that clears field-usable accuracy.

A patent trail, not a fresh idea

Read the three documents in order and the lineage is exact. In 1962 Hough patents a way to find lines and curves by voting in a parameter space rather than in the image. In 1997 a Schlumberger team led by Saito files a patent that applies exactly that idea to borehole images: a tilted planar feature cuts the cylindrical borehole as an ellipse, the unrolled wall renders that ellipse as a sinusoid, and a Hough transform in the sinusoid's parameter space recovers its dip and azimuth. The 1999 grant makes the three targets explicit. Fractures and beddings are planar, so each is one sinusoid; an elliptical vug is a closed curve, handled in the same voting framework. The mechanics of why a planar feature is a sinusoid are the subject of our earlier primer on borehole image logs, so we will not re-derive them here beyond the fitting form the whole field agrees on:

y=Asin(0.0175x+φ)+offsety = A\,\sin(0.0175\,x + \varphi) + \text{offset}

where the amplitude encodes dip, the phase encodes azimuth, and the offset is depth. This is not our formulation and it is not new. It is the object the 1999 patent already parameterises. When we later built a classical baseline to compare against, it was a Canny edge detector at thresholds 200 and 255 with aperture 7, feeding a probabilistic Hough line transform at a vote threshold of 30 with a minimum line length of 30. That baseline is a direct descendant of the patented method. We did not choose it because it was clever; we chose it because it was the incumbent, the line every serious attempt at this problem had drawn since 1999.

Why the incumbent stalled

If the idea was patented in 1999, why was fracture picking still a manual job in 2021, done by an interpreter dragging sinusoids across a screen well by well? The Hough approach has a structural ceiling that no parameter sweep fixes. It votes for the single most-supported curve in a window, which means it degrades exactly where the geology gets interesting. In a densely fractured zone, three or four sinusoids overlap in the same interval; a vote-based method collapses them into one dominant curve or, worse, finds nothing at all. Our own archive holds a preserved failure case where the Hough baseline returned zero lines on a section a geologist reads as clearly fractured. That is not a tuning error. It is what a voting method does when the signal it wants to isolate is buried under other signals of the same shape.

The classical stack also has no notion of "not measured." Borehole-image pads cover only part of the wall, so the unrolled image arrives with structural gaps, and edge-and-vote pipelines happily draw lines across them. Every downstream fix (path-opening morphology to recover overlapping sinusoids, unsupervised Hough auto-scribbles to bootstrap labels) was an attempt to prop up the same 1999 method, and we have written those attempts up separately. They bought accuracy at the margins. None of them moved the interpretation from "assistive" to "field-usable," which is the only line a client cares about.

THE PROBLEM STATEMENT DID NOT CHANGE IN 24 YEARS85%2023 fracture sensitivity at 8 cmSame three feature classes, same sinusoids, from a 1962 patent onwardTHE UNCHANGING TARGET · EACH PLANAR FEATURE IS ONE SINUSOIDfracturesbeddingselliptical vugsIdentical geometry under every era: a tilted plane cuts the borehole as an ellipse, unrolled to a sine.1962ROOT PATENTUS 3,069,654Hough1999THE PROBLEM, STATEDUS 5,960,371Saito et al. (Schlumberger)2023FIELD-USABLEGeoBFDT / DETRsupervised transformerACCURACY LEDGER · ONLY 2023 CROSSES FIELD-USABLE050100field-usable band (80%, reading aid)method onlymethod patented85%at 8 cmdip 90% at 3 degazimuth 92% at 15 degDEPTH-MATCH THRESHOLDthe sensitivity claim is threshold-conditioned, not one headline number8 cm85%3 cm65%sourced: US 3,069,654 (1962); US 5,960,371 (1997/1999, 13 claims); GeoBFDT/DETR 2023, 14 wells
One problem statement, three eras. Above the spine, the target that never changed: three feature classes (fractures, beddings, elliptical vugs), each a single sinusoid in the unrolled borehole image, drawn identically under 1962, 1999 and 2023. The spine marks the patent trail: Hough's root patent US 3,069,654 (1962), Schlumberger's US 5,960,371 (filed 1997, granted 1999, 13 claims, targeting the same fractures, beddings and elliptical vugs via a Hough method), and the 2023 supervised transformer trained on 14 vertical wells. The accuracy ledger below carries the argument: a labelled field-usable band at 80 percent that only the 2023 result clears. The orange bar is the sole element that argues, the 2023 fracture sensitivity, and the toggle flips its depth-match threshold between 8 cm (about 85 percent) and 3 cm (about 65 percent) so the accuracy claim reads as threshold-conditioned rather than a single number. Dip accuracy of about 90 percent at 3 degrees and azimuth accuracy of about 92 percent at 15 degrees sit in their own lane. The patent numbers, claim count, well count and 2023 metrics are sourced; the 80 percent field-usable line is a labelled reading aid, not a measured threshold.

What actually changed in 2023

The break was not a new problem statement. It was a new method applied to the old one. We treated each sinusoid as an object and had a detection transformer regress its parameters directly, as a set, rather than painting a pixel mask and fitting curves afterward. The anchor-free, mask-free adaptation of object queries to sinusoid picking is its own subject, covered in a companion piece. The accuracy that switch bought is the number the patent trail was waiting on for twenty-four years.

Trained on fourteen vertical wells, the transformer reached roughly 85 percent fracture sensitivity at an 8 cm depth-matching threshold, with dip accuracy near 90 percent at 3 degrees and azimuth accuracy near 92 percent at 15 degrees. Tighten the depth threshold to 3 cm and fracture sensitivity falls to about 65 percent, worth stating plainly: the sensitivity claim is conditioned on how strict a depth match you demand, not a single headline. A 3 cm threshold is close to the raster floor of the logs, so a model that holds two thirds of its picks at that tolerance is operating near the physical limit of the data. At 8 cm, the tolerance an interpreter actually works to, the same model clears the bar a supervised interpreter would sign off on. The overlapping-fracture case that broke the Hough baseline is where set prediction earns its keep: a fixed set of learned queries can emit three sinusoids in the same window because none of them competes for a single vote.

The lesson the patent teaches

It is tempting to read this as "AI finally solved an old problem," but the patent trail argues something more specific. The problem was correctly specified in 1999. Saito and colleagues knew the target was three feature classes rendered as sinusoids, and they wrote a method to recover them. They were not wrong about the problem. They were limited by a class of algorithm that cannot separate overlapping instances of the shape it is voting for. The transformer did not out-think the 1999 patent about what to detect. It replaced a voting rule with a supervised set predictor, and that substitution is what carried the accuracy across the field-usable line.

For anyone deciding where to spend effort on subsurface automation, that is the reusable point. When a task has been "almost automated" for decades, the bottleneck is rarely the problem statement, which someone competent usually nailed years ago. It is the method's structural ceiling. The question to ask of a long-stalled interpretation task is not "what is the target," since a patent from the 1990s probably already answers that, but "what class of model does the target require, and has it existed until now." For borehole-image fracture picking, the answer arrived in 2023, twenty-four years after the problem was filed.

Limitations

The 2023 metrics come from models trained on fourteen vertical wells from a single carbonate operator; sensitivity at a 3 cm threshold sits near the raster resolution of the logs, so it should be read as a floor on that dataset rather than a portable benchmark, and generalisation to other basins, tool types, or highly deviated wells is not established by this work. The field-usable band shown in the instrument at 80 percent is a labelled reading aid to separate assistive from sign-off-grade accuracy, not a measured or contractual threshold. The patent history is drawn from the granted documents themselves; our reading of why the Hough approach stalled reflects our own baseline experiments on this dataset and the preserved failure cases in that archive, not an exhaustive survey of every classical borehole-image method published since 1999.

References

[1] Hough, P. V. C. "Method and means for recognizing complex patterns." US Patent 3,069,654, granted 18 December 1962. https://patents.google.com/patent/US3069654A

[2] Saito, T. et al. (Schlumberger Technology Corporation). "Method of determining dips and azimuths of fractures from borehole images." US Patent 5,960,371, filed 4 September 1997, granted 28 September 1999. https://patents.google.com/patent/US5960371A

[3] Carion, N. et al. "End-to-End Object Detection with Transformers (DETR)." ECCV 2020. https://arxiv.org/abs/2005.12872

Go to Top

© 2026 Copyright. Earthscan