国产自产21区,亚洲97,免费毛片网,国产啪视频,青青青国产在线观看,国产毛片一区二区三区精品

電性源瞬變電磁地空探測方法
【專利摘要】本發明涉及一種電性源瞬變電磁地空探測方法，包括：采用接地長導線作為發射源，利用直升機或者無人機搭載接收線圈進行數據采集；采用全域視電阻率方法對觀測數據進行處理，完成對深部地質目標的定性解釋，獲得地質目標的產狀、走向、規模等概況信息；采用瞬變電磁波場變換方法、三維偏移成像技術及逆逆合成孔徑方法對地空數據進行解釋處理，完成對深部目標體的精細探測，獲得深部地質結構的深度、形狀、電性結構等詳細信息。本發明能夠實現快速、大規模的電磁探測，能夠對深部地質目標體進行詳細勘查，能夠獲得深層地質目標體的準確、詳細的地質信息。
【專利說明】電性源瞬變電磁地空探測方法

【技術領域】
[0001] 本發明涉及地球物理勘探領域，特別是涉及了瞬變電磁探測方法與解釋技術。

【背景技術】
[0002] 隨著國民經濟的飛速發展，礦產資源消耗和需求進入了持續高增長階段，供需矛 盾日趨嚴重。東部平原地區利于勘探，開發時間較早，由于長期的開發，資源枯竭的問題已 經顯現，礦產資源供給不足已成為制約我國經濟發展的瓶頸。為了保持國民經濟的快速發 展，迫切需要尋找礦產資源供給的接替區和資源戰略儲備基地。首先，青藏高原橫跨全球三 大成礦域中的古亞洲和特提斯-喜馬拉雅兩大成礦域。尤其是特提斯-喜馬拉雅成礦域自 中生帶以來，經歷了多島弧碰撞造山演化歷程。在喜馬拉雅隆升階段，受印度板塊向北持續 俯沖的影響，在青藏高原形成了沿雅魯藏布江一帶分布的R底斯火山-巖楽巖帶，伴隨著 大規模的成礦作用，是最有找礦潛力的地區。其次，我國擁有漫長的海岸線和廣闊的淺海， 淺海礦產資源主要是石油與天然氣和各類濱海砂礦，最近還發現一種極富發展前景的天然 氣水合物可燃冰。青藏高原地區和海洋大陸架由于其豐富的礦產資源優勢，必將成為我國 礦產資源供給的接替區和資源戰略儲備基地。面對青藏高原地形復雜、高差懸殊、氣象多變 等特殊的自然條件，目前僅完成了大比例的基礎地質調查，詳細的地質勘查還屬空白。海洋 條件惡劣，在海上作業，海面的風浪襲擊，并在工作船航行過程中還有遇到暗礁的危險，在 這樣惡劣的環境中，船體在海水中不停地顛簸起伏，儀器裝置的海底投放和回收都極其困 難，和陸地想象的完全不同，甚至甲板上實驗人員的安全也受到威脅，實驗過程如果沒有很 好地計劃和組織，很難順利完成。為了克服上述種種困難，迫切需要我們提出合適的工作方 式以及新的解釋方法以適應這些地區的詳細地質勘探。
[0003] 航空電磁法（AirborneElectromagneticMethod簡稱AEM)，是航空物探常用的 測量方法之一。AEM具有速度快、成本低、通行性好、可大面積覆蓋、可用于海域等優勢，尤 其是在運積層或植被發育的覆蓋地區，它更具有一般勘探手段不可比擬的優勢，因此可以 用于高原海洋等地的勘察。但是航空瞬變電磁發射源受飛行器載重限制，使得航空瞬變電 磁發射磁矩有限，因此航空瞬變電磁的探測深度受限。目前應用較多的是地面中心回線 方式瞬變電磁法裝置，由于這種裝置簡單，易于解釋，得到廣泛的應用。當邊長較大時（大 于300米），為了提高工作效率可以采用大定源裝置進行測量，但是大定源解釋方法較中心 回線復雜，目前已經提出了一些大定源瞬變電磁的數據解釋方法。目前使用的儀器如V-5、 V8、⑶P-32、ROTEM、SIROTEM、PEM等大多采用地面大回線裝置形式。電性源瞬變電磁法利 用接地點電極通以脈沖電流而在地下建立起一次脈沖磁場。最早于60年代初發展起來的 是該方法的遠區工作方式，即長偏移距離法LOTEM(LongOffsetTransientEM)，俄羅斯稱 其為固定源建場法，工作時間屬于波區，即收發距離大于4倍勘探目的層的深度，實現探測 nX102-nX103m的深層勘探。瞬變曲線響應簡單，與頻率域方法基本上相似。時間域地空電 磁法（timedomaingroundairborneelectromagneticmethod)是融合地面TEM與航空 TEM優勢的一種新型勘探方法，通常將發射系統放置于地面，并鋪設幾公里的長接地導線或 大定源回線源，然后將接收系統、傳感器安裝在直升機、無人機或飛艇上進行飛行測量.此 方法不僅具有航空TEM的空間分辨率高、野外布線方便快捷、探測高效等優勢，還具有地面 TEM的大發射磁矩、信噪比高、勘探深度大的優勢，是深部礦產資源快速勘查的重要技術手 段·
[0004] 最早的半航空系統出現在上世紀70年代初期，名為TURAIR系統。它屬于頻率域 電磁系統，采用兩個分開的接收機確定振幅比和相位差。進入90年代，半航空瞬時電磁系 統FLAIRTEM系統和TerraAir系統先后問世。1997年12月，Fugro公司用TerraAir系統 進行了實驗，將其與航空TEM系統（GEOTEM)和地面TEM系統（PR0TEM37)進行實測對比。
[0005] 國外從上世紀90年代初開始研究地空電磁勘探技術，Tohru，將基于直升機的地 空電磁探測系統成功應用于mount bandai火山結構勘查；Hisatoshi等成功應用直升機地 空電磁系統探測到了 800m深處的地下水資源.我國地空電磁探測研究起步較晚，吉林大學 的研究人員于2009年開始研究地空電磁探測方法，由于直升機勘探費用高、飛行員及飛機 在飛行過程中存在風險，所以，研究小組研發了基于無人飛艇的地空電磁勘探系統.基于 飛艇的地空電磁勘探方法不僅可以解決直升機航空飛行勘探困難的問題，而且適用于地形 復雜的山區資源探測，在我國深部礦產資源探測和地質普查應用中具有廣闊的發展前景.


