本發(fā)明涉及水力壓裂微地震監(jiān)測(cè)
技術(shù)領(lǐng)域：
，特別是微地震事件定位方法。
背景技術(shù)：
：油氣田的微地震監(jiān)測(cè)一種主要用于油氣田開(kāi)發(fā)的新的地震方法，它是利用水力壓裂、油氣開(kāi)采、常規(guī)注水等石油工程作業(yè)而產(chǎn)生的地震波。常規(guī)的井中微震事件定位前提是，具有較高的信噪比，微地震事件的同向軸具有良好的連續(xù)性，在此前提下利用微地震事件初至的匹配程度作為目標(biāo)函數(shù)進(jìn)行定位。這存在兩個(gè)問(wèn)題：1)初至拾取的人為誤差問(wèn)題2)信噪比較低時(shí)微地震事件同向軸沒(méi)有良好的連續(xù)性是無(wú)法拾取初至。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是針對(duì)現(xiàn)有技術(shù)存在的問(wèn)題，提供一種微地震事件定位方法。該方法利用任意兩道間的波形互相關(guān)，計(jì)算拾取初至?xí)r差，同時(shí)利用所有檢波器對(duì)間的拾取初至差與計(jì)算初至差間的差的加權(quán)疊加作為目標(biāo)函數(shù),對(duì)微地震事件進(jìn)行優(yōu)化迭代定位。本發(fā)明通過(guò)如下方法措施來(lái)實(shí)現(xiàn)：一種微地震事件定位方法,首先確定工區(qū)的層狀速度模型和檢波器的位置,確定微地震事件信號(hào)，利用所有檢波器對(duì)間的拾取初至差與計(jì)算初至差間的差的加權(quán)疊加作為目標(biāo)函數(shù),對(duì)微地震事件進(jìn)行優(yōu)化迭代定位。上述方案中優(yōu)化迭代定位過(guò)程包含如下步驟：1)確定初始的微地震位置p0(x0,y0,z0)，和迭代終止閾值dth；2)在可能的微地震事件位置范圍內(nèi)應(yīng)用一定網(wǎng)格劃分空間；3)根據(jù)初始位置p0進(jìn)行優(yōu)化定位，其定位結(jié)果為p1(x1,y1,z1)；4)計(jì)算初始位置p0與優(yōu)化位置p1之間距離d01，其計(jì)算公式如下d01＝[(x1-x0)2+(y1-y0)2+(z1-z0)2]1/2公式(1)；5)如果d01>dth,令p0＝p1，轉(zhuǎn)到步驟3)；否則結(jié)束迭代微地震事件位置pe＝p1；上述方案進(jìn)一步包括：步驟3)中的優(yōu)化定位過(guò)程包含以下步驟：以初始的微地震位置p0和劃分的網(wǎng)格點(diǎn)為基礎(chǔ)，應(yīng)用網(wǎng)格搜索算法針對(duì)每一個(gè)空間網(wǎng)格點(diǎn)(xs,ys,zs)計(jì)算目標(biāo)函數(shù)；選取目標(biāo)函數(shù)最小的網(wǎng)格點(diǎn)的位置作為定位結(jié)果p1。進(jìn)一步，假設(shè)包含N級(jí)檢波器的觀(guān)測(cè)系統(tǒng)，任意兩級(jí)檢波器對(duì)為Ri,Rj，所述的空間網(wǎng)格點(diǎn)(xs,ys,zs)的目標(biāo)函數(shù)建立包含如下步驟：計(jì)算空間網(wǎng)格點(diǎn)(xs,ys,zs)對(duì)應(yīng)的兩級(jí)檢波器間的拾取初至的差為Δtpij和計(jì)算初至的差為Δtcsij；網(wǎng)格點(diǎn)(xs,ys,zs)的目標(biāo)函數(shù)應(yīng)用下式計(jì)算F(xs,ys,zs)=Σi=1NΣj>iNλij|Δtpij-Δtcsij|]]>公式(2)其中λij為檢波器對(duì)Ri,Rj對(duì)應(yīng)的權(quán)值。進(jìn)一步，檢波器對(duì)Ri,Rj間的拾取初至的差Δtpij和互相關(guān)系數(shù)rij應(yīng)用兩道間的微地震記錄進(jìn)行互相關(guān)運(yùn)算的方法求取，其計(jì)算步驟如下：分別選取檢波器對(duì)Ri,Rj的微地震信號(hào)分別為Si(h1)和Sj(h2)，其中1≤h1≤n1，1≤h2≤n2采樣點(diǎn)Si(1)和Sj(1)對(duì)應(yīng)的時(shí)刻分別為tpi和tpj，設(shè)微地震信號(hào)的采樣間隔為Δt；分別將信號(hào)Si(h1)和Sj(h2)延拓為SI(h3)和SJ(h4)，其中1≤h3≤n1+n2-1，1≤h4≤n1+n2-1，且SI(h3)＝Si(h3)1≤h3≤n1公式(3)SI(h3)＝0n1+1≤h3≤n1+n2-1SJ(h4)＝01≤h4≤n1-1公式(4)SJ(h4)＝Si(h4)n1≤h4≤n1+n2-1對(duì)Si(h1)和Sj(h2)進(jìn)行互相關(guān)運(yùn)算用下式表示XRij(τ)=Σk=1τSI(k)SI(n1+n2-τ)1≤τ≤n1+n2-1]]>公式(5)考慮到微地震信號(hào)Si(h1)和Sj(h2)極性可能相反，對(duì)互相關(guān)函數(shù)取絕對(duì)值即：XRij(τ)＝|XRij(τ)|公式(6)設(shè)τ＝τmax時(shí)XRij(τ)取得最大值，則Δtpij按照下式計(jì)算Δtpij＝tpj-tpi+Δt*(τmax-n1-n2+1)j≥i公式(7)Δtpij＝tpi-tpj-Δt*(τmax-n1-n2+1)j<i兩道間的最大互相關(guān)系數(shù)rij按照下式計(jì)算rij=XRij(τmax)/Σk=1τmaxSI(k)2Σk=1τmaxSJ(n1+n2-τ)2]]>公式(8)其中互相關(guān)系數(shù)rij代表檢波器對(duì)Ri,Rj間的拾取初至的差Δtpij的可信程度。空間網(wǎng)格點(diǎn)(xs,ys,zs)對(duì)應(yīng)檢波器對(duì)Ri,Rj間的計(jì)算初至的差Δtcij應(yīng)用以下步驟實(shí)現(xiàn)：通過(guò)射線(xiàn)追蹤方法空間網(wǎng)格點(diǎn)位置(xs,ys,zs)到對(duì)應(yīng)兩級(jí)檢波器對(duì)Ri,Rj的初至旅行時(shí)tcsi和tcsj；應(yīng)用以下公式計(jì)算檢波器對(duì)間的計(jì)算初至的差Δtcsij＝tcsi-tcsjj≥i公式(9)Δtcsij＝tcsi-tcsjj<i檢波器對(duì)Ri,Rj的權(quán)重系數(shù)λij的確定過(guò)程包含如下步驟：根據(jù)確定的層狀速度模型，分別應(yīng)用射線(xiàn)追蹤方法計(jì)算微地震初始位置p0到檢波器對(duì)Ri,Rj位置的射線(xiàn)線(xiàn)路經(jīng)Rayi和Rayj；設(shè)射線(xiàn)路徑Rayi和Rayj所經(jīng)過(guò)的點(diǎn)分別為(pi1,pi2,pi3..piNi)和(pj1,pj2,pj3..pjNj)，其中pi1和pj1微地震初始位置p0，piNi和pjNj分別為檢波器對(duì)Ri,Rj所在的位置，(pi2,pi3..piNi-1)為射線(xiàn)路徑Rayi與層狀介質(zhì)的分界面的交點(diǎn)，(pj2,pj3..pjNj-1)為射線(xiàn)路徑Rayj與層狀介質(zhì)的分界面的交點(diǎn),則射線(xiàn)路徑Rayi和Rayj在p0位置處的射線(xiàn)方向用如下單位向量表示，Vi＝[vixviyviz]公式(10)Vj＝[vjxvjyvjz]其中vix＝(xi2-x0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2viy＝(yi2-y0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2公式(11)viz＝(zi2-z0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2vjx＝(xj2-x0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2vjy＝(yj2-y0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2公式(12)vjz＝(zj2-z0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2其中(x0,y0,z0)為p0所對(duì)應(yīng)的坐標(biāo)位置,(xi2，yi2，zi2)為pj2點(diǎn)所對(duì)應(yīng)的坐標(biāo)位置,(xj2，yj2，zj2)為pj2點(diǎn)所對(duì)應(yīng)的坐標(biāo)位置；檢波器對(duì)Ri,Rj對(duì)應(yīng)的射線(xiàn)線(xiàn)路徑Rayi和Rayj在p0處的射線(xiàn)差異可以用LΔVij表示，其計(jì)算公式如下LΔVij=||ΔVij||2]]>公式(13)其中ΔVij＝Vi-Vj公式(14)||ΔVij||2代表的是向量ΔVij的第二范數(shù)，射線(xiàn)差異代表定位結(jié)果對(duì)Ri,Rj間的拾取初至差與計(jì)算初至差的差的誤差的敏感程度，越大敏感程度越大，反之亦然；檢波器對(duì)Ri,Rj對(duì)應(yīng)的權(quán)值λij應(yīng)用互相關(guān)系數(shù)rij與射線(xiàn)路徑差異的模表示，用下式計(jì)算λij=rij*LΔVij]]>公式(15)。本發(fā)明利用任意兩道間的波形互相關(guān)，計(jì)算拾取初至?xí)r差，同時(shí)利用所有檢波器對(duì)間的拾取初至差與計(jì)算初至差間的差的加權(quán)疊加作為目標(biāo)函數(shù),對(duì)微地震事件進(jìn)行優(yōu)化迭代定位。該方法可以有效的克服低信噪比同向軸不連續(xù)的問(wèn)題；同時(shí)利用“利用所有檢波器對(duì)間的拾取初至差與計(jì)算初至差間的差的加權(quán)疊加作為目標(biāo)函數(shù)”，對(duì)于包含N個(gè)建波器的觀(guān)測(cè)系統(tǒng)可以提供N(N-1)/2種檢波器的組合，從而提供了大量的信息冗余，提高微地震事件的定位精度。合理的權(quán)重系數(shù)選擇對(duì)于定位結(jié)果至關(guān)重要，本發(fā)明中主要從以下兩方面考慮一是：計(jì)算拾取初至?xí)r差時(shí)的最大互相關(guān)時(shí)的相關(guān)系數(shù)，用來(lái)體現(xiàn)該時(shí)差的可信程度，相關(guān)系數(shù)越高，可信度越高；二是：利用兩級(jí)檢波器的在“微震事件點(diǎn)位置處”射線(xiàn)路徑的差異，射線(xiàn)路徑的差異可以體現(xiàn)道間時(shí)差對(duì)定位結(jié)果的敏 感程度，射線(xiàn)路徑差異越大，敏感程度越低，權(quán)重越大。附圖說(shuō)明圖1本發(fā)明處理流程圖。圖2正演觀(guān)測(cè)系統(tǒng)及定位過(guò)程和結(jié)果。其中圖2-a正演觀(guān)測(cè)系統(tǒng)及定位過(guò)程和結(jié)果；標(biāo)注1為檢波器，標(biāo)注2為網(wǎng)格掃描過(guò)過(guò)程中的第一個(gè)點(diǎn)(-200,200,2800)，標(biāo)注3為可能的微地震區(qū)域，標(biāo)注4為初始微地震位置p0(100,100,2900)。圖2-b是圖2-a震源點(diǎn)附近放大圖；標(biāo)注5為首次迭代p0為(100,100,2900)時(shí)的定位結(jié)果(0，0，2990)，標(biāo)注6為最終定位結(jié)果(0,0,2998)，標(biāo)注7為事件點(diǎn)的真實(shí)位置，標(biāo)注8為應(yīng)用常規(guī)定位方法的定位結(jié)果(-12,10,3014)。圖3正演速度模型。圖4三種不同子波。圖5正演微地震記錄及拾取的初至；標(biāo)注9為拾取的初至。圖6利用互相關(guān)法得到的任意檢波器對(duì)拾取初至的差。圖7利用互相關(guān)法計(jì)算任意檢波器對(duì)拾取初至的差時(shí)的互相關(guān)系數(shù)。圖8網(wǎng)格點(diǎn)(-200,200,2800)對(duì)應(yīng)的任意兩級(jí)檢波器間的計(jì)算初至的差。圖9相對(duì)于初始點(diǎn)位置p0計(jì)算得到的任意兩級(jí)檢波器射線(xiàn)路徑的差異。圖10相對(duì)于首次迭代初始點(diǎn)位置p0為(100,100,2900)時(shí)計(jì)算得到的任意兩級(jí)檢波器權(quán)重系數(shù)。圖11初次迭代初始點(diǎn)位置p0為(100,100,2900)時(shí)計(jì)算所有網(wǎng)格點(diǎn)的目標(biāo)函數(shù)和最小目標(biāo)函數(shù)值點(diǎn)；標(biāo)注10為目標(biāo)函數(shù)最小值點(diǎn)。具體實(shí)施方式下文特舉出較佳實(shí)施例，對(duì)本發(fā)明作詳細(xì)說(shuō)明如下。結(jié)合附圖1，本實(shí)施例的具體流程為：一種微地震事件定位方法,首先確定工區(qū)的層狀速度模型和檢波器的位置,確定微地震事件信號(hào)，利用所有檢波器對(duì)間的拾取初至差與計(jì)算初至差間的差的加權(quán)疊加作為目標(biāo)函數(shù),對(duì)微地震事件進(jìn)行優(yōu)化迭代定位。其中，優(yōu)化迭代定位過(guò)程包含如下步驟：1)確定初始的微地震位置p0(x0,y0,z0)，和迭代終止閾值dth；2)在可能的微地震事件位置范圍內(nèi)應(yīng)用一定網(wǎng)格劃分空間；3)根據(jù)初始位置p0進(jìn)行優(yōu)化定位，其定位結(jié)果為p1(x1,y1,z1)；4)計(jì)算初始位置p0與優(yōu)化位置p1之間距離d01，其計(jì)算公式如下d01＝[(x1-x0)2+(y1-y0)2+(z1-z0)2]1/2公式(1)；5)如果d01>dth,令p0＝p1，轉(zhuǎn)到步驟3)；否則結(jié)束迭代微地震事件位置pe＝p1；步驟3)中的優(yōu)化定位過(guò)程包含以下步驟：以初始的微地震位置p0和劃分的網(wǎng)格點(diǎn)為基礎(chǔ)，應(yīng)用網(wǎng)格搜索算法針對(duì)每一個(gè)空間網(wǎng)格點(diǎn)(xs,ys,zs)計(jì)算目標(biāo)函數(shù)；選取目標(biāo)函數(shù)最小的網(wǎng)格點(diǎn)的位置作為定位結(jié)果p1。進(jìn)一步，假設(shè)包含N級(jí)檢波器的觀(guān)測(cè)系統(tǒng)，任意兩級(jí)檢波器對(duì)為Ri,Rj，所述的空間網(wǎng)格點(diǎn)(xs,ys,zs)的目標(biāo)函數(shù)建立包含如下步驟：計(jì)算空間網(wǎng)格點(diǎn)(xs,ys,zs)對(duì)應(yīng)的兩級(jí)檢波器間的拾取初至的差為Δtpij和計(jì)算初至的差為Δtcsij；網(wǎng)格點(diǎn)(xs,ys,zs)的目標(biāo)函數(shù)應(yīng)用下式計(jì)算F(xs,ys,zs)=Σi=1NΣj<iNλij|Δtpij-Δtcsij|]]>公式(2)其中λij為檢波器對(duì)Ri,Rj對(duì)應(yīng)的權(quán)值。進(jìn)一步，檢波器對(duì)Ri,Rj間的拾取初至的差Δtpij和互相關(guān)系數(shù)rij應(yīng)用兩道間的微地震記錄進(jìn)行互相關(guān)運(yùn)算的方法求取，其計(jì)算步驟如下：分別選取檢波器對(duì)Ri,Rj的微地震信號(hào)分別為Si(h1)和Sj(h2)，其中1≤h1≤n1，1≤h2≤n2采樣點(diǎn)Si(1)和Sj(1)對(duì)應(yīng)的時(shí)刻分別為tpi和tpj，設(shè)微地震信號(hào)的采樣間隔為Δt；分別將信號(hào)Si(h1)和Sj(h2)延拓為SI(h3)和SJ(h4)，其中1≤h3≤n1+n2-1，1≤h4≤n1+n2-1，且SI(h3)＝Si(h3)1≤h3≤n1公式(3)SI(h3)＝0n1+1≤h3≤n1+n2-1SJ(h4)＝01≤h4≤n1-1公式(4)SJ(h4)＝Si(h4)n1≤h4≤n1+n2-1對(duì)Si(h1)和Sj(h2)進(jìn)行互相關(guān)運(yùn)算用下式表示XRij(τ)=Σk=1τSI(k)SI(n1+n2-τ)1≤τ≤n1+n2-1]]>公式(5)考慮到微地震信號(hào)Si(h1)和Sj(h2)極性可能相反，對(duì)互相關(guān)函數(shù)取絕對(duì)值即：XRij(τ)＝|XRij(τ)|公式(6)設(shè)τ＝τmax時(shí)XRij(τ)取得最大值，則Δtpij按照下式計(jì)算Δtpij＝tpj-tpi+Δt*(τmax-n1-n2+1)j≥i公式(7)Δtpij＝tpi-tpj-Δt*(τmax-n1-n2+1)j<i兩道間的最大互相關(guān)系數(shù)rij按照下式計(jì)算rij=XRij(τmax)/Σk=1τmaxSI(k)2Σk=1τmaxSJ(n1+n2-τ)2]]>公式(8)其中互相關(guān)系數(shù)rij代表檢波器對(duì)Ri,Rj間的拾取初至的差Δtpij的可信程度?？臻g網(wǎng)格點(diǎn)(xs,ys,zs)對(duì)應(yīng)檢波器對(duì)Ri,Rj間的計(jì)算初至的差Δtcij應(yīng)用以下步驟實(shí)現(xiàn)：通過(guò)射線(xiàn)追蹤方法空間網(wǎng)格點(diǎn)位置(xs,ys,zs)到對(duì)應(yīng)兩級(jí)檢波器對(duì)Ri,Rj的初至旅行時(shí)tcsi和tcsj；應(yīng)用以下公式計(jì)算檢波器對(duì)間的計(jì)算初至的差Δtcsij＝tcsi-tcsjj≥i公式(9)Δtcsij＝tcsi-tcsjj<i檢波器對(duì)Ri,Rj的權(quán)重系數(shù)λij的確定過(guò)程包含如下步驟：根據(jù)確定的層狀速度模型，分別應(yīng)用射線(xiàn)追蹤方法計(jì)算微地震初始位置p0到檢波器對(duì)Ri,Rj位置的射線(xiàn)線(xiàn)路經(jīng)Rayi和Rayj；設(shè)射線(xiàn)路徑Rayi和Rayj所經(jīng)過(guò)的點(diǎn)分別為(pi1,pi2,pi3..piNi)和(pj1,pj2,pj3..pjNj)，其中pi1和pj1微地震初始位置p0，piNi和pjNj分別為檢波器對(duì)Ri,Rj所在的位置，(pi2,pi3..piNi-1)為射線(xiàn)路徑Rayi與層狀介質(zhì)的分界面的交點(diǎn)，(pj2,pj3..pjNj-1)為射線(xiàn)路徑Rayj與層狀介質(zhì)的分界面的交點(diǎn),則射線(xiàn)路徑Rayi和Rayj在p0位置處的射線(xiàn)方向用如下單位向量表示，Vi＝[vixviyviz]公式(10)Vj＝[vjxvjyvjz]其中vix＝(xi2-x0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2viy＝(yi2-y0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2公式(11)viz＝(zi2-z0)/[(xi2-x0)2+(yi2-y0)2+(zi2-z0)2]1/2vjx＝(xj2-x0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2vjy＝(yj2-y0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2公式(12)vjz＝(zj2-z0)/[(xj2-x0)2+(yj2-y0)2+(zj2-z0)2]1/2其中(x0,y0,z0)為p0所對(duì)應(yīng)的坐標(biāo)位置,(xi2，yi2，zi2)為pj2點(diǎn)所對(duì)應(yīng)的坐標(biāo)位置,(xj2，yj2，zj2)為pj2點(diǎn)所對(duì)應(yīng)的坐標(biāo)位置；檢波器對(duì)Ri,Rj對(duì)應(yīng)的射線(xiàn)線(xiàn)路徑Rayi和Rayj在p0處的射線(xiàn)差異可以用LΔVij表示，其計(jì)算公式如下LΔVij=||ΔVij||2]]>公式(13)其中ΔVij＝Vi-Vj公式(14)||ΔVij||2代表的是向量ΔVij的第二范數(shù)，射線(xiàn)差異代表定位結(jié)果對(duì)Ri,Rj間的拾取初至差與計(jì)算初至差的差的誤差的敏感程度，越大敏感程度越大，反之亦然；檢波器對(duì)Ri,Rj對(duì)應(yīng)的權(quán)值λij應(yīng)用互相關(guān)系數(shù)rij與射線(xiàn)路徑差異的模表示，用下式計(jì)算λij=rij*LΔVij]]>公式(15)。下面結(jié)合附圖2-11，對(duì)上述實(shí)施例的實(shí)施方案做進(jìn)一步說(shuō)明。建立如圖2所示觀(guān)測(cè)系統(tǒng)，應(yīng)用如圖3所示的速度模型，正演模擬過(guò)程中每級(jí)檢波器隨機(jī)選用不同的震源子波如圖4所示，同時(shí)添加按照信噪比為1添加隨機(jī)噪音最終獲得如圖5所示的微地震記錄；利用已知的速度模型(圖3所示)和檢波器位置(圖2-a所示)對(duì)該微震記錄進(jìn)行重新定位1)以偏離微震事件點(diǎn)173米的位置(100,100,2900)為初始位置p0(如圖2-a標(biāo)注4)，設(shè)置迭代終止閾值dth＝3米；2)以震源位置為中心將其周?chē)?00*200*200三維空間(如圖2-a標(biāo)注3)作為可能的微地震事件范圍，以2*2*2米的網(wǎng)格進(jìn)行劃分；3)利用化相關(guān)法計(jì)算出任意檢波器對(duì)Ri,Rj間的拾取初至的差Δtpij(如圖6所示)，同時(shí)計(jì)算出任意檢波器對(duì)互相關(guān)系數(shù)rij(如圖7所示)；4)根據(jù)射線(xiàn)追蹤方法，利用公式9計(jì)出步驟1)劃分網(wǎng)格中任意一點(diǎn)(xs,ys,zs)相對(duì)于檢波器對(duì)Ri,Rj間的計(jì)算初至的差Δtcsij(如圖8所示為網(wǎng)格掃描過(guò)程中第一個(gè)網(wǎng)格點(diǎn)(-200,-200,2800)(如圖2-a標(biāo)注2)對(duì)應(yīng)的所有檢波器對(duì)計(jì)算初至的差)；5)利用公式13計(jì)算任意檢波器對(duì)Ri,Rj在p0位置處的射線(xiàn)路徑差的差異(如圖9所示為首次迭代p0為(100,100,2900)時(shí)，所有檢波器對(duì)在p0位置處的射線(xiàn)路徑差的差異)；6)根據(jù)公式15，利用步驟3)中計(jì)算的任意檢波器對(duì)互相關(guān)系數(shù)和步驟5)中計(jì)算的計(jì)算任意檢波器對(duì)Ri,Rj在p0位置處的射線(xiàn)路徑差的差異計(jì)算任意檢波器對(duì)的權(quán)重系數(shù)λij(如圖10所示為首次迭代p0為(100,100,2900)時(shí)，任意檢波器對(duì)Ri,Rj的權(quán)重系數(shù))。7)根據(jù)公式2，利用步驟3)中計(jì)算出的任意檢波器對(duì)Ri,Rj間的 拾取初至的差Δtpij和步驟4)中計(jì)算出的劃分網(wǎng)格中任意一點(diǎn)相對(duì)于檢波器對(duì)Ri,Rj間的計(jì)算初至的差Δtcsij以及步驟6)中計(jì)算得到的任意檢波器對(duì)Ri,Rj的權(quán)重系數(shù)λij，計(jì)算出所有掃描網(wǎng)格點(diǎn)的目標(biāo)函數(shù)值(如圖11所示為首次迭代初始點(diǎn)位置p0為(100,100,2900)時(shí)所有網(wǎng)格點(diǎn)的目標(biāo)函數(shù)值)8)選取步驟7中目標(biāo)函數(shù)值中最小的點(diǎn)對(duì)應(yīng)的位置為當(dāng)前迭代定位結(jié)果p1。(如圖2-b所示標(biāo)注5為首次迭代p0為(100,100,2900)時(shí)的定位結(jié)果(0，0，2990))；9)根據(jù)公式1計(jì)算初始位置p0與當(dāng)前迭定位結(jié)果p1的距離d01(首次迭代p0為(100,100,2900)時(shí)d01＝10米)；10)如果d01>dth(首次迭代d01(10米)>dth(3米))，則令p0＝p1轉(zhuǎn)到步驟5)繼續(xù)計(jì)算；否則結(jié)束迭代過(guò)程，微震事件的定位結(jié)果pms＝p1(本次定位共迭代3次，最終定位結(jié)果為pms為(0,0,2998)，如圖2-b中標(biāo)注6所示)；同時(shí)根據(jù)微地震事件拾取的初至(如圖5標(biāo)注9所示)根據(jù)常規(guī)的震源定位方法(以計(jì)算初至和拾取初至的匹配程度作為目標(biāo)函數(shù)利用網(wǎng)格搜索算法進(jìn)行定位)計(jì)算得到的定位結(jié)果為(-12,10,3014)(如圖2-b標(biāo)注8所示)，本發(fā)明方法能夠很好的克服初至拾取困難造成的定位不準(zhǔn)確問(wèn)題，提高了定位精度。當(dāng)前第1頁(yè)1 2 3