電磁感應和誘導極化效應的量化分解方法
2023-07-10 07:23:26
專利名稱：電磁感應和誘導極化效應的量化分解方法
技術領域：
本發明屬於電法勘探研究領域，是一種地下巖層性質的測定方法，同時還對地層的電學性質參數進行測量，並分解測得參數。本發明用於陸上及海洋中的可控電磁場源電法勘探，適用於尋找並圈定油氣藏的綜合性油氣勘探工作。它對次生蝕變巖在油氣遷移作用下的回應進行分離，並以其為基礎將測得信號分解成某些分量。本發明可保證一系列用於分層確定介質的地電參數值、誘導極化和電磁感應過程特徵的技術方法。
背景技術：
已知發明《地電勘探方法(方案)》，專利號2235347RU，發表日期2004年8月27 日，國際專利分類號G01V3/06。根據該發明，在被研究介質的深層周期性依次通過矩形電流脈衝，並在每兩次脈衝之間保持間歇，從而激發電磁場。在整個過程中，於每次電流脈衝的結尾和兩次脈衝之間的間隙測量軸向和正交電勢差的瞬時值，確定電勢差。按公式計算η 個歸一化參數，根據極化導電介質中的偶極子源場強的數學物理差分波動方程運算反解問題，建立其幾何結構和電學參數最接近被研究介質的介質模型根據該方程中包含的電物理參數，建立該模型的時間剖面圖。利用該發明可發現並圈定油氣藏、評估其飽和質量，並可分解電導率參數和誘導極化參數。然而，該發明不能解決分別量化研究電磁感應過程和誘導極化過程的重要問題，未考慮所建立模型的電物理參數的等值性，不能確定這些參數的 fn息各裡ο已知發明《使用聚焦電流的海洋地電勘探方法(方案)》，專利號RU 2284555，發表日期2006年9月27日，國際專利分類號G01V3/06。根據該發明，在被研究介質的深層激發電磁場，測量第一和第二電勢差的瞬時值，使沿剖面的軸向合成電勢差等於零並從等式中求得因數，計算三個獨立的歸一化電學參數集合，進行反解數學運算，求得三個電物理參數電導率、誘導極化度和誘導極化電勢差的衰減時間常數，並按這些參數建立三個時間剖面圖。該發明的本質在於分解電導率參數和誘導極化參數，作為補充，還可確定與前兩個參數同樣重要的第三個參數——誘導極化電勢差的衰減時間常數。然而，該發明僅能獲得有關被研究介質所有結構元素的總和信息，這是因為電磁場源的被測電流I在介質中的空間分布不受任何監測，沒有在真實存在的三維非均質介質中關於該分布的信息。該發明不能解決分別量化研究電磁感應過程和誘導極化過程的重要問題，未考慮所建立模型電物理參數的等值性，不能確定這些參數的信息容量。已知發明《地電勘探方法(方案)》，專利號2231089，發表日期2003年7月8日， 國際專利分類號G01V3/06。根據該發明，在被研究介質的深層周期性依次通過矩形電流脈衝，並在每次脈衝之後保持間歇，從而激發電磁場。在觀測點測量電勢差，使用歸一化參數值和電化學極化導電介質中的偶極子源電場場強的數學物理差分波動方程，運算反解問題，得出每個介質單元所固有的電物理參數電導率、誘導極化和誘導極化電勢差的衰減時間常數，並按這些參數建立時間剖面。利用該發明可發現並圈定油氣藏、評估其飽和質量，並可分解電導率參數和誘導極化參數。然而，該發明不能將地質介質單元所固有的誘導極化與瞬變電動力學過程徹底分解開來，該電動力學過程同所示地質巖石的剖分單元的電導率相關。該發明亦不能解決分別量化研究電磁感應過程和誘導極化過程的重要問題，未考慮所建立模型電物理參數的等值性，不能確定這些參數的信息容量。已知專利《地電勘探信號的分析方法》，專利號2253137RU，發表日期2005年1月 10日，國際專利分類號G01V3/38。根據該發明，將每個測量信號同時用至少兩個調製器轉換，然後找出最小的歩距間隔。該發明通過對信號(固定分量)和幹擾(隨機分量)的同時分析，能夠提高準確性和信息容量。然而，該發明的目的是利用該方法在計算機中對比增量調製器集成電路的輸出信號和相應通道的輸入信號。而這種分析方法的使用，卻必須具備至少兩個量子化歩距差別很大的增量調節器、一個接收電極系統、一套磁場傳感器或者具有足夠寬頻帶性和空間角度選擇性的任何其他傳感器。該方法不能用於分析在使用誘導極化過程尋找和勘探油氣藏時獲得的參數。已知發明《充滿井內流體的井筒周圍地質構造的電阻率測定方法》，專利號 2209451RU，發表日期2003年7月27日，國際專利分類號G01V3/38。該發明包括獲取某些圖表、選擇模擬剖面、引入附加模擬的剖面、重複修正模擬剖面直至圖表和相應模擬圖之間的差值小於選定的閾值、將每個模擬剖面作為值的關係函數進行修正等。該發明考慮了井內流體對井筒周圍構造的入侵，完善了地質構造電阻率的測定方法。然而，該發明不能發現並圈定油氣藏、評估其飽和質量，並分解電導率參數和誘導極化參數。此外，該發明亦不能解決分別量化研究電磁感應過程和誘導極化過程的重要問題，未考慮所建立模型電物理參數的等值性，不能確定這些參數的信息容量。同本發明在地電性徵的處理方法上最接近的是《地下構造的性質測定方法》發明， 專利號2294M7RU，發表日期2007年2月27日，國際專利分類號G01V3/38。該發明包括構造參數模型的擬定、模擬電學性質的測量、將模擬測量同已測得電學性質相比較、以各種測量為基礎確定構造的物理/石油物理性質以保證聯合反演、利用電學性質模型和各種測量評估構造的物理/石油物理性質，還包括使用綜合模型。該發明可保證在應用電力對鑽井進行剖分時能夠實施更完善的數據處理方法，並且更準確地確定Rt值以評價構造特徵，在此過程中使用將初始信息引入反演圖的構造層系測定。初始信息即穩定未精確設置的「可能性擴展」任務的幾何信息。然而，該發明不能發現並圈定油氣藏、評估其飽和質量，並分解電導率參數和誘導極化參數。此外，該發明亦不能解決分別量化研究電磁感應過程和誘導極化過程的重要問題，未考慮所建立模型的電物理參數的等值性，不能確定這些參數的信息容量。該發明有另一個目的，即當介質模型可通過迭代反演方法進行簡單明確時實現該方法，不將測得信號分解成某些分量。而在我們的發明中，通過反解運算得出的介質模型用於將測得信號AU(t)分解成分量EM (t)和IPg (t)。
發明內容
當介質中通過電流時，會產生一些由不同物理現象所決定的效應。此時會發出電磁感應，其特點是符合波動理論的感應場，並且發生由介質中的電化學過程所決定的誘導極化。到目前為止，對這兩種效應的分解辦法是增加記錄介質中激發脈衝瞬變的時間，而電磁感應和極化正是由該脈衝所引起的。然而，隨著記錄時間的增加，電磁感應場在測量之前幾經衰減到了幾乎可以忽略不計的微小數值。相反，如果縮短測量時間，在激發脈衝之後的初始時間內，極化效應幾乎沒有。在測量沉積蓋層較厚的地臺時，由於沉積蓋層的導納較大(數百或數千西門子)，電磁波衰減會比較慢。該問題的存在及其解決方法的缺失在J.維特的書中有提到(《大地電磁學》，James R. Wait，莫斯科，地址出版社，1987年，235 頁)。本方法的技術成果是它可以解決分別量化研究電磁感應場和誘導極化的重要問題，考慮了所建立模型的電物理參數等值性，能夠確定其信息容量。本方法能夠一量化測定直流電的電磁感應場強和誘導極化場強，誤差不超過0. 5% ；一保證電磁場分解從1微秒至數秒的較大時間間隔，包括針對高導納的地質剖分。由於函數AU總是解析性的，即總是連續的，並且其二階導數亦為連續，因而本方法是可用的。此外，由於部分誘導極化場由感應電流引起，而部分感應場反過來又由極化場引起，所以電磁感應場和誘導極化場是非相加性的。我們將只討論直流電在通過介質的時間段內所引起的誘導極化。對所聲明技術成果的保證，是因為電磁感應和誘導極化效應的量化分解方法包括一測量極化介質的瞬變過程；一水平層狀極化介質範圍內的電磁測量模擬；一從瞬變過程中分離出電動力學分量；一從瞬變過程中分離出誘導極化的電流分量；一評估誘導極化電流分量的局域變化，對其進行地質解釋。本方法的特點是記錄幾個取決於介質參數的、在介質表面測量的輸入函數，這些函數對電磁感應場和誘導極化場有不同的依賴關係，同時其中一個函數(例如Ps)同DU(t) 相比，可增大電磁感應和誘導極化的比值，而另一個函數(例如P1)同DU(t)相比，則可縮小電磁感應和誘導極化的比值，還有一個函數(例如 Φ)則是瞬變場的時間導數和空間導數的組合。然後，同時反演所有輸入函數，得出剖面的地電模型，並且此時表示誘導極化的介質參數數量明顯過剩，因此得出的介質模型不是唯一的，而是位於某個等值區域內。在得出的模型中使每一層的極化率都等於零，通過正解問題計算出電磁感應場EM ；使同一模型中的波數等於零，正解問題計算出誘導極化的電流分量場IP。此外，舉例而言，建立參數 P1，將其作為同Δυα)相比提高誘導極化/電動力比值的篩選參數，以瞬變場第二軸向電勢差Δ2υα)同瞬變場第一軸向電勢差Δυα)在斷電制式下的比值計算P1，計算公式為 P1 (t) = A2U(t)/AU(t) Itiffle。ff ；或者，舉例而言，建立參數1^，將其作為同AU(t)相比降低誘導極化/電動力比值的篩選參數，以瞬變場第二軸向電勢差八2U(t)同瞬變場第一軸向電勢差AU(t)在通電制式下的比值計算Ps，計算公式為Ps(t) = A2U(t)/AU(t) Itime 。n;再或者建立參數Dcp(t)/作為按公式Dct(t) = It (O-P1U)計算的參數It (t)同參數 P1 (t)的差值，並且It (t)作為時間導數Δ2υα)同Δυα)的比值，其計算公式為It(t)= At(A2U(t))/At(AU(t))。
在建立同一記錄點測得的輸入函數之後，將它們最小化，為此在輸入時給出參數的最小化程序—DU(t) = AU(t)/AUQ，式中Δ Uq-電流通過時測量的電勢差AU;-P1 (t) = Δ勺(t) / Δ U (t)，在斷電制式下進行計算；—Ps(t) = A2U(t)/AU(t)，在通電制式下進行計算;—DcHt) = It (O-P1U)，式中It(t) = At(A2U(t))/At(AU(t))，並對它們在水平分層極化介質範圍內進行反演，並且此時表示模型極化性質的地電參數數量明顯過剩。同時作為一種特殊情況，正解波數等於零求出IP場，正解極化率等於零求出EM場。完成所有運算之後，檢查解的可重複性。為解決所設置的效應分解任務，需要完成以下幾個步驟。第I階段記錄幾個可在介質表面測量並建立的輸入函數。函數的設計應使它們對電磁感應和誘導極化有不同的依賴關係。其中一個同電壓降曲線DU相比，可增加電磁感應同誘導極化的比值，而另一個則縮小該比值。這些函數是瞬變場同其空間導數的組合。此外，還有一個函數是時間導數和空間導數的組合。第II階段1.同時反演(反解問題)所有輸入函數，建立一個剖面模型。同時逐層求解介質參數電阻率P，極化率η，弛豫時間τ，冪c。極化率、弛豫時間和冪根據Cole-Cole模型確定。利用這些參數對模型進行參數化時應使之明顯過剩，意即獲得等值性的某些解。2.在建立的模型中使每層的極化率等於零，正解問題求得電磁感應場EM。3.在同一模型中使波數等於零，正解問題求得直流電的誘導極化場IP。重複運算(第II階段)顯示計算出的EM場和IP場重合，並精確至反解運算的解。 所得出的模型是不相同的，但場是重合的。其中，IP——使剖面任意點的波數等於零，對該模型進行正解運算得出的誘導極化場；EM——使剖面任意點的極化率等於零，對該模型進行正解運輸得出的電磁感應場。
本方法用圖表進行說明和證明。圖1是AU1和AU2的曲線，即偶極子軸向裝置在接收線OM和ON上的電勢差。這個圖形顯示存在著IP場和EM場不同組合的無窮集。對合成模型計算了衰減曲線DU並進行了兩次不同零點逼近的反演，精確至0.1%。得到兩個不同的等值(按DU參數)解，其中每一個符合顯著不同的IP場和EM場的值。圖2顯示的是在瞬時過程不同時期(不同相)的EM和IP電流分布間隔初期(A)、 過渡期⑶和晚期(C)。此時EM場服從集膚效應，於瞬態過程中在下半空間內擴散，而由其引起的渦流則趨向於在下半空間內呈均勻分布。IP場則相反，總是呈空間不均勻分布，並且其強度隨著同電流源(發射源)的距離增大而減弱。根據所設置問題的涵義，需要為瞬變場引入可在同一點記錄的轉換參數，並且這些轉換參數對誘導極化場和感應場具有不同的函數依賴。最自然的解決方法——引入瞬變場的歸一化空間導數(更準確來講是有限差分)，做為對瞬變場本身的補充。這是因為感應場和誘導極化場的空時結構是不同的。在圖3中顯示的是a).測量裝置在輪船移動時記錄DU(t) ,P1 (t)、Ps(t)、D(p⑴參數的示意圖；b).瞬變激發電流在振蕩線AB上的形狀；c).在接收線上測量的信號A2U(t) η AU(t)。圖4顯示的是DU (t)曲線在地球表面上隨時間的衰減，包括存在(B)和不存在(A) 極化、在接收線M0(距離為IOOOm(I))和接收線0N(距離為2400米(2))上的情況。圖5顯示的是參數P1 (t)在地球表面上的曲線，包括存在( 和不存在(1)極化的情況。條件是 η=1%，Ρ=50Ω·π^Ριι=0%。圖6顯示的是參數Ps在地球表面上的曲線，相對於條件η = 4%和η = 0，其中 η-在Cole-Cole模型術語中的極化係數。圖7顯示的是在η =0和η = 1 %條件下計算的Dip曲線，其中η -在Cole-Cole 模型術語中的極化係數。曲線相對於Dcp = 0的直線建立。除了空間性質外還使用時間導數，從而獲得更詳細的瞬時過程特徵。圖8顯示的是三維參考模型。利用它計算了符合實地測量所得到的參數綜合曲線 (DU (t)、P1 (t)、Ps (t)、Dcp(t))，以及EM場和IP場的值。並且給出了與其等效的一維模型。圖9顯示的是EM場和IP場，它們同圖8中所給出的模型相吻合，精確度為0. 15%。圖10顯示的是在北古利亞耶夫斯卡規區(巴倫支海大陸架)進行實地測量時的觀測曲線及與其相符的模型。給出了等效地電模型之一。圖11顯示的是在北古利亞耶夫斯卡規區(巴倫支海大陸架)進行實地測量時的觀測曲線及與其相符的模型。給出了與圖10不同的等效地電模型之一。圖12顯示的是從等效模型中得出的EM場和IP場的衰減曲線。它們彼此間的差值要小於測量誤差，即小於0. 5%。此時IP場從在大範圍時間間隔內比其高2個數量級的 EM場中分離出來。圖13顯示的是已知油氣田的IP場剖面等值線和油氣藏沿IP場的劃定。
具體實施例方式已知對地球表面進行測量獲得的一系列輸入數據。數據的組成應使它們對誘導極化和電磁感應有不同的函數依賴關係。以這些數據為基礎，在表示誘導極化過程的剖面參數數量過剩的條件下反演得出地電模型。如果在這個模型中，令波數為零正解得出IP場，令極化率為零正解得出EM場，則此時的IP場和EM場符合可重複性條件。這就意味著對不同等效模型進行IP場和EM場計算所得出的偏差，不會超過實際觀測誤差(約0. 5% )。為實現本方法，1.在電流脈衝切斷制式下測量時建立參數Pl (見圖1)，將其作為同AU(t)相比增大誘導極化/電動力比值的篩選參數。一為此以瞬變場第二軸向電勢差Δ2υα)同瞬變場第一軸向電勢差Δυα)在斷電制式下的比值計算Ρ1，BP =P1U) = A2U(t)/AU(t)f。
在圖1中顯示了八仏和—偶極子軸向裝置在接收線OM和ON上的電勢差。 參數P1對瞬變場的空間不均勻性作出反應。電磁場以擴散方程表示，因此在增大衰減時間時會趨向於在介質中均勻分布。誘導極化場正比於極化電流密度，因此它總是呈空間不均勻性(圖2)[列格依多P. Yu.，《油氣地球物理中研究極化剖面的差分歸一化地電勘探理論和技術》。地質與礦物學正博士學位答辯論文。1998年]。2.在電流脈衝接通制式下測量時建立參數1^，將其作為同AU(t)相比降低誘導極化/電動力比值的篩選參數。一為此計算瞬變場第二軸向電勢差A2U同第一軸向電勢差AU在通電制式下的比值。在瞬變初期的誘導極化/電動力比值很小，而在晚期階段I^s參數基本上由阻抗的分布來決定，極化率對它的影響很小。因此I3S參數可被視為同瞬變場AU相比較能夠縮小誘導極化/電動力比值的篩選參數(圖6)。3.建立參數Dcp，以便更詳細地表示瞬變過程，因為這個參數以時間導數A2U(t) 和Δυα)為基礎計算。一它是相對於作為瞬變時間對數的假定單位的參數差，即參數It同參數P1的差值=Dcp = It-P"—其中，It=At(A2U)/At(U)，即時間導數Δ2_同Δ_的比值。誘導極化場和電磁感應場的衰退時間一般是不同的，因此有意識地使用時間參數，以更詳細地描述瞬變過程(圖7)。4.對在同一點記錄的參數進行最小化。一為此在輸入時給出參數的最小化程序(圖3)；—DU(t) = Δυα)/Δυ。，式中Δ Uq-電流通過時測量的電勢差AU;-P1 (t) = Δ 2U/ Δ U (t)，在斷電制式下進行計算；-Ps (t) = Δ 2U/ Δ U (t)，在通電制式下進行計算；— Φ (t) = It-Pl，式中:It = Δ t ( Δ 2U) / Δ t ( Δ U) ο在模型的地電參數數量明顯過剩的條件下，對它們在水平分層極化介質範圍內進行反演，以儘可能準確地接近對非均質各向異性介質和梯度介質記錄的瞬變場曲線。5.按建立的介質模型計算IP場和EM場。一使模型在介質所有區域的極化率等於零，正解問題求出電磁感應場ΕΜ，—使模型在介質所有區域的波數等於零，正解問題求出誘導極化場IP。6.檢查解的可重複性。一對EM場和IP場的重複性檢查顯示，它們實際上對所有種類的剖面都是可重複的，並且誤差不超過測量誤差(約0. 5% )(見圖10、13)。感應場和極化場的分解方法，可以用以下示例進行展示。在三維模型中獲得符合實地測量的參數綜合曲線和EM場與IP場的值(圖8)。對於綜合曲線，已在水平分層極化介質範圍內建立了等值模型，以其為基礎計算了 EM場和IP場(圖8)。對於從三維模型和水平分層模型中所得的場值進行了比較。它們之間的誤差 (0. 15% )不超過三維模型綜合曲線與水平分層極化介質模型曲線之間的偏離誤差。由於偏離很小，曲線之間的差別幾乎看不出來。實踐當中，對電磁感應效應和誘導極化效應的量化分解，在各種不同的地質和地球物理條件下進行了實驗。對已知的北古利亞耶夫斯卡規區(巴倫支海大陸架)的油氣田進行了電法勘探， 測量了參數DU (t)、P1 (t)、Ps (t)、D(p⑴。為場曲線建立了水平分層極化模型，逼近誤差小於 0.5% (圖11、12)。此時，為同一個測量標點建立了幾個等效模型，並且計算出的EM場和 IP場的圖形偏離小於0.5% (圖13)。在此條件下，成功分離出了比EM場小兩個數量級的 IP場。通過建立IP場的等值線，清除地顯示出了油氣藏(圖14)。其中Ro-阻抗，Ω · m ；C-在Cole-Cole模型中的冪；H-層的厚度，m;η-Cole-Cole模型的極化率係數，％t-Cole-Cole模型的弛豫時間，S。例如以下參數的模型Ro = 10 Ω · m，η = 5%, t = 0. 5s, c = 0. 5, H = IOOm ；Ro = 50 Ω · m，η = 2%, t = 0. 5s, c = 0. 5, H = IOOm ；Ro = 20 Ω · m，η = 2%, t = Is, c = 0. 5, H = 200m ；Ro = 5Ω · m，η = 5%, t = Is, c = 0. 5, H = 300m ；Ro = 200 Ω · m，η = 7%, t = Is, c = 0. 5 ；Ro = 500 Ω · m，η = 0，H = 450m ；Ro = 5 Ω · m，η = 0，H = 800m ；Ro = 1000Ω · m，η = 0，H =⑴。對該模型計算了符合實地測量所得的輸入數據值，並計算了 IP場和EM場，為巴倫支海大陸架北古利亞耶夫斯卡規油氣田建立了曲線圖(見圖13)。換句話說，根據所設置問題的涵義，為瞬變場引入了轉換參數，這些轉換參數理論上可在同一點記錄得到，並且對誘導極化場和感應場具有不同的函數依賴。同時引入瞬變場的歸一化空間導數(更準確來講是有限差分)做為對瞬變場本身的補充，因為感應場和誘導極化場的空時結構是不同的。 從而得出了所聲明的技術成果。
權利要求
電磁感應和誘導極化效應的量化分解方法
1.一種電磁感應和誘導極化效應的量化分解方法，包括對極化介質進行瞬變場測量， 在水平分層極化介質範圍內模擬電磁測量，從瞬變場中分離出電磁分量，從瞬變場中分離出誘導極化的電流分量，評估誘導極化電流分量的局域變化，對其進行地質解釋；其特點是記錄若干個依賴於介質參數的、在介質表面測得的輸入函數，並且這些函數同電磁感應場和誘導極化場有著不同的依賴關係；其中一個函數同DU(t)相比，可提高電磁感應同誘導極化的比值；另外一個函數同DU(t)相比，可降低電磁感應同誘導極化的比值；第三個函數則是瞬變場時間導數和空間導數的組合；然後對同一記錄點錄得的所有輸入函數進行同時反演，包括DU(t)，從而建立介質剖面的地電模型；此時表示誘導極化的介質參數數量明顯過剩，因此建立的介質模型不是唯一的，而是位於某個等值區域內；在建立的模型中使每一層的極化率等於零，計算正解問題得出電磁感應場EM ；同時在同一模型中使波數等於零，正解問題求得誘導極化的電流分量場IP。
2.根據權利要求1所述的量化分解的方法，其特點是建立一組輸入函數，並且其中一個函數為Ps，通過與DU(t)相比增大電磁感應同誘導極化的比值進行計算；另一個函數是 P1,通過與DU (t)相比縮小電磁感應同誘導極化的比值進行計算；第三個函數是DcK通過組合瞬變場時間導數和空間導數進行計算。
3.根據權利要求1所述的量化分解的方法，其特點是建立參數P1，將其作為同AU(t) 相比增大誘導極化/電動力比值的篩選參數，以瞬變場第二軸向電勢差八2U(t)同第一軸向電勢差Δυα)在斷電制式下的比值進行計算，計算公式為=P1U) = A2U(t)/AU(t) ItiffleOff0
4.根據權利要求1所述的量化分解的方法，其特點是建立參數1^，將其作為同Δυα) 相比縮小誘導極化/電動力比值的篩選參數，以瞬變場第二軸向電勢差同第一軸向電勢差Δυα)在通電制式下的比值進行計算，計算公式為Ps(t) = A2U(t)/AU(t) Itimeon0
5.根據權利要求1所述的量化分解的方法，其特點是建立參數Dcp(t)，作為按公式 Φω = It (O-P1U)計算的參數It (t)同參數P1U)的差值，並且It(t)作為時間導數 A2U(t)同 AU(t)的比值，其計算公式為It(t) = At(A2U(t))/At(AU(t))0
6.根據權利要求1所述的量化分解的方法，其特點是對在同一記錄點測得的輸入函數同時進行最小化，為此在輸入時給出參數的最小化程序—DU(t) = Δυα)/Δυ。，式中AUq-電流通過時測量的電勢差AU;—Ρι⑴=Δ AU(t)，在斷電制式下進行計算;-Ps (t) = Δ 2U(t)/AU(t)，在通電制式下進行計算；—D(Ht) =It(t)-Pl(t)，S*:It(t) = At(A2U(t))/At(AU(t))，並對它們在水平分層極化介質範圍內進行反演，而且此時表示模型極化性質的地電參數數量明顯過剩。
7.根據權利要求1所述的量化分解的方法，其特點是使波數等於零，計算對模型進行最小化得到的正解問題，從而算出IP場。
8.根據權利要求1所述的量化分解的方法，其特點是使極化率等於零，計算對模型進行最小化得到的正解問題，從而算出EM場。
9.根據權利要求1 8所述的量化分解的方法，其特點是檢查解的可重複性。
全文摘要
本發明涉及在用受控制的電磁場源進行陸上和海上電法勘探期間判定地下地層的性質的方法，該方法用於石油和天然氣勘探作業，通過隔離來自在油氣運移的作用下已遭受次生變化的巖石的響應以便將測量信號劃分成任何成分來勘探並測繪石油和天然氣礦床。本發明使得提供這樣一批處理方法成為可能，這批處理方法用來逐層判定地電參數值和介質與電磁感應的誘導極化的處理特性。該方法包含根據表面面積及其地質判讀來評估誘導極化的流電成分的變化，同時記錄幾個依賴於在介質表面測量的介質參數的輸入函數，並根據所述參數形成所述函數使得所述函數以不同方式取決於電磁感應場和誘導極化場。
文檔編號G01V3/06GK102265188SQ200880132434
公開日2011年11月30日 申請日期2008年11月17日 優先權日2008年10月23日
發明者伊萬·約葉維奇·別斯捷列夫, 葉夫根尼·烏拉第米羅維奇·阿根科沃, 彼得·約葉維奇·樂哲都, 斯維特拉娜·約葉維娜·格瑞娜, 由瑞·亞歷山大羅維奇·達維登科, 艾琳娜·歐蓮格瓦娜·庫德裡耶夫塞娃, 謝爾蓋·亞歷山大羅維奇·伊萬諾夫 申請人:西伯利亞地球物理學研究和生產有限責任公司