h，上延高度分別為100、200、300、400、500米。
[0106] 對每一尺度的網(wǎng)格化TMI異常數(shù)據(jù)Th進行TDR計算；
[0107] 分別針對每一個上延高度，對計算得到的TDR，進行基于水平梯度的多方向邊緣檢 巧。。針對每個TDR計算結(jié)果進行64個方向的邊緣檢測，各方向a取值分別為kn/(2w)， 其中k= 0, 1, 2,…，（2"-1)，n= 6,得到各尺度的邊緣。
[0108] 隨后，采用骨骼算法對計算得到的各尺度的邊緣進行細化處理，得到各個尺度的 磁構(gòu)造格架圖。
[0109] 圖2為將網(wǎng)格化?I異常數(shù)據(jù)上延300米上延高度，邊緣經(jīng)骨骼算法細化提取的 相應一個尺度的單像素寬度磁構(gòu)造格架?？蒞看出，根據(jù)本發(fā)明的自動提取得到的單像素 點寬度的構(gòu)造格架更接近于實際地質(zhì)填圖識別構(gòu)造特點，便于進行地質(zhì)解釋。另外圖面信 息更清晰，便于不同尺度構(gòu)造格架的疊加分析。
[0110] 圖3為將網(wǎng)格化TMI異常數(shù)據(jù)分別向上延拓100、200、300、400、500米上延高度 后經(jīng)64個方向邊緣檢測得到的各尺度的構(gòu)造格架進行顏色漸變疊加形成的綜合構(gòu)造格架 圖。將上延后的各個高度對應于相應的源區(qū)深度，可將該綜合構(gòu)造格架圖用于表征研究區(qū) 域不同深度的構(gòu)造格架信息。用從灰白色到黑色漸變色代表了從低到高的向上延拓高度或 從淺入深的切割深度，綜合構(gòu)造格架圖反映了不同切割深度的構(gòu)造信息。
[0111] 隨后，計算各尺度格架上各邊緣點處的TDRh的梯度的模W表征構(gòu)造埋藏深度。圖 4為上延300米上延高度格架的邊緣點處的梯度的模所表征的構(gòu)造埋藏深度，從灰白到黑 色漸變顏色反應了水平梯度的模所代表的構(gòu)造埋藏深度逐漸減小。該圖反映了相應一個尺 度下構(gòu)造埋藏深度變化。
[0112] 隨后，將所述多個尺度的表征相對埋藏深度的構(gòu)造格架圖疊置生成綜合埋藏深度 構(gòu)造格架圖。圖5為將網(wǎng)格化的TMI異常數(shù)據(jù)分別上延100、200、300、400、500米上延高 度后計算得到的表征相對埋藏深度構(gòu)造格架圖進行疊加得到的綜合埋藏深度構(gòu)造格架圖。 從圖中可W看到，在不同尺度上，反映出主要構(gòu)造帶埋藏深度變化較大，通常埋藏深的構(gòu)造 即上覆地層較厚，對應的?I異常幅值較低，為隱伏構(gòu)造，因此該方法有助于隱伏構(gòu)造的識 別。
[0113] 隨后，對每一尺度的網(wǎng)格化?I異常數(shù)據(jù)Th計算S維解析信號ASh。計算得到的 解析信號ASh值與磁傾角值大小無關，可W用來指示磁異常源區(qū)位置和磁場強度。
[0114] 該實例中，各尺度格架上各邊緣點處ASh值用于指示構(gòu)造磁性強弱。圖6為上延 300米上延高度格架的邊緣點處采用ASh值指示磁性強弱的構(gòu)造格架圖，從灰白色到黑色漸 變顏色表征了解析信號ASh值所指示的構(gòu)造帶磁性逐漸增強。
[0115] 隨后，將所述多個尺度的磁性強弱構(gòu)造格架圖疊置生成綜合磁性強弱構(gòu)造格架 圖，如圖7所示，各尺度或深度邊緣的ASh值采用灰白色到黑色漸變顏色的疊加顯示，突出 顯示不同深度范圍的構(gòu)造格架的磁性變化信息。
[0116] 在不考慮源區(qū)埋藏深度的情況下，=維解析信號ASh的幅值大小指示了異常源的 磁性強弱，磁性強的構(gòu)造通常與成礦關系密切。
[0117] 另外，同時觀察圖3和圖5,可W看到具有切割深度和埋藏深度信息的構(gòu)造格架圖 能夠清楚的反映出區(qū)域中深度大延展長的主干構(gòu)造、深度較淺延展相對短的次級構(gòu)造、各 構(gòu)造的覆蓋厚度W及互相的交割關系，因此根據(jù)本發(fā)明的方法得到的構(gòu)造格架圖可W幫助 本領域技術人員認識所探測的區(qū)域構(gòu)造格局。
[0118] 沿磁性較強、深度較大的構(gòu)造帶附近、不同方向構(gòu)造的交匯部位W及構(gòu)造帶轉(zhuǎn)折 彎曲部位是發(fā)現(xiàn)潛在金屬礦床的重要位置，根據(jù)本發(fā)明的方法利用磁測數(shù)據(jù)可W快速準確 提取得到具有切割深度、埋藏深度、磁性強弱、主次關系和交割關系等信息的構(gòu)造格架圖， 幫助勘探人員準確快速地發(fā)現(xiàn)潛在的金屬礦床。
[0119] W上借助優(yōu)選實施例對本發(fā)明進行了詳細說明，但是本發(fā)明不限于此。本發(fā)明提 出的低紳度地區(qū)磁構(gòu)造格架自動提取方法同樣適用于中高紳度地區(qū)的磁構(gòu)造格架的自動 提取，另外同樣適用于重力位場構(gòu)造格架的自動提取。
[0120] 本技術領域技術人員可W根據(jù)本發(fā)明的原理進行各種修改。因此，凡按照本發(fā)明 原理所作的修改，都應當理解為落入本發(fā)明的保護范圍。
【主權(quán)項】
1. 一種構(gòu)造格架自動提取方法，該方法包括以下步驟： 對來自待研宄區(qū)域的磁法或重力測量數(shù)據(jù)進行預處理，得到總磁場強度（TMI)異常數(shù) 據(jù)或布格重力異常數(shù)據(jù)； 將所得到的TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)網(wǎng)格化，并將網(wǎng)格化的TMI異常數(shù)據(jù)或 布格重力異常數(shù)據(jù)向上延拓多個預定高度，得到多個尺度的網(wǎng)格化TMI異常數(shù)據(jù)或布格重 力異常數(shù)據(jù)Th，h為向上延拓后的高度； 分別利用每一尺度的網(wǎng)格化TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)Th計算各尺度的TMI異 常數(shù)據(jù)或布格重力異常數(shù)據(jù)的傾斜導數(shù)TDRh; 分別針對每一尺度的網(wǎng)格化TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)的傾斜導數(shù)，進行基于 水平梯度的多方向邊緣檢測，得到各尺度的磁或重力異常源體邊緣； 采用形態(tài)學骨骼算法將計算得到的各尺度的磁或重力異常源體邊緣分別細化為單像 素寬度，得到多個尺度的構(gòu)造格架圖。2. 如權(quán)利要求1所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括，將計 算得到的所述多個尺度的構(gòu)造格架圖疊置生成綜合構(gòu)造格架圖。3. 如權(quán)利要求2所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括，將網(wǎng) 格化TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)向上延拓多個預定高度后提取的邊緣對應于不同 深度的構(gòu)造，對得到的各深度的構(gòu)造格架圖進行疊加得到反映不同切割深度信息的綜合構(gòu) 造格架圖。4. 如權(quán)利要求1所述的構(gòu)造格架自動提取方法，其特征在于，所述分別針對每一尺度 的網(wǎng)格化TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)的傾斜導數(shù)進行基于水平梯度的多方向邊緣 檢測，包括以下步驟：該水平梯度的輻角為：貝1J，對于高度h方向a的磁或重力異常源體邊緣點為模沿輻角方向 有局部極大值的點； 針對每一方向a，將傾斜導數(shù)IDRh的水平梯度的模的局部極大值點沿梯度的垂直方向 連接，得到的曲線構(gòu)成邊緣； 針對同一高度h，以多個不同的方向a計算邊緣，對計算得到的各邊緣求并集得到相 應尺度的磁或重力異常源體邊緣， 其中，所述多個不同方向a取值為分別1^3!八21^)，1^ = 0，1，2，"％(211-1)，11為大于或 等于2的整數(shù)。5. 如權(quán)利要求4所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括： 分別將各尺度構(gòu)造格架圖上的每一邊緣點處的水平梯度的模表征該尺度構(gòu)造格架圖 中該邊緣點處的構(gòu)造埋藏深度，得到多個尺度的表征構(gòu)造埋藏深度的構(gòu)造格架圖。6. 如權(quán)利要求5所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括將所 述多個尺度的表征構(gòu)造埋藏深度的構(gòu)造格架圖疊置生成綜合埋藏深度構(gòu)造格架圖。7. 如權(quán)利要求1所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括， 分別基于每一尺度的網(wǎng)格化TMI異常數(shù)據(jù)或布格重力異常數(shù)據(jù)Th計算三維解析信號 ASh，得到各邊緣點的ASh值，由此得到多個尺度的表征邊緣點處磁性或密度強弱的構(gòu)造格 架圖。8. 如權(quán)利要求7所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括將所 述多個尺度的表征邊緣點磁性或密度強弱的構(gòu)造格架圖疊置生成綜合磁性強弱構(gòu)造格架 圖或綜合密度強弱構(gòu)造格架圖。9. 如權(quán)利要求1所述的構(gòu)造格架自動提取方法，其特征在于，該方法進一步包括，在進 行邊緣檢測前對計算得到的傾斜導數(shù)進行去除噪聲的處理。10. 如權(quán)利要求1所述的構(gòu)造格架自動提取方法，其特征在于，該方法適用于低煒度地 區(qū)的磁構(gòu)造格架自動提取；優(yōu)選地，該方法適用于磁傾角為±30°之間的區(qū)域的磁法測量 數(shù)據(jù)；進一步優(yōu)選地，該方法適用于磁傾角為±20°之間的區(qū)域的磁法測量數(shù)據(jù)。
【專利摘要】本發(fā)明提供一種構(gòu)造格架自動提取方法。該方法包括以下步驟：對網(wǎng)格化磁或重力異常數(shù)據(jù)進行多個尺度下基于傾斜導數(shù)的水平梯度的多方向邊緣檢測，分別得到各尺度的邊緣；采用形態(tài)學骨骼算法將計算得到的各尺度的邊緣細化；通過分別賦予每個邊緣點切割深度、埋藏深度和磁性或密度強弱屬性，得到多個尺度的各種類型的構(gòu)造格架圖，并通過將相同類型的各尺度構(gòu)造格架圖疊置生成各類綜合構(gòu)造格架圖。根據(jù)本發(fā)明的方法，可以實現(xiàn)對控制礦床形成的地質(zhì)構(gòu)造的識別和定性解釋，并實現(xiàn)金屬礦床的靶區(qū)定位。本發(fā)明特別適用于低緯度地區(qū)的磁構(gòu)造格架的自動提取。
【IPC分類】G01V3/38
【公開號】CN104965232
【申請?zhí)枴緾N201510303316
【發(fā)明人】曹殿華
【申請人】中國地質(zhì)科學院礦產(chǎn)資源研究所
【公開日】2015年10月7日
【申請日】2015年6月4日