【發明內容】

[0006] 本發明所要解決的技術問題是提供電性源瞬變電磁地空探測方法及解釋技術，能 夠對深部地質目標體進行探測,并獲得深層地質目標體的詳細地質信息。
[0007] 為了解決上述問題，本發明公開了長導線源瞬變電磁地空探測方法，包括：
[0008] 采用長導線源瞬變電磁地空探測方法，對深部地質目標體進行探測，獲得瞬變電 磁的觀測數據；
[0009] 采用多點數據合成的方法對瞬變電磁的觀測數據進行處理及解釋，完成對深部地 質目標體的精細探測，獲得深部地質目標體的信息；
[0010] 采用瞬變電磁虛擬波場連續速度分析及成像方法，對瞬變電磁的觀測數據進行處 理及解釋，完成對深部地質目標體的精細探測，獲得深部地質目標體的信息。
[0011] 優選的，所述深部地質目標體的信息為深部地質目標體的位置信息、大小信息或 形狀信息。
[0012] 優選的，所述采用長導線源瞬變電磁地空探測方法，對深部地質目標體進行探測， 獲得瞬變電磁的觀測數據的步驟，包括：
[0013] 通過地面長導線源AB向空中發射電磁場；
[0014] 在高度范圍為50-150米的空中設置觀測點，對深部地質目標體進行探測，采集數 據，獲得瞬變電磁的觀測數據；
[0015] 其中，所述對深部地質目標體進行探測，獲得瞬變電磁的觀測數據的過程是通過 長導線源地空探測的工作裝置來完成的。
[0016] 優選的，所述地面長導線源AB的長度為3-5千米；所述地面長導線源AB的發射功 率不小于100千瓦。
[0017] 優選的，所述采用多點數據合成的方法對瞬變電磁的觀測數據進行處理及解釋， 完成對深部地質目標體的精細探測，獲得深部地質目標體的信息的步驟，包括：
[0018] 對瞬變電磁的觀測數據，采用逆合成孔徑雷達算法的方式，進行相關加權疊加，獲 得瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體；
[0019] 對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行克?；舴蚱瞥上瘢@ 得深層地質目標體的數字圖像；
[0020] 對所述深層地質目標體的數字圖像，進行數字圖像的分析，獲得深部地質目標體 的信息。
[0021] 優選的，所述對瞬變電磁的觀測數據，采用逆合成孔徑雷達算法的方式，進行相關 疊加，獲得瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體的步驟，包括：
[0022] 設計半空間中賦存高阻塊狀異常體的三維模型，獲得半空間中賦存高阻塊狀異常 體的三維模型；
[0023] 設計半空間中賦存低阻塊狀異常體的三維模型，獲得半空間中賦存低阻塊狀異常 體的三維模型；
[0024] 依據半空間中賦存低阻塊狀異常體的三維模型及半空間中賦存低阻塊狀異常體 的三維模型，采用逆合成孔徑雷達算法的方式，對瞬變電磁的觀測數據進行波場變換，獲得 瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體。
[0025] 優選的，所述對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行克?；舴?偏移成像，獲得深層地質目標體的數字圖像的步驟，包括：
[0026] 采用三維邊界元方式，對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行 克希霍夫偏移成像，獲得深層地質目標體的數字圖像。
[0027] 優選的，所述采用瞬變電磁虛擬波場連續速度分析及成像方法，對瞬變電磁的觀 測數據進行處理及解釋，完成對深部地質目標體的精細探測，獲得深部地質目標體的信息 的步驟，包括：
[0028] 采用數據插值的方式，獲得瞬變電磁的觀測數據的虛擬速度數據體的數據量；
[0029] 采用速度建模的方式，獲得瞬變電磁的觀測數據的連續速度分析體成像圖像；
[0030] 通過上述瞬變電磁虛擬波場連續速度分析及成像方法，對瞬變電磁的觀測數據進 行處理及解釋，完成對深部地質目標體的精細探測，獲得深部地質目標體的信息。
[0031] 優選的，所述采用數據插值的方式，獲得瞬變電磁的觀測數據的虛擬速度數據體 的數據量的步驟，包括：
[0032] 采用全局距離加權插值的方式，對散亂數據進行分析，獲得瞬變電磁的觀測數據 的最初虛擬速度數據體的數據量；
[0033] 采用近點線性插值的方式，對領域對角線的數據進行分析，獲得瞬變電磁的觀測 數據的精細虛擬速度數據體的數據量。
[0034] 優選的，所述采用速度建模的方式，獲得瞬變電磁的觀測數據的連續速度分析體 成像圖像的步驟，包括：
[0035] 采用射線方法，對瞬變電磁的觀測數據進行速度建模，獲得瞬變電磁的觀測數據 的連續速度分析體成像圖像；
[0036] 采用波動分析方法，對瞬變電磁的觀測數據進行速度建模，獲得瞬變電磁觀測數 據的連續速度分析體成像圖像；
[0037] 其中，所述波動分析方法為深度聚焦分析方法及剩余曲率分析方法。
[0038] 與現有技術相比，本發明具有以下優點：
[0039] 首先，本發明提供了長導線源瞬變電磁地空探測方法，是一種新型的瞬變電磁測 深法，能夠對深層地質目標體進行準確勘探，從而獲得深層地質目標體的準確、詳細的地質 信息。
[0040] 其次,本發明提供的長導線源瞬變電磁地空探測方法,相對于地面瞬變電磁系統 而言，具有方便、高效的優勢。尤其是在測量條件較為復雜的地區，如地勢起伏的山區，方 便、高效的優勢更加明顯。
[0041] 再者，本發明提供的長導線源瞬變電磁地空探測方法，相對于航空瞬變電磁系統 而言，具有信噪比更高、空間分辨率更好的優勢。
[0042] 另外，本發明提供的長導線源瞬變電磁地空探測方法，隨著導體埋藏加深，地面 TEM系統的晚期信噪比優勢將減弱，而本發明所提供的方法，可以克服隨著導體埋藏加深， 地面TEM系統的晚期信噪比優勢將減弱的缺陷。
[0043] 總之，本發明提供了長導線源瞬變電磁地空探測方法，是一種新型的瞬變電磁測 深法。

