該發(fā)明屬于石油勘探地球物理壓裂地震數(shù)據(jù)處理方法領(lǐng)域，尤其涉及水力壓裂震源定位。
背景技術(shù)：
：在我國，新油田的發(fā)現(xiàn)越來越困難，老油田開發(fā)利用率不高，低滲透油田或超低滲透油田儲量豐富，另外頁巖氣的開采力度也越來越大。水力壓裂是一種比較實用的增產(chǎn)方法，這個過程中就涉及到壓裂監(jiān)測的問題。微震監(jiān)測水力壓裂項目實施時，隨著壓力增加，當(dāng)超過巖石承受極限時破裂，產(chǎn)生地震波。破裂形成的裂縫在地下，很難通過直觀手段進(jìn)行觀測得到裂縫空間位置以及形態(tài)，具體表現(xiàn)在裂縫的高度、長度、寬度、延伸、方位，而這些裂縫參數(shù)正是壓裂技術(shù)人員和采油公司想知道的，通過監(jiān)測以上參數(shù)可以保證油氣開產(chǎn)過程中的安全性，也可以即時改變壓裂參數(shù)和方案，最終實現(xiàn)穩(wěn)產(chǎn)和高產(chǎn)的目的。微震監(jiān)測是一種有效的裂縫監(jiān)測技術(shù)，與其他裂縫監(jiān)測方法相比，監(jiān)測結(jié)果準(zhǔn)確，信息豐富。國內(nèi)外學(xué)者對水力壓力微震定位方法已做了大量研究，比較常用的技術(shù)是根據(jù)射孔和測井資料建立符合靶區(qū)的地質(zhì)模型，然后求解出微震位置。微地震反演定位方面，在實際生產(chǎn)實踐中，目前技術(shù)主要利用直達(dá)縱波的初至信息進(jìn)行震源定位，在某些特定的觀測系統(tǒng)下，會給定位帶來很大誤差。反演方法選取上，有一類方法通過震源和檢波器位置關(guān)系建立方程組或?qū)δ繕?biāo)函數(shù)進(jìn)行線性化，然后通過求解方程組來獲取微震的真實震源位置，由于初至拾取誤差和檢波器排列等影響，結(jié)果往往不能收斂或定位精度差；另外一種比較有代表性并且運用于商業(yè)化的微震定位軟件主要選取直接或者快速網(wǎng)格搜索算法，該方法當(dāng)監(jiān)測區(qū)域范圍比較大時，計算效率低，而且定位精度與網(wǎng)格大小有關(guān)。因此，有必要對微地震震源定位方法進(jìn)行深入研究，提高定位精度。技術(shù)實現(xiàn)要素：本發(fā)明的目的在于為了克服以上所述存在的問題，提供一種適用于水力壓裂微震震源定位的方法，以提高水力壓裂過程中微震源的定位精度。本發(fā)明的適用于水力壓裂微震震源定位的方法，包括以下步驟：(1)建立微震監(jiān)測的觀測系統(tǒng)，其中坐標(biāo)和線道號參數(shù)通過現(xiàn)場獲取，同時利用常規(guī)聲波測井或全波列測井得到聲波資料和密度曲線，結(jié)合測井資料在層狀介質(zhì)假設(shè)下構(gòu)建初始速度模型；(2)在建立微震監(jiān)測的觀測系統(tǒng)后，分析研究區(qū)域地層條件和射孔記錄，根據(jù)背景干擾頻率低、微震信號頻率高、持續(xù)時間短的特點確定背景干擾和微震信號，對壓裂過程中的記錄進(jìn)行常規(guī)的預(yù)處理，然后采用頻率域濾波進(jìn)行壓噪處理，最后用長短時窗能量比法對去噪后有效微震信號進(jìn)行初步波形識別；(3)根據(jù)確定的射孔位置和時間，采用人工干預(yù)方法識別射孔記錄，利用初至自動拾取和人機(jī)交互方法對直達(dá)縱橫波的初至走時進(jìn)行拾取，在步驟(1)構(gòu)建的初始速度模型約束下改變每層的厚度和速度，對速度模型進(jìn)行校正；而其他在壓裂過程中被記錄識別到的有效事件用于震源的反演定位，利用同樣的方式獲取縱橫波初至走時；(4)在利用縱波走時的基礎(chǔ)上，由測井獲取得到了橫波速度，此時加入橫波走時信息，結(jié)合校正后的速度模型，利用縱波走時和橫波走時同時約束進(jìn)行微震震源定位，反演采用全局搜索和局部搜索相結(jié)合的方法，包括：①根據(jù)水力壓裂段的位置坐標(biāo)，確定震源定位搜索范圍；②初始化震源參數(shù)，輸入檢波器坐標(biāo)和步驟(3)中拾取到的縱橫波走時信息，在搜索范圍內(nèi)，隨機(jī)產(chǎn)生一個點作為反演初始震源點，并正演計算縱波走時和橫波走時，公式如下：T=Σi=1m+1li/vi=Σi=1m+1[(xi-xi-1)2+(yi-yi-1)2+(zi-zi-1)2]1/2/vi---(1)]]>其中m表示檢波器(xm+1,ym+1,zm+1)和假設(shè)震源點(x0,y0,z0)與建立地層模型交點個數(shù)，(xi,yi,zi)(i～m)為上述交點的空間坐標(biāo)，vi表示射線段對應(yīng)的地層橫波或者縱波速度；③建立并計算目標(biāo)函數(shù)：F=Σi=1n((Tscal(i)-Tpcal(i))-(Tsobs(i)-Tpobs(i)))2---(2)]]>其中，Tscal(i)、Tpcal(i)分別為根據(jù)已知模型和②中假設(shè)點正演計算得到的橫波走時和縱波走時，Tsobs(i)、Tpobs(i)分別為通過步驟(3)獲取到的各檢波器對應(yīng)的橫波和縱波時間值，n為檢波器個數(shù)，在目標(biāo)函數(shù)中加入橫波信息可以提高定位結(jié)果的精度；④然后判斷根據(jù)公式(2)計算的目標(biāo)函數(shù)值是否滿足預(yù)先給定的精度，如果滿足，則轉(zhuǎn)到⑦，否則，進(jìn)入下一步；⑤設(shè)定模擬退火搜索迭代終止條件，采用能找到全局極小值的模擬退火法進(jìn)行全局搜索，利用公式(1)計算走時，利用式(2)計算目標(biāo)函數(shù)，每次迭代過程中都判斷是否滿足模擬退火終止條件，如果滿足，則進(jìn)入下一步；⑥設(shè)定模式搜索算法的迭代終止條件，以模擬退火的輸出結(jié)果為初始輸入模型，然后用模式搜索算法進(jìn)行局部搜索，采用基于正交基的模式搜索法，模式為：Pk＝BCk(3)其中，k表示迭代次數(shù)，B為模式的基矩陣，為一非奇異矩陣，Ck為一個生成矩陣，Pk為計算得到的模式矩陣,通過當(dāng)前點和模式計算新的搜索點的目標(biāo)函數(shù)值，如果新搜索值的目標(biāo)函數(shù)小于前面所有點的目標(biāo)函數(shù)，則通過(3)式重新修改模式，繼續(xù)新的搜索，直到進(jìn)行完所有的模式搜索，最后判斷定位結(jié)果精度是否滿足，如果滿足，則進(jìn)入⑦，否則，返回到步驟⑤中，重新進(jìn)行全局搜索；⑦定位結(jié)果是否滿足實際地質(zhì)條件，如果滿足，則輸出微震震源位置。上述方案還包括：步驟(4)確定震源定位搜索范圍為：深度范圍選取上下100米，而平面范圍為400米。上述方案進(jìn)一步包括：步驟(1)所述的測井資料利用已有的常規(guī)聲波測井得到聲波資料。步驟(2)所述的射孔記錄為采用聚能器進(jìn)入井眼預(yù)定層位進(jìn)行爆炸開孔讓井下地層內(nèi)流體進(jìn)入孔眼時儀器記錄到的振幅強(qiáng)的微震信號；步驟(2)中采用的壓噪處理為帶通濾波算法。步驟(3)中初至自動拾取為長短時窗比和基于信息準(zhǔn)則相結(jié)合的方法。步驟(4)的橫波速度從全波列測井、偶極橫波測井或者研究區(qū)域經(jīng)驗公式獲得。其中，長短時窗比的方法為公式(4)，STAi=STAi-1+CF(i)-CF(i-Nsta)NstaLTAi=LTAi-1+CF(i-Nsta-1)-CF(i-Nsta-Nlta-1)Nlta---(4)]]>式中，STAi和LTAi分別表示微震信號在i個采樣點的短時平均值和長時平均值，CF(i)為記錄在i時刻的特征函數(shù)值，Nsta和Nlta分別表示為短時平均時窗和長時平均時窗所包含的記錄點數(shù)。AIC的拾取方法為公式(5)：A(i)＝ilg{var(x[1,i])}+(N-i-1)lg{var(x[i+1,N])}(5)式中，A(i)為要求的值，用于判斷信號到達(dá)的初至。var(x[1,i])表示信號從第一個采樣點到第i個采樣點的方差。N表示選取信號的總采樣點數(shù)。步驟(4)的正演計算縱波走時和橫波走時采用迭代法射線追蹤算法是，首先根據(jù)檢波器、震源和地層關(guān)系，假定初始射線路徑；然后利用費馬原理推導(dǎo)出線性方程組的系數(shù)值，設(shè)定迭代終止條件，最后通過求解線性方程組，不斷迭代計算節(jié)點，直到滿足終止條件。本發(fā)明克服了單獨利用縱波進(jìn)行水力壓裂微地震震源定位時，橫向定位精度受初至拾取誤差影響比較大的影響，利用了縱橫波走時聯(lián)合約束的反演方法進(jìn)行微震定位，提高了反演的精度。另外，在反演的方法上，本發(fā)明采用全局搜索算法和局部搜索算法相結(jié)合的方法，避免了單獨利用全局搜索算法計算量大和局部搜索算法容易陷入局部極值的影響，兼顧了定位精度和計算效率。無論是全局搜索還是局部反演算法，采用的都是完全非線性反演，特別適合微震定位這種非線性問題。附圖說明圖1為水力壓裂微震源定位技術(shù)流程圖；圖2為正演縱波加入[-11]ms誤差擾動下縱波目標(biāo)函數(shù)等值線圖；圖3為正演縱橫波加入[-11]ms誤差擾動下縱橫波目標(biāo)函數(shù)等值線圖；圖4為某實際工區(qū)中測井資料和射孔記錄約束建立的速度模型；圖5為某記錄壓噪和有效事件識別結(jié)果；圖6縱橫波初至拾取結(jié)果；圖7反演結(jié)果走時曲線和拾取走時曲線對比。具體實施方式下面從理論和實際應(yīng)用兩個方面具體說明應(yīng)用步驟和效果。首先結(jié)合附圖1，對本發(fā)明的流程、參數(shù)等做進(jìn)一步說明：(1)建立微震監(jiān)測的觀測系統(tǒng)，獲取測井資料和其他地震資料，建立初始速度模型；(2)調(diào)研分析研究區(qū)域地層條件和射孔記錄，確定背景干擾和微震信號特點，然后對壓裂過程中的記錄進(jìn)行壓噪處理以及有效微震信號進(jìn)行波形識別；(3)利用初至自動拾取和人工干預(yù)等方法對直達(dá)波的初至走時進(jìn)行拾取，其中射孔記錄的走時用于速度模型的優(yōu)化，其他有效事件用于震源的反演定位；(4)在利用縱波走時的基礎(chǔ)上，加入橫波走時，利用縱波走時和橫波走時同時約束進(jìn)行微震震源定位，反演采用全局搜索和局部搜索相結(jié)合的方法，包括:①根據(jù)水力壓裂段的位置坐標(biāo)，確定震源定位搜索范圍，提高反演效率的同時減少解的多解性；②輸入檢波器坐標(biāo)和步驟(3)中拾取到的縱橫波走時信息，在搜索范圍內(nèi)，隨機(jī)產(chǎn)生一個點作為反演初始震源點，并正演計算縱波走時和橫波走時，公式如下：T=Σi=1m+1li/vi=Σi=1m+1[(xi-xi-1)2+(yi-yi-1)2+(zi-zi-1)2]1/2/vi---(1)]]>其中m表示檢波器(xm+1,ym+1,zm+1)和假設(shè)震源點(x0,y0,z0)與建立地層模型交點個數(shù)，(xi,yi,zi)(i～m)為上述交點的空間坐標(biāo)，vi表示射線段對應(yīng)的地層橫波或者縱波速度；③建立并計算目標(biāo)函數(shù)：F=Σi=1n((Tscal(i)-Tpcal(i))-(Tsobs(i)-Tpobs(i)))2---(2)]]>其中，Tscal(i)、Tpcal(i)分別為根據(jù)已知模型和②中假設(shè)點正演計算得到的橫波走時和縱波走時，Tsobs(i)、Tpobs(i)分別為通過步驟(3)獲取到的各檢波器對應(yīng)的橫波和縱波時間值，n為檢波器個數(shù)。在目標(biāo)函數(shù)中加入橫波信息可以提高定位結(jié)果的精度；④然后判斷根據(jù)公式(2)計算的目標(biāo)函數(shù)值是否滿足預(yù)先給定的精度，如果滿足，則轉(zhuǎn)到⑦，否則，進(jìn)入下一步；⑤設(shè)定模擬退火搜索迭代終止條件，采用能找到全局極小值的模擬退火法進(jìn)行全局搜索，利用公式(1)計算走時，利用式(2)計算目標(biāo)函數(shù)，每次迭代過程中都判斷是否滿足模擬退火終止條件，如果滿足，則進(jìn)入⑦，否則，進(jìn)入下一步；⑥設(shè)定模式搜索算法的迭代終止條件，以模擬退火的輸出結(jié)果為初始輸入模型，然后用模式搜索算法進(jìn)行局部搜索，采用基于正交基的模式搜索法，模式為：Pk＝BCk(3)其中，k表示迭代次數(shù)，B為模式的基矩陣，為一非奇異矩陣，Ck為一個生成矩陣。通過當(dāng)前點和模式計算新的搜索點的目標(biāo)函數(shù)值，如果新搜索值的目標(biāo)函數(shù)小于前面所有點的目標(biāo)函數(shù)，則通過(3)式重新修改模式，繼續(xù)新的搜索，直到進(jìn)行完所有的模式搜索，最后判斷定位結(jié)果精度是否滿足，如果滿足，則進(jìn)入⑦，否則，返回到步驟⑤中，重新進(jìn)行全局搜索；⑦定位結(jié)果是否滿足實際地質(zhì)條件，如果滿足，則輸出微震震源位置。所述的測井資料為利用已有的常規(guī)聲波測井或全波列測井得到聲波資料。所述的射孔記錄為采用特殊聚能器材進(jìn)入井眼預(yù)定層位進(jìn)行爆炸開孔讓井下地層內(nèi)流體進(jìn)入孔眼時儀器記錄到的振幅強(qiáng)的微震信號。所述的壓噪算法為帶通濾波算法，所述的事件識別為長短時窗能量比法，在經(jīng)過背景噪音壓制處理后，能夠提高原始記錄的信噪比，更加有利于微震有效事件的識別。所述的初至自動拾取為長短時窗比和基于信息準(zhǔn)則(AIC)相結(jié)合的方法，其中，長短時窗比的方法為公式(4)，AIC的拾取方法為公式(5)：STAi=STAi-1+CF(i)-CF(i-Nsta)NstaLTAi=LTAi-1+CF(i-Nsta-1)-CF(i-Nsta-Nlta-1)Nlta---(4)]]>其中，STAi和LTAi分別表示微震信號在i個采樣點的短時平均值和長時平均值，CF(i)為記錄在i時刻的特征函數(shù)值，Nsta和Nlta分別表示為短時平均時窗和長時平均時窗所包含的記錄點數(shù)。A(i)＝ilg{var(x[1,i])}+(N-i-1)lg{var(x[i+1,N])}(5)其中，A(i)為要求的值，用于判斷信號到達(dá)的初至。var(x[1,i])表示信號從第一個采樣點到第i個采樣點的方差。N表示選取信號的總采樣點數(shù)。所述的橫波速度從全波列測井、偶極橫波測井或者研究區(qū)域經(jīng)驗公式獲得。所述的正演走時采用迭代法射線追蹤算法：首先根據(jù)檢波器、震源和地層關(guān)系，假定初始射線路徑；然后利用費馬原理推導(dǎo)出線性方程組的系數(shù)值；設(shè)定迭代終止條件，最后通過求解線性方程組，不斷迭代計算節(jié)點，直到滿足終止條件?；谏鲜龇桨傅幕A(chǔ)上，本次選用的理論模型中，檢波器橫坐標(biāo)都為0m，深度從-2510m至-2660m，間隔為10m，共16級檢波器，理論震源的位置為[300,-2722]m。當(dāng)單純利用縱波反演時，反演結(jié)果受初至影響較大。圖2為正演縱波走時加入[-11]ms誤差擾動下縱波目標(biāo)函數(shù)等值線圖。從圖2中可以比較直觀看到，目標(biāo)函數(shù)的最小點基本在[350,-2750]m附近，圖中用星形符號表示，然后利用所述方法反演方法但是沒有加入橫波信息時，通過P波走時單獨反演時，結(jié)果為[348,-2743]m。圖3為同時利用縱橫波走時正演計算得到的目標(biāo)函數(shù)值，同樣對縱橫波加入[-11]ms的走時拾取誤差，此時，就可以完全按照所述方法進(jìn)行縱橫波聯(lián)合反演，反演的目標(biāo)函數(shù)定義為：Σi=1n((tscal(i)-tpcal(i))-(tsobs(i)-tpobs(i)))2---(2)]]>其中tscal(i)和tpcal(i)分別為正演計算得到的縱橫波走時，tsobs(i)和tpobs(i)分別為觀測到的縱橫波走時。利用縱橫波走時聯(lián)合反演，結(jié)果為[299,-2728]m，在圖3中用星形符號表示，可以看到，在相同拾取誤差下，反演結(jié)果受初至影響較小，接近理論值，本發(fā)明的方法定位精度高。具體實施實例為使本發(fā)明的目的、技術(shù)方案更加清楚明白，下面結(jié)合附圖和實例對本發(fā)明進(jìn)行詳細(xì)說明，圖1為本發(fā)明的技術(shù)流程圖。主要步驟如下：(1)建立微震監(jiān)測的觀測系統(tǒng)，獲取測井資料和其他地震資料，建立初始速度模型；(2)調(diào)研分析研究區(qū)域地層條件和射孔記錄，確定背景干擾和微震信號特點，然后對壓裂過程中的記錄進(jìn)行壓噪處理以及有效微震信號進(jìn)行波形識別。根據(jù)該地區(qū)的微震信號分析和濾波參數(shù)實驗后，發(fā)現(xiàn)40-400Hz的帶通濾波參數(shù)能夠去除低頻干擾影響，保留有效的微震信號；(3)利用初至自動拾取和人工干預(yù)等方法對直達(dá)波的初至走時進(jìn)行拾取，其中射孔記錄的走時用于速度模型的優(yōu)化，其他有效事件用于震源的反演定位，圖4為利用射孔記錄和測井資料約束共同建立的速度模型，圖5為實際地震數(shù)據(jù)中通過壓噪處理和波形識別后的某個有效微震記錄，可以看到，該記錄信噪比高，有利于初至的準(zhǔn)確拾取；(4)在利用縱波走時的基礎(chǔ)上，加入橫波走時，利用縱波走時和橫波走時同時約束進(jìn)行微震震源定位，反演采用全局搜索和局部搜索相結(jié)合的方法:①根據(jù)水力壓裂段的位置坐標(biāo)，確定震源定位搜索范圍，提高反演效率的同時減少解的多解性，本實例中，根據(jù)監(jiān)測范圍，橫向上的坐標(biāo)范圍為[150600]m，縱向上的深度范圍為[26502750]m；②輸入檢波器坐標(biāo)和步驟(3)中拾取到的縱橫波走時信息，其中走時信息如圖6所示，在搜索范圍內(nèi)，隨機(jī)產(chǎn)生一個點作為反演初始震源點，并正演計算縱波走時和橫波走時，公式如下：T=Σi=1m+1li/vi=Σi=1m+1[(xi-xi-1)2+(yi-yi-1)2+(zi-zi-1)2]1/2/vi---(1)]]>其中m表示檢波器(xm+1,ym+1,zm+1)和假設(shè)震源點(x0,y0,z0)與建立地層模型交點個數(shù)，(xi,yi,zi)(i～m)為上述交點的空間坐標(biāo)，vi表示射線段對應(yīng)的地層橫波或者縱波速度；③建立并計算目標(biāo)函數(shù)：F=Σi=1n((Tscal(i)-Tpcal(i))-(Tsobs(i)-Tpobs(i)))2---(2)]]>其中，Tscal(i)、Tpcal(i)分別為根據(jù)已知模型和②中假設(shè)點正演計算得到的橫波走時和縱波走時，Tsobs(i)、Tpobs(i)分別為通過步驟(3)獲取到的各檢波器對應(yīng)的時間值，n為檢波器個數(shù)。在目標(biāo)函數(shù)中加入橫波信息有兩方面作用，提高定位結(jié)果的精度和減少反演發(fā)震時刻信息；④然后判斷根據(jù)公式(2)計算的目標(biāo)函數(shù)值是否滿足預(yù)先給定的精度，如果滿足，則轉(zhuǎn)到⑦，否則，進(jìn)入下一步；⑤設(shè)定模擬退火搜索迭代終止條件，采用能找到全局極小值的模擬退火法進(jìn)行全局搜索，利用公式(1)計算走時，利用式(2)計算目標(biāo)函數(shù)，每次迭代過程中都判斷是否滿足模擬退火終止條件，如果滿足，則進(jìn)入⑦，否則，進(jìn)入下一步；⑥設(shè)定模式搜索算法的迭代終止條件，以模擬退火的輸出結(jié)果為初始輸入模型，然后用模式搜索算法進(jìn)行局部搜索，采用基于正交基的模式搜索法，模式為：Pk＝BCk(3)其中，k表示迭代次數(shù)，B為模式的基矩陣，為一非奇異矩陣，Ck為一個生成矩陣。通過當(dāng)前點和模式計算新的搜索點的目標(biāo)函數(shù)值，如果新搜索值的目標(biāo)函數(shù)小于前面所有點的目標(biāo)函數(shù)，則通過(3)式重新修改模式，繼續(xù)新的搜索，直到進(jìn)行完所有的模式搜索，最后判斷定位結(jié)果精度是否滿足，如果滿足，則進(jìn)入⑦，否則，返回到步驟⑤中，重新進(jìn)行全局搜索，在本次實例中，采用基于最大正交基的模式搜索方法，最大搜索次數(shù)設(shè)為1000次；⑦定位結(jié)果是否滿足實際地質(zhì)條件，如果滿足，則輸出微震震源位置。所述的測井資料為利用已有的常規(guī)聲波測井或全波列測井得到聲波資料。所述的射孔記錄為采用特殊聚能器材進(jìn)入井眼預(yù)定層位進(jìn)行爆炸開孔讓井下地層內(nèi)流體進(jìn)入孔眼時儀器記錄到的振幅強(qiáng)的微震信號。所述的壓噪算法為帶通濾波算法，所述的事件識別為長短時窗能量比法，在經(jīng)過背景噪音壓制處理后，能夠提高原始記錄的信噪比，更加有利于微震有效事件的識別。所述的初至自動拾取為長短時窗比和基于信息準(zhǔn)則(AIC)相結(jié)合的方法，其中，長短時窗比的方法為公式(4)，AIC的拾取方法為公式(5)：STAi=STAi-1+CF(i)-CF(i-Nsta)NstaLTAi=LTAi-1+CF(i-Nsta-1)-CF(i-Nsta-Nlta-1)Nlta---(4)]]>其中，STAi和LTAi分別表示微震信號在i個采樣點的短時平均值和長時平均值，CF(i)為記錄在i時刻的特征函數(shù)值，Nsta和Nlta分別表示為短時平均時窗和長時平均時窗所包含的記錄點數(shù)。A(i)＝ilg{var(x[1,i])}+(N-i-1)lg{var(x[i+1,N])}(5)其中，A(i)為要求的值，用于判斷信號到達(dá)的初至。var(x[1,i])表示信號從第一個采樣點到第i個采樣點的方差。N表示選取信號的總采樣點數(shù)。圖6為根據(jù)所述方法拾取到的微震縱橫波走時信息。所述的橫波速度從全波列測井、偶極橫波測井或者研究區(qū)域經(jīng)驗公式獲得。所述的正演走時采用迭代法射線追蹤算法：首先根據(jù)檢波器、震源和地層關(guān)系，假定初始射線路徑；然后利用費馬原理推導(dǎo)出線性方程組的系數(shù)值；設(shè)定迭代終止條件，最后通過求解線性方程組，不斷迭代計算節(jié)點，直到滿足終止條件。定位效果如圖7所述，實線部分表示利用所述方法拾取到的走時信息，其中前面的一條實斜線代表縱波走時信息，后面的一條實斜線代表橫波走時信息，虛線表示利用所述方法拾取到的走時信息，其中前面的一條虛斜線代表縱波走時信息，后面的一條虛斜線代表橫波走時信息，可以看到，二者基本重合，說明了本發(fā)明的有效性。當(dāng)前第1頁1 2 3