Generaliserad β-VDR-baserad derivatberäkning för robust käntlektavgränsning och djupbestämning från potentifältdata

Att se dolda strukturer under våra fötter

Mineraler, grundvatten, geotermiska resurser och olja ligger ofta djupt under markytan, dolda för direkt observation. Geofysiker använder subtila variationer i jordens gravitations- och magnetfält för att kartlägga dessa begravda strukturer utan borrning. Men att översätta de känsliga mätningarna till tydliga bilder av sprickor, dagbrott och kontaktytor är komplicerat, eftersom vanliga bearbetningsmetoder tenderar att förstärka brus tillsammans med signalen. Den här studien presenterar ett sätt att skärpa underjordiska bilder samtidigt som bruset hålls i schack, vilket gör kartorna över undergrunden mer tillförlitliga för både vetenskap och prospektering.

Varför kanter under marken är viktiga

När bergarter med olika egenskaper eller täthet möts bildas ”kanter” i undergrunden—förkastningar, kontakter och intrång som ofta styr var vätskor, värme och malmer ackumuleras. Gravitions- och magnetundersökningar kan indirekt upptäcka dessa gränser som små variationer, så kallade potentialfältsanomalier. För att bestämma källornas lägen och djup beräknar tolkarna matematiska derivator av data, vilka framhäver var fältet ändras mest kraftigt. Tyvärr fungerar dessa derivator som filter som förstärker högfrekvent innehåll, så även små mängder slumpmässigt brus kan överväldiga de intressanta signalerna. Befintliga lösningar fungerar antingen bara när data är mycket rena eller kräver tung beräkning som är svår att tillämpa på dagens stora undersökningar.

Ett smartare sätt att ta skillnader

En tidigare metod, känd som β-VDR, erbjöd redan ett stabilare sätt att beräkna vertikala derivator genom att skickligt kombinera versioner av data som har uppåträckts—matematiskt projicerats till högre höjder för att släta ut brus. β-VDR ger renare vertikala derivator än standardfilter i Fourier-domänen, men hade två stora nackdelar. För det första beroende den fortfarande på mer sköra ändlig differens-formler för de horisontella komponenterna, vilket ledde till en obalans: vertikala derivator var robusta, horisontella var det inte. För det andra krävde den ursprungliga proceduren fem separata omgångar tunga Fourier-beräkningar, vilket gjorde den långsam och kostsam för stora rutnät.

Att balansera vertikal och sidovy

Författarna omformulerade β-VDR-idén till ett kompakt filter i frekvensdomänen som uppnår samma effekt med endast en framåt- och en invers Fouriertransform istället för fem. Detta minskar teoritiskt beräkningstiden med ungefär en faktor fem. De utvidgade sedan samma stabiliserande logik till horisontella derivator och skapade en matchande familj filter som de kallar β-HDR. Tillsammans bildar vertikala β-VDR och horisontella β-HDR ett enhetligt schema, β-VDR-with-β-HDR, som behandlar alla derivatriktningar konsekvent. I enkla termer slätar metoden bort brus precis tillräckligt i varje riktning samtidigt som den bevarar skarpa övergångar som markerar verkliga geologiska gränser.

Att testa metoden

För att kontrollera att den nya metoden är både korrekt och användbar genomförde teamet omfattande datorexperiment. De började med syntetiska modeller—idealiserade underjordiska block med kända former, djup och fysikaliska egenskaper—och genererade deras gravitations- och magnetiska svar. Genom att lägga till olika nivåer av slumpmässigt brus återskapade de den typ av röriga data som förekommer i verkliga undersökningar. Med hjälp av en standard teknik för kantförstärkning kallad totalgradienten, som beror på både vertikala och horisontella derivator, jämförde de fyra alternativ: traditionella Fourier-filter, en metod kallad ISVD, den ursprungliga β-VDR kombinerad med konventionella horisontella skillnader, och nya β-VDR-with-β-HDR. Den nya metoden reproducerade standardresultaten när ingen stabilisering användes, vilket bekräftade matematiken. Under brusiga förhållanden stack den tydligt ut: kanter förblev skarpa, falska toppar var sällsynta och uppskattade djup höll sig nära de verkliga värdena även när andra metoder havererade.

Från testmodeller till en verklig sedimentär bassäng

Författarna tillämpade sedan sin teknik på högupplösta aeromagnetiska data från den nigerianska delen av Chadbassängen, en region med tjocka sediment där förkastningar och intrång påverkar geotermisk och kolvätepotential. Utan att använda den vanliga förhandsutjämningen beräknade de de stabiliserade derivatorna och totalgradienten, och uppskattade därefter positioner och djup för magnetiska källor både i profil och i fulla 3D-vyer. Lösningarna överensstämde med kända regionala tendenser och avslöjade koherenta förkastningslika strukturer och intrång, inklusive grunda strukturer och djupare linjament som kan styra vätskeflöde. Viktigt är att djupuppskattningarna från 2D-profiler och 3D-rutnät överensstämde nära, vilket tyder på att resultaten inte är artefakter av metoden.

Tydligare bilder av underjorden för svåra data

För en icke-specialist är slutsatsen att detta arbete erbjuder ett bättre ”skärpningsfilter” för brusiga gravitions- och magnetkartor över undergrunden. Genom att omdesigna hur vertikala och horisontella skillnader beräknas extraherar β-VDR-with-β-HDR-metoden kanter och djup för begravda strukturer mer pålitligt, även när mätningarna är kontaminerade av betydande brus. Eftersom den också är mer effektiv att beräkna kan den tillämpas på stora moderna datamängder. Det betyder tydligare, mer trovärdiga bilder av vad som ligger under våra fötter—understödjande säkrare borrningsbeslut, förbättrade geotermiska bedömningar och en djupare förståelse av jordens dolda arkitektur.

Citering: Falade, S.C., Falade, A.H. Generalizing β-VDR-based derivative computation for robust source edge detection and depth estimation from potential field data. Sci Rep 16, 5672 (2026). https://doi.org/10.1038/s41598-026-36635-7

Nyckelord: gravitions- och magnetundersökningar, kantdetektion, djupbestämning, brusrobusta derivator, Chadbassängens geologi