【專利附圖】

【附圖說明】
[0044] 圖1是本發明長導線源瞬變電磁地空探測方法實施例的流程圖；
[0045] 圖2是本發明實施例中瞬變電磁法原理示意圖；
[0046] 圖3是本發明實施例中瞬變電磁法工作機理示意圖；
[0047] 圖4是本發明實施例中長導線源地空探測方法工作原理示意圖；
[0048] 圖5是本發明實施例中長導線源地空探測方法工作裝置示意圖；
[0049] 圖6是本發明實施例中地井及井中瞬變電磁法的工作裝置示意圖；
[0050] 圖7是本發明實施例中航空瞬變電磁法的工作裝置示意圖；
[0051] 圖8是本發明實施例中瞬變電磁煙圈示意圖；
[0052] 圖9是本發明實施例中觀測曲線額波場轉換示意圖；
[0053] 圖10是本發明實施例中虛擬子波斷面對比示意圖；
[0054] 圖11是本發明實施例中子波寬度壓縮前后對比示意圖；
[0055] 圖12是本發明實施例中逆合成孔徑示意圖；
[0056] 圖13是本發明實施例中模型三維示意圖；
[0057] 圖14是本發明實施例中模型計算結果示意圖；
[0058] 圖15是本發明實施例中速度模型迭代算法的示意圖；
[0059] 圖16是本發明實施例中三維波場速度成像效果示意圖。

【具體實施方式】
[0060] 為使本發明的上述目的、特征和優點能夠更加明顯易懂，下面結合附圖和具體實 施方式對本發明作進一步詳細的說明。
[0061] 本發明的核心構思之一在于，發明提供了長導線源瞬變電磁地空探測方法，是一 種新型的瞬變電磁測深法，能夠對深層地質目標體進行準確勘探，從而獲得深層地質目標 體的準確、詳細的地質信息。
[0062] 參照圖1，示出了本發明長導線源瞬變電磁地空探測方法實施例的流程圖，具體可 以包括：
[0063] 步驟101、采用長導線源瞬變電磁地空探測方法，對深部地質目標體進行探測，獲 得瞬變電磁的觀測數據。
[0064] 本發明采用長導線源地空探測方法工作裝置形式，對深部地質目標體進行探測， 獲得觀測數據。
[0065]參照圖2，示出了本發明實施例中瞬變電磁法原理示意圖。
[0066] 瞬變電磁場法（TransientElectromagneticField，簡稱TEM)是一種建立在電磁 感應原理基礎上的時間域人工源電磁探測方法。它是利用階躍形波電磁脈沖激發，利用不 接地回線向地下發射一次場，在一次場斷電后，測量由地下介質產生的感應二次場隨時間 的變化，來達到尋找各種地質目標的一種地球物理勘探方法。
[0067]從圖2中，可以看出：電磁場感應渦流程產生的過程及發射信號與接收信號的關 系。
[0068]參照圖3，示出了本發明實施例中瞬變電磁法工作機理示意圖。
[0069] 瞬變電磁法是利用接地電極或者不接地回線通以脈沖電流，在地下建立起一次脈 沖磁場，在一次磁場間隙期間，利用探測線圈觀測二次渦流場，根據觀測信號來判斷地下介 質電性變化情況。
[0070] 從圖3中，可以看出：
[0071] 圖3A為發送方波電流信號的示意圖，圖3A可以清晰地說明發送方波電流信號情 況；
[0072] 圖3B為發送電流在大地中建立的磁場，即一次場的示意圖，圖3B可以清晰地說明 發送電流在大地中建立的磁場，即一次場的電流信號情況。
[0073] 圖3C為一次場消失后，接收線圈的自感信號的示意圖，圖3C可以清晰地說明一次 場消失后，接收線圈的自感信號情況。
[0074] 瞬變電磁場法的研究工作主要包括地面探測瞬變電磁法、地井及井中瞬變電磁法 及航空瞬變電磁法三個代表性裝置形式。
[0075] 所述對深部地質目標體進行探測，獲得瞬變電磁的觀測數據的過程是通過長導線 源地空探測的工作裝置來完成的。
[0076] 所述深部地質目標體的信息為深部地質目標體的位置信息、大小信息或形狀信 肩、。
[0077] 在本發明的一種優選實施例中，所述步驟101，具體可以包括：
[0078] 子步驟111、通過地面長導線源AB向空中發射電磁場；
[0079] 子步驟121、在高度范圍為50-80米的空中設置觀測點，對深部地質目標體進行探 測，采集數據，獲得瞬變電磁的觀測數據；
[0080] 其中，所述對深部地質目標體進行探測，獲得瞬變電磁的觀測數據的過程是通過 長導線源地空探測的工作裝置來完成的。
[0081] 參照圖4,示出了本發明實施例中長導線源地空探測方法工作原理示意圖。
[0082] 從圖4中，可以看出：本發明長導線源地空探測方法的工作原理。
[0083] 參照圖5,示出了本發明實施例中長導線源地空探測方法工作裝置示意圖。
[0084] 從圖5中，可以看出：本發明長導線源地空探測方法的工作裝置設計與以往不同。
[0085] 所述地面長導線源AB的長度為3-5千米；
[0086] 所述地面長導線源AB的發射功率為30千瓦。
[0087] 對本發明長導線源地空探測方法工作裝置做說明如下：
[0088] (I)AB為地面長導線源，測量在空中。
[0089] ⑵航線與AB平行。
[0090] (3)發射AB的長度3-5KM，發射功率30KW-50KW。
[0091] (4)觀測點在空中；可以是無人機、直升機等搭載接收裝置進行測量；飛機高度50 米-100米。
[0092] 參照圖6,示出了本發明實施例中地井及井中瞬變電磁法的工作裝置示意圖。
[0093] 從圖6中，可以看出：地井及井中瞬變電磁法的工作裝置的布置情況。
[0094] 井中TEM方法探測的地質目的在于探測分布于鉆孔附近的深部導電礦體，并獲得 礦體形態，產狀及位置等信息。發送回線通常布置地面，接收線圈（探頭）沿鉆孔井軸逐點 移動觀測磁場分量的微分參量。也可以把發送回線和接收回線同時放在井中，兩個線圈按 照一定的排列方式，沿井壁逐點移動進行觀測。當勘查區有彼此相靠近的多個鉆孔條件下， 一般只敷設一個大發送回線，從不同鉆孔中觀測到的異常變化規律可獲得地下隱伏導體的 位置等方面的信息。在僅有單個鉆孔的情況下，需要在地面敷設發送回線，根據Tx位于不 同方位上所觀測到的異常變化規律再去反演有關參數。也可以把發送線圈和接收線圈同時 都發置在井壁內，進行偶極形式觀測。
[0095] 參照圖7,示出了本發明實施例中航空瞬變電磁法的工作裝置示意圖。
[0096] 從圖7中，可以看出：航空電磁法的工作裝置的布置情況。
[0097] 發射線圈和接收線圈都在空中。航空TEM系統的發送線圈安裝于機身，接受線圈 及前置放大器安裝在吊艙之中，吊艙用電纜拖拽在飛機的后下部，即：發射和接收裝置都在 空中。飛行高度一般為150m。航空TEM方法主要應用于大面積范圍內快速普查良導電礦體 及地質填圖，在我國開展不多。
[0098] 步驟102、采用多點數據合成的方法對瞬變電磁的觀測數據進行處理及解釋，完成 對深部地質目標體的精細探測，獲得深部地質目標體的信息。
[0099] 在本發明的一種優選實施例中，所述步驟102,具體可以包括：
[0100] 子步驟112、對瞬變電磁的觀測數據，采用逆合成孔徑雷達算法的方式，進行相關 加權疊加，獲得瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體。
[0101] 在本發明的另一種優選實施例中，所述子步驟112,具體可以包括：
[0102] 子步驟A1、設計半空間中賦存高阻塊狀異常體的三維模型，獲得半空間中賦存高 阻塊狀異常體的三維模型；
[0103] 子步驟A2、設計半空間中賦存低阻塊狀異常體的三維模型，獲得半空間中賦存低 阻塊狀異常體的三維模型；
[0104] 子步驟A3、依據半空間中賦存低阻塊狀異常體的三維模型及半空間中賦存低阻塊 狀異常體的三維模型，采用逆合成孔徑雷達算法的方式，對瞬變電磁的觀測數據進行波場 變換，獲得瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體。
[0105] 子步驟122、對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行克?；舴蚱?移成像，獲得深層地質目標體的數字圖像；
[0106] 在本發明的另一種優選實施例中，所述子步驟122,具體可以包括：
[0107] 采用三維邊界元方式，對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行 克?；舴蚱瞥上瘢@得深層地質目標體的數字圖像。
[0108] 子步驟132、對所述深層地質目標體的數字圖像，進行數字圖像的分析，獲得深部 地質目標體的信息。
[0109] 為了使本領域的技術人員更好地理解本發明，下面詳細介紹采用多點數據合成 的方法對瞬變電磁的觀測數據進行處理及解釋，完成對深部地質目標體的精細探測，獲得 深部地質目標體的信息的過程。
[0110] (1)采用煙圈解釋方法，從觀測曲線得到視電阻率斷面剖面圖。
[0111] 參照圖8,示出了本發明實施例中瞬變電磁煙圈示意圖。
[0112] 從圖8中，可以看出：本發明實施例中瞬變電磁煙圈的工作原理。
[0113] 在電導率為〇和磁導率為的均勻大地上，敷設輸入階躍電流的回線，當發送 回線中電流突然斷開時，在下半空間中就要被激勵起感應渦流場以維持在斷開電流以前存 在的磁場，此瞬間的電流集中于Tx附近的地表，并按f4規律衰減（r為中心至觀測點的距 離）。隨后，面電流開始擴散到下半空間中，在切斷電流后的任一晚期時間里，感應渦流呈多 個層殼的環帶形，并形成一系列與發送回線同形狀并且向下及向外擴散的"電流環"，通常 稱之為"煙圈"。
[0114] 大地感應渦流在地表面產生的電磁場可近似地用圓形電流環表示。這些電流環就 像由發射回線吹出的"煙圈"，其半徑隨著時間增大而擴大，其深度隨時間延長而加深。
[0115] 當計算均勻半空間的瞬變電磁響應時，可以用一個鏡像電流環來代替。
[0116] 在層狀介質中，仍然保持同樣的"煙圈"效應，只是"煙圈"的傳播將逐漸局限于導 電地層中。
[0117] 隨著時間的延長，渦流場將向下及向外擴展。渦流場極大值將沿從中心起始與地 面成銳角的錐形斜面向下及向外傳播。
[0118] M.N.Nabighian指出，感應渦流場在地表引起的磁場為整個"環帶"各個渦流層的 總效應，這種效應可以用一個簡單的電流環相等效，當發送電流切斷以后某個時刻的地下 等效電流環時，它為一系列與發送回線同形狀并且向下及向外擴散的電流環，可以把它看 作一系列的二次發送線圈，可以計算某時刻的半徑，深度及電流，最終計算出某時刻的響應 值，以及其隨時間的變化規律。 i=--(1)
[0119] 它的等效電流為：4泥
[0120] 它的半徑表達式為(2)
[0121] 它的深度表達式(3)
[0122] 式中c2 = 8/JI-2 = 0· 546479
[0123] 由于tanΘ=d/a= 1.07,Θ= 47°，故"煙圈"將沿47°傾斜錐面擴展，向下傳 播的速度為 Cd2
[0124] ⑷ 〇tV辱丨，
[0125] 其中t為傳播時間，〇為介質電導率，μ^為真空中的磁導率。
[0126] 計算均勻半空間的瞬變電磁響應時，可以把"煙圈"看作一系列的二次發送線圈， 很容易地計算出在某時刻沿地面測線的響應值，以及在某個測點的響應值隨時間變化的規 律。在層狀介質中，仍然保持同樣的"煙圈"效應，只是"煙圈"的傳播將逐漸局限于導電地 層中。
[0127] (2)瞬變電磁場的波場變換，把電阻率斷在轉換成了擬地震剖面分布情況。
[0128] 參照圖9,示出了本發明實施例中觀測曲線額波場轉換示意圖。
[0129] 目前已經實現了波場轉換，并可以成功地提取出虛擬地震子波。
[0130] 其中，圖9中的（a)為現場觀測數據曲線示意圖；
[0131] 圖9中的（b)為現場觀測數據的加噪聲5%以后的曲線示意圖；
[0132] 圖9中的（C)為現場觀測數據曲線的波場轉換曲線示意圖；
[0133] 圖9中的（d)為現場觀測數據加噪聲以后的波場轉換曲線示意圖。
[0134] 從圖9中，可以看出：本發明實施例中觀測曲線額波場轉換情況。
[0135] 瞬變電磁場與虛擬波場之間存在如下對應關系
[0136] Hm (r，t)=--1__ [re- (r.r)dτ ('入/ 2V/Tt'
[0137] 式中，Hm(r，t)為時域瞬變響應擴散場，U(r，τ)為以波速^傳播的虛擬波場， 自變量τ是時間平方根的量綱。
[0138] 通過上述技術，本發明實現了從電阻率剖面到虛擬子波斷面的轉換。
[0139] 參照圖10,示出了本發明實施例中虛擬子波斷面對比示意圖。
[0140] 從圖10中，可以看出：本發明實施例中虛擬子波斷面對比的分布情況。
[0141] (3)壓縮子波寬度
[0142] 參照圖11，示出了本發明實施例中子波寬度壓縮前后對比示意圖；
[0143] 從圖11中，可以看出：本發明實施例中子波寬度壓縮前后的分布情況；經過壓縮 處理后，波形明顯得到改善。
[0144] 經波場變換獲得的虛擬子波存在嚴重的波形展寬現象，使得計算得到的虛擬波場 數值分辨率降低，嚴重影響著TEM成像的空間分辨能力。為此，對虛擬波場的離散數據求取 反褶積，消除波場變換的波形展寬效應。
[0145] 設U(r，τ)為實際轉換出來的虛擬子波，本發明通過一個反褶積濾波因子h(t)， 進行如下計算，求取一個寬度得到壓縮后的新的子波U' (r，τ)
[0146] υ'(τ,τ) = 人1，4 -τ)，
[0147] (6)
[0148] 其中，反褶積濾波因子h(t)可用最小平方反褶積來求得。
[0149] 隨著瞬變電磁波場變換算法的提出，實現了由具有擴散特征的瞬變場向虛擬波場 的轉變，這就為實現航空瞬變電磁法的逆合成孔徑成像創造了條件。瞬變電磁逆合成孔徑 成像技術是借用逆合成孔徑雷達成像的思路[3?8]，就是利用機載真實孔徑發射線圈與 目標的相對運動，把尺寸較小的真實天線孔徑用數據處理的方法合成一較大的等效孔徑的 發射線圈，使其分辨能力更高、穿透能力更強。對于航空瞬變電磁法而言其觀測方式與機載 逆合成孔徑雷達十分相似，完全可以借助逆合成孔徑雷達的成像思想，實現虛擬波場條件 下的航空瞬變電磁逆合成孔徑成像。
[0150] 雖然逆合成孔徑技術在雷達成像中早已廣泛應用，但與逆合成孔徑雷達信號不 同，瞬變電磁場滿足的是擴散方程，場的變化特征是隨著時間衰減的，顯然不利于相關疊加 處理，也不利于合成成像。為此，必須利用瞬變電磁場與虛擬波場間存在的數學上的關系式 進行波場變換，通過優化算法提取出這種虛擬波場信號。數字模擬和模型試驗的結果都證 明了：相鄰位置上同一地質體的反射回波具有較好的相關性，因此根據不同位置信號的相 關系數生成不同的權值函數，相鄰各列信號在做相關疊加時以權函數進行加權，將重建的 地質異常體信號加強，從而提高信噪比，達到突出弱異常的目的，進而提高了分辨率，加大 了勘探深度。在分析逆合成孔徑雷達算法的基礎上結合瞬變電磁信號的特點，對采樣信號 進行相關加權疊加形成瞬變電磁逆合成孔徑數據體，對該數據體進行克?；舴蚱瞥上瘢?得到地質體高清晰度的數字圖像。
[0151] 通過模型數據的計算，并進行逆合成孔徑成像處理，說明了該方法的有效性；對以 往山區地面數據的再處理，并與原來的處理結果相比較，說明了該方法在提高瞬變電磁分 辨力方面的優越性?？梢?，該方法的成功應用，將對用高分辨的航空瞬變電磁法代替海面和 山區地面勘探具有重要意義。
[0152] 由于瞬變電磁場經過波場變換，已經把原來的感應場轉換成了波場，每一點的數 據相當于變成了自激自收的波動場。而前人已經通過實驗分析確定瞬變電磁場在多激勵源 情況下存在場的相關疊加性，基于瞬變電磁場的上述特點，我們采用相關疊加的方法來進 行逆合成孔徑。其合成示意圖如圖1所示。
[0153] 參照圖12,示出了本發明實施例中逆合成孔徑示意圖；
[0154] 從圖12中，可以看出：本發明實施例中逆合成孔徑的工作原理。
[0155] 首先選取一個中心點，取為第i點，此點的波場值可表示為UOvτ)，其中ri為i 點到-N，"·，Ν內某點的距離，τ為相對時移量。然后我們選定2N+1個測點的長度為逆合 成孔徑的長度，即選取i點左右兩側從-N到N的測點分別與中心點做相關，其歸一化的互 相關系數

【權利要求】
1. 電性源瞬變電磁地空探測方法，其特征在于所述方法包括： 采用接地長導線源AB作為發射裝置，向地下發送一次場，利用直升機或者無人機搭載 接收線圈進行數據采集，獲得瞬變電磁的觀測數據； 利用全區視電阻率定義對觀測數據進行處理，獲得虛擬波場的速度信息，并完成對深 部地質目標的定性解釋，獲得深部地質體的概況信息； 對瞬變電磁數據進行波場變換，采用瞬變電磁虛擬波場連續速度分析及虛擬波場偏移 成像方法對瞬變電磁觀測數據進行處理解釋，完成對深部地質目標的精細探測，并且為了 提高數據的信噪比，采用多點數據合成的方法對瞬變電磁虛擬波場數據進行處理，實現深 部弱信息的提取，獲得深部地質目標體的詳細信息。
2. 如權利要求1所述的方法，其特征在于： 所述地質目標體的詳細信息為地質目標體的電性信息、深度信息、規模信息與形狀信 息。
3. 如權利要求1所述的方法，其特征在于， 地空瞬變電磁探測系統采用接地長導線源作為發射，利用直升機或者無人機搭載接收 線圈進行數據采集，其步驟包括： 通過在地面鋪設接地長導線源AB向地下發射一次電磁場； 在高度范圍為50-150米的空中設置觀測點，利用直升機或無人機搭載接收線圈進行 數據采集，獲得瞬變電磁的觀測數據； 其中， 所述對地質目標體進行探測，獲得觀測數據的過程，是通過電性源瞬變電磁探測裝置 完成的。
4. 如權利要求1或3所述的方法，其特征在于： 所述接地長導線源AB的長度為1-4千米； 所述接地長導線源AB的發射功率為不小于10千瓦。
5. 如權利要求1或3所述方法，其特征在于： 利用直升機或無人機搭載接收線圈進行數據采集，直升機或無人機飛行高度在保證飛 行安全情況下為30-150米， 飛行范圍與電性源尺度有關，垂向與軸向方向均不大于3倍電性源尺度。
6. 根據權利要求1所述的方法，其特征在于： 采用全域視電阻率方法對觀測數據進行處理，完成對深部地質目標的定性解釋，獲得 地質目標的產狀、走向、規模等概況信息，并獲得虛擬波場速度信息，包括： 利用反函數原理實現電性源瞬變電磁的視電阻率定義，通過時深轉換獲得深度信息， 最終完成初步的視電阻率成像，并獲得不同深度虛擬波場的傳播速度。
7. 根據權利要求6所述的方法，其特征在于： 所述利用反函數原理實現電性源瞬變電磁的視電阻率定義，通過時深轉換獲得深度信 息，完成初步的視電阻率成像，包括： 通過偶極疊加原理推導均勻半空間空中位置電性源瞬變電磁場響應表達解析式，通過 反函數原理實現視電阻率的定義； 通過濾波方法實現多分量層狀模型電磁響應的計算，利用電阻率定義公式計算層狀模 型視電阻形態，驗證定義方法； 利用視電阻率值與時間參數實現時深轉換。
8. 根據權利要求1所述的方法，其特征在于： 所述對瞬變電磁觀測數據進行波場變換，將瞬變電磁擴散場數據轉換為虛擬波場函 數，并利用相關疊加方法實現虛擬波場的逆合成孔徑處理，結合速度分析，對數據進行克希 霍夫偏移成像，獲得深層地質目標體的圖像，具體步驟，包括： 采用預調件正則化共軛梯度法實現波場反變換； 采用相關疊加方法，對虛擬波場進行逆合成孔徑處理； 采用三維邊界元方式，對瞬變電磁觀測數據的瞬變電磁逆合成孔徑數據體，進行克希 霍夫偏移成像，獲得深層地質目標體的圖像。
9. 根據權利要求1所述的方法，其特征在于，所述采用瞬變電磁虛擬波場連續速度分 析及成像方法，對瞬變電磁的觀測數據進行處理及解釋，完成對深部地質目標體的精細探 測，獲得深部地質目標體的信息的步驟，包括： 利用全域視電阻率定義方法，獲得瞬變電磁的觀測數據的虛擬速度數據； 利用速度建模的方式，獲得瞬變電磁的觀測數據的連續速度分析體成像圖像； 通過上述瞬變電磁虛擬波場連續速度分析及成像方法，對瞬變電磁的觀測數據進行處 理及解釋，完成對深部地質目標體的精細探測，獲得深部地質目標體的信息。
10. 根據權利要求8所述的方法，其特征在于，所述采波場變換方法，及利用預調件正 則化共軛梯度發實現擴散場到虛擬波場的變換，使得變換后函數滿足波動方程，包括： 利用預調件子對反變換系數矩陣進行預調件處理，降低矩陣條件數； 采用正則化共軛梯度迭代實現擴散場到波場的轉換。
11. 根據權利要求8所述的方法，其特征在于，所述采用速度分析的方式，獲得瞬變電 磁的觀測數據的連續速度分析體成像圖像的步驟，包括： 采用視電阻率定義的方式，對瞬變電磁的觀測數據計算相應電阻率，根據電阻率與虛 擬波動速度之間的關系獲得瞬變電磁虛擬波場的初始速度； 采用近點線性插值方法，擴大速度數據體，實現連續速度分析，使速度數據滿足后續波 場延拓需要； 其中，所述波場延拓指利用克?；舴蚍e分實現波動方程的向下延拓成像深。
【文檔編號】G01V3/165GK104237956SQ201410081433
【公開日】2014年12月24日 申請日期:2014年3月6日 優先權日:2014年3月6日
【發明者】李貅, 戚志鵬, 張瑩瑩 申請人:長安大學