專利名稱：基于方位角約束的微地震事件定位方法
技術(shù)領(lǐng)域：
本發(fā)明涉及石油天然氣地震勘探領(lǐng)域的微地震震源定位技術(shù)，更具體地，涉及一種基于方位角約束掃描疊加的微地震事件定位方法。
背景技術(shù)：
由于需要比較快速地得到壓裂結(jié)果，對(duì)壓裂效果進(jìn)行評(píng)價(jià)，進(jìn)而對(duì)開采方案進(jìn)行及時(shí)調(diào)整，故需要對(duì)水力壓裂微地震監(jiān)測(cè)資料進(jìn)行實(shí)時(shí)的處理解釋。目前所采用的微地震監(jiān)測(cè)資料處理方法主要以人工識(shí)別、初至反演為基礎(chǔ)，由于微地震監(jiān)測(cè)資料數(shù)據(jù)記錄時(shí)間長、數(shù)據(jù)量大，采用人工識(shí)別并對(duì)微地震有效事件進(jìn)行初至拾取會(huì)耗費(fèi)大量的人力與時(shí)間。傳統(tǒng)方法處理一次水力壓裂監(jiān)測(cè)數(shù)據(jù)得到裂縫分布圖像需要比較長的時(shí)間，不能滿足實(shí)時(shí)監(jiān)測(cè)的需求；另一方面，傳統(tǒng)方法無法適用于開發(fā)、監(jiān)測(cè)等今后微地震監(jiān)測(cè)的發(fā)展方向。盡管實(shí)時(shí)監(jiān)測(cè)的方法國內(nèi)外已有人研究(主要是疊加方法)，但是目前研究的疊加方法精度較低、多解問題嚴(yán)重，不能較好地應(yīng)用到實(shí)際微地震實(shí)時(shí)監(jiān)測(cè)中。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種基于方位角約束的微地震事件定位方法，所述方法包括根據(jù)測(cè)井資料及射孔資料在井口下包含微地震源的地下空間中建立速度模型；對(duì)井口下包含微地震源的地下空間建立三維坐標(biāo)系統(tǒng)并進(jìn)行離散網(wǎng)格化，并通過速度模型來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差；從實(shí)際采集的微地震資料中識(shí)別微地震事件；對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離，以獲得微地震事件的相應(yīng)震相的走時(shí)和走時(shí)差；通過使用微地震事件的相應(yīng)震相的走時(shí)和走時(shí)差以及每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差對(duì)微地震事件進(jìn)行掃描疊加，以獲得微地震事件所對(duì)應(yīng)的至少一個(gè)網(wǎng)格點(diǎn)；使用微地震事件所對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)和微地震事件在三維坐標(biāo)系統(tǒng)的三個(gè)分量中的兩個(gè)水平分量來計(jì)算微地震事件的方位角，并根據(jù)微地震事件的方位角以及微地震事件對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)來確定微地震事件發(fā)生的位置。對(duì)井口下包含微地震源的地下空間建立速度模型的步驟可包括根據(jù)測(cè)井資料對(duì)井口下包含微地震源的地下空間劃分層位，建立初始水平層狀速度模型；通過使用射孔記錄的走時(shí)差以及已知的射孔點(diǎn)位置，根據(jù)初始水平層狀速度模型，通過反演對(duì)初始水平層狀速度模型進(jìn)行修正以獲得精細(xì)速度模型?？赏ㄟ^使用三分量檢波器來實(shí)際采集微地震資料。對(duì)井口下包含微地震源的地下空間建立三維坐標(biāo)系統(tǒng)的步驟可包括以井口或第一級(jí)三分量級(jí)檢波器作為坐標(biāo)原點(diǎn)、正東為X方向、正北為Y方向、垂直向下為Z方向來建立三維坐標(biāo)系統(tǒng)。計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差的步驟可包括使用射線追蹤方法來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差。從實(shí)際采集的微地震資料中識(shí)別微地震事件的步驟可包括對(duì)微地震資料進(jìn)行如下操作對(duì)微地震資料設(shè)置預(yù)定門限值；選擇微地震資料中的較強(qiáng)的三道數(shù)據(jù)以及合適的窗口長度；沿選擇的三道數(shù)據(jù)在所述窗口內(nèi)進(jìn)行能量大小計(jì)算和偏振分析，以獲得并判斷窗口的前時(shí)窗中的數(shù)據(jù)的能量與窗口的后時(shí)窗中的數(shù)據(jù)的能量的比值是否達(dá)到預(yù)定門限值；如果達(dá)到預(yù)定門限值，則從微地震資料中的較強(qiáng)的三道數(shù)據(jù)中拾取出微地震事件，并對(duì)拾取出的微地震事件截取合適大小，如果未達(dá)到預(yù)定門限值，則進(jìn)行窗口滑動(dòng)并再次進(jìn)行能量大小計(jì)算和偏振分析，直到對(duì)三道數(shù)據(jù)完成能量大小計(jì)算和偏振分析以及微地震事件的拾取。對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離的步驟可包括當(dāng)微地震事件是雙軸掃描事件時(shí)，如果第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差大于第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差，則第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)S波，第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)P波，反之亦然。對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離的步驟可包括當(dāng)微地震事件是單軸掃描 事件時(shí)，確定與被掃描的微地震事件相鄰的微地震事件的走時(shí)和位置，當(dāng)被掃描的微地震事件的走時(shí)同相軸的走時(shí)初至記錄與相鄰的微地震事件的走時(shí)同相軸的走時(shí)初至記錄的趨勢(shì)相同時(shí)，如果速度模型為P波的速度模型，則單軸事件為P波，如果速度模型為S波的速度模型，則單軸事件為S波。對(duì)微地震事件進(jìn)行掃描疊加的步驟可包括沿所有假設(shè)的微地震事件的走時(shí)設(shè)置預(yù)定窗口，對(duì)微地震事件的能量進(jìn)行掃描疊加，以獲得至少一個(gè)非均一化互相關(guān)系數(shù)；選擇所述至少一個(gè)非均一化互相關(guān)系數(shù)中的最大非均一化互相關(guān)系數(shù)所對(duì)應(yīng)的至少一個(gè)網(wǎng)格點(diǎn)作為微地震事件所對(duì)應(yīng)的網(wǎng)格點(diǎn)。將在接下來的描述中部分闡述本發(fā)明另外的方面和/或優(yōu)點(diǎn)，還有一部分通過描述將是清楚的，或者可以經(jīng)過本發(fā)明的實(shí)施而得知。


通過下面結(jié)合附圖進(jìn)行的詳細(xì)描述，本發(fā)明的上述和其它目的和特點(diǎn)將會(huì)變得更加清楚，其中圖I是示出根據(jù)本發(fā)明示例性實(shí)施例的基于方位角約束掃描疊加的微地震事件定位方法的流程圖；圖2是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的建立精細(xì)速度模型的流程圖；圖3是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的步驟S130的流程圖；圖4是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的步驟S150的流程圖。
具體實(shí)施例方式現(xiàn)在，詳細(xì)描述本發(fā)明的示例性實(shí)施例，其示例在附圖中表示，其中，相同的標(biāo)號(hào)始終表示相同的部件。
圖I是示出根據(jù)本發(fā)明示例性實(shí)施例的基于方位角約束掃描疊加的微地震事件定位方法的流程圖。參照?qǐng)D1，在步驟S110，根據(jù)測(cè)井資料及射孔資料對(duì)井口下包含微地震源的地下空間建立速度模型。優(yōu)選地，這里可對(duì)井口下包含微地震源的地下空間建立精細(xì)速度模型，下面將參照?qǐng)D2來詳細(xì)描述。但應(yīng)了解，采用現(xiàn)有技術(shù)中的普通速度模型同樣能夠?qū)嵤┍景l(fā)明的技術(shù)方案。圖2是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的建立精細(xì)速度模型的流程圖。如圖2所示，在操作S210，根據(jù)測(cè)井資料對(duì)井口下包含微地震源的地下空間劃分層位，建立初始水平層狀速度模型，其中，所述測(cè)井資料可以是諸如測(cè)井聲波時(shí)差、自然電位等。在操作S220，通過使用射孔記錄的走時(shí)差以及已知的射孔點(diǎn)位置，根據(jù)初始水平層狀速度模型，通過遺傳算法反演方法對(duì)初始水平層狀速度模型進(jìn)行修正，以獲得精細(xì)速度模型。這里，遺傳算法反演方法屬于本領(lǐng)域中的現(xiàn)有技術(shù)，為了不模糊本發(fā)明的主題，將不在這里進(jìn)行詳細(xì)描述。 在步驟S120，對(duì)井口下包含微地震源的地下空間建立三維坐標(biāo)系統(tǒng)并進(jìn)行離散網(wǎng)格化，并通過精細(xì)速度模型來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差。這里，僅作為示例，可以以井口或第一級(jí)三分量檢波器作為坐標(biāo)原點(diǎn)、正東為X方向、正北為Y方向、垂直向下為Z方向來建立三維坐標(biāo)系統(tǒng)。優(yōu)選地，這里可使用射線追蹤方法來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差，所述射線追蹤方法可包括如下步驟給定射線的某個(gè)出射角；計(jì)算射線與每個(gè)層位的交點(diǎn)；當(dāng)射線到達(dá)三分量檢波器所在的層位時(shí)，確定能否到達(dá)所述三分量檢波器；如果能夠到達(dá)所述三分量檢波器，則記錄射線路徑；否則，給定射線的一個(gè)新的出射角，并再次進(jìn)行計(jì)算交點(diǎn)與判斷能否到達(dá)所述三分量檢波器的步驟，直到找到可以準(zhǔn)確到達(dá)所述三分量檢波器的出射角，并且記錄射線路徑；通過彎曲法對(duì)射線進(jìn)行微調(diào)，輸出微調(diào)的射線路徑；通過使用精細(xì)速度模型和微調(diào)的射線路徑來獲得假設(shè)的微地震事件的走時(shí)和走時(shí)差。應(yīng)該了解，在具體實(shí)現(xiàn)中，為了通用性，可把層位設(shè)定為變量，通過確定三分量檢波器和震源的所在層位，來控制所述射線追蹤方法中重復(fù)判斷出射角的步驟的起始點(diǎn)和終止點(diǎn)。應(yīng)該理解，使用射線追蹤方法來計(jì)算假設(shè)的微地震事件的走時(shí)和走時(shí)差僅是示例，本領(lǐng)域技術(shù)人員完全可采用其它的方法來計(jì)算假設(shè)的微地震事件的走時(shí)和走時(shí)差。在步驟S130，從實(shí)際采集的微地震資料中識(shí)別微地震事件。這里，可通過使用三分量檢波器來實(shí)際采集微地震資料。在進(jìn)行微地震事件的識(shí)別時(shí)，主要根據(jù)實(shí)際采集的微地震資料中的有效數(shù)據(jù)(SP，微地震事件)與噪聲在能量、偏振性質(zhì)和走時(shí)特點(diǎn)等方面的差異來進(jìn)行，也就是說，微地震資料中可包括微地震事件、噪聲等。下面將參照?qǐng)D3來詳細(xì)描述圖I中的步驟S130。圖3是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的步驟S130的流程圖。如圖3所示，可對(duì)微地震資料進(jìn)行如下操作在操作S310，對(duì)微地震資料設(shè)置預(yù)定門限值，這里，可采用微地震資料的信噪比的3/4作為門限值，但應(yīng)理解這僅是示例，本領(lǐng)域技術(shù)人員可采用其它的值作為門限值，諸如，為適應(yīng)微地震資料的變化，可根據(jù)能量比-偏振度響應(yīng)值按照如下等式I來設(shè)置動(dòng)態(tài)門限值H(t)H(t) =Em(t_ τ ) + a Ev (t_ τ ) (I)
其中，Em(t)是響應(yīng)包絡(luò)面的期望值，Ev(t)是包絡(luò)面的標(biāo)準(zhǔn)差，用來消除信號(hào)波動(dòng)帶來的影響，其中，通過對(duì)微地震資料進(jìn)行希爾伯特(Hilbert)變換來獲得包絡(luò)面；α是用于調(diào)整標(biāo)準(zhǔn)差的權(quán)系數(shù)，τ為延遲樣點(diǎn)數(shù)。在操作S320，選擇微地震資料中的較強(qiáng)的三道數(shù)據(jù)以及合適的窗口長度。這里，窗口的前時(shí)窗(BTA)的長度可較大，窗口的后時(shí)窗(ATA)的長度可較小。在操作S330,沿選擇的三道數(shù)據(jù)在所述窗口內(nèi)進(jìn)行能量大小計(jì)算和偏振分析，以獲得BTA中的數(shù)據(jù)的能量與ATA中的數(shù)據(jù)的能量的比值，并判斷該比值是否達(dá)到預(yù)定門限值。如果達(dá)到預(yù)定門限值，則在操作S340，從所述微地震資料中的較強(qiáng)的三道數(shù)據(jù)中拾取出微地震事件，并對(duì)拾取出的微地震事件截取合適大小。如果未達(dá)到預(yù)定門限值，則在操作S350，進(jìn)行窗口滑動(dòng)，并再 次進(jìn)行操作S330，直到對(duì)所述三道數(shù)據(jù)完成能量大小計(jì)算和偏振分析以及微地震事件的拾取。應(yīng)該理解，圖3所示的識(shí)別微地震事件的步驟僅是示例，本領(lǐng)域技術(shù)人員完全可采用其它的方法，諸如地震波反演能量聚焦方法，來進(jìn)行微地震事件的識(shí)別。 在步驟S140，對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離，并獲得微地震事件的相應(yīng)震相(即，主要(P)波、次要(S)波)的走時(shí)和走時(shí)差。優(yōu)選地，當(dāng)微地震事件是雙軸掃描事件(即，同時(shí)含有P波和S波)時(shí)，通過相應(yīng)震相到時(shí)時(shí)差的差異特征來進(jìn)行微地震震相的識(shí)別，即，如果第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差大于第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差，則第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)S波，第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)P波，反之亦然；當(dāng)微地震有效事件是單軸掃描事件(即，僅含有P波和S波中的一個(gè))時(shí)，首先確定與被掃描的微地震事件相鄰的微地震事件的走時(shí)和位置，當(dāng)被掃描的微地震事件的走時(shí)同相軸的走時(shí)初至記錄與相鄰的微地震事件的走時(shí)同相軸的走時(shí)初至記錄的趨勢(shì)相同時(shí)，如果精細(xì)速度模型為P波的速度模型，則單軸事件為P波，如果精細(xì)速度模型為S波的速度模型，則單軸事件為S波。應(yīng)該理解，上述對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離的方法僅是示例，本領(lǐng)域技術(shù)人員完全可采用其它的方法進(jìn)行微地震震相識(shí)別分離。在步驟S150，通過使用微地震事件的相應(yīng)震相的走時(shí)和走時(shí)差以及每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差對(duì)微地震事件進(jìn)行掃描疊加，以獲得微地震事件所對(duì)應(yīng)的至少一個(gè)網(wǎng)格點(diǎn)。下面將參照?qǐng)D4來詳細(xì)描述圖I中的步驟S150中的微地震事件的掃描疊加操作。圖4是示出圖I的基于方位角約束掃描疊加的微地震事件定位方法的步驟S150的流程圖。如圖4所示，在操作S410，沿所有假設(shè)的微地震事件的走時(shí)設(shè)置預(yù)定窗口，對(duì)微地震事件的能量進(jìn)行掃描疊加，以獲得至少一個(gè)非均一化互相關(guān)系數(shù)?？赏ㄟ^如下等式2來獲得非均一化互相關(guān)系數(shù)E(x,, V1 ,ζ,) = Σ ΣS(χ.,少',，-/，G(rec), j)(2)
red其中，對(duì)一個(gè)微地震事件記錄在對(duì)直達(dá)波部分利用初至進(jìn)行校正之后得到數(shù)據(jù)[xi (j)，j=l,…，M, j=l,…，N]，N是總的三分量檢波器個(gè)數(shù)，M是所選取預(yù)定窗口內(nèi)采樣點(diǎn)的個(gè)數(shù)，E為非歸一化互相關(guān)系數(shù)，G為互相關(guān)預(yù)定窗口范圍內(nèi)的微地震事件，預(yù)定窗口長度可以是I到M之間的一個(gè)值，S為互相關(guān)函數(shù)，i為微地震事件對(duì)應(yīng)的網(wǎng)格號(hào)，j為三分量檢波器號(hào)，可以是I到N之間的一個(gè)值，Xi, Yi, Zi分別為微地震事件的網(wǎng)格坐標(biāo)。在操作S420，選擇所述至少一個(gè)非均一化互相關(guān)系數(shù)中的最大非均一化互相關(guān)系數(shù)(即，最大能量)所對(duì)應(yīng)的至少一個(gè)網(wǎng)格點(diǎn)作為微地震事件所對(duì)應(yīng)的網(wǎng)格點(diǎn)。由于微地震事件可能對(duì)應(yīng)多個(gè)網(wǎng)格點(diǎn)，因此在步驟S160，使用微地震事件所對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)和微地震事件在三維坐標(biāo)系統(tǒng)的三個(gè)分量中的兩個(gè)水平分量來計(jì)
算微地震事件的方位角，并根據(jù)微地震事件的方位角和兩水平分量的關(guān)系tana = +_，以及
Λ
微地震事件對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)，來確定微地震事件的符合所述方位角的發(fā)生的位置，其中，α表示微地震事件的方位角，X、Y分別表示兩個(gè)水平分量。這里，可通過偏振分析方法來計(jì)算微地震事件的方位角，但這僅是示例，本領(lǐng)域技術(shù)人員完全可采用其它的方·法來計(jì)算微地震事件的方位角。應(yīng)該了解，上述步驟S11(TS160僅描述了對(duì)一個(gè)微地震事件確定其位置的處理。然而，如果從實(shí)際采集的微地震資料中識(shí)別出多個(gè)微地震事件，可分別對(duì)每個(gè)微地震事件執(zhí)行上述步驟S11(TS160。通過本發(fā)明的基于方位角約束掃描疊加的微地震事件定位方法，使用有約束的掃描疊加來實(shí)現(xiàn)自動(dòng)定位，較好地解決了定位的多解性問題；此外，本發(fā)明的微地震事件定位方法不需要進(jìn)行人為的初至拾取而對(duì)微地震事件進(jìn)行空間定位，并且對(duì)直井與斜井均可使用所述微地震事件定位方法。雖然已經(jīng)參照特定示例性實(shí)施例示出和描述了本發(fā)明，但是本領(lǐng)域的技術(shù)人員將理解，在不脫離范圍由權(quán)利要求及其等同物限定的本發(fā)明的精神和范圍的情況下可作出形式和細(xì)節(jié)上的各種改變。
權(quán)利要求
1.一種基于方位角約束的微地震事件定位方法，所述方法包括 根據(jù)測(cè)井資料及射孔資料在井口下包含微地震源的地下空間中建立速度模型； 對(duì)井口下包含微地震源的地下空間建立三維坐標(biāo)系統(tǒng)并進(jìn)行離散網(wǎng)格化，并通過速度模型來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差； 從實(shí)際采集的微地震資料中識(shí)別微地震事件； 對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離，以獲得微地震事件的相應(yīng)震相的走時(shí)和走時(shí)差; 通過使用微地震事件的相應(yīng)震相的走時(shí)和走時(shí)差以及每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差對(duì)微地震事件進(jìn)行掃描疊加，以獲得微地震事件所對(duì)應(yīng)的至少一個(gè)網(wǎng)格占. 使用微地震事件所對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)和微地震事件在三維坐標(biāo)系統(tǒng)的三個(gè)分量中的兩個(gè)水平分量來計(jì)算微地震事件的方位角，并根據(jù)微地震事件的方位角以及微地震事件對(duì)應(yīng)的所述至少一個(gè)網(wǎng)格點(diǎn)來確定微地震事件發(fā)生的位置。
2.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，在井口下包含微地震源的地下空間中建立速度模型的步驟包括 根據(jù)測(cè)井資料對(duì)井口下包含微地震源的地下空間劃分層位，建立初始水平層狀速度模型； 通過使用射孔記錄的走時(shí)差以及已知的射孔點(diǎn)位置，根據(jù)初始水平層狀速度模型，通過反演對(duì)初始水平層狀速度模型進(jìn)行修正以獲得精細(xì)速度模型。
3.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，通過使用三分量檢波器來實(shí)際采集微地震資料。
4.如權(quán)利要求3所述的基于方位角約束的微地震事件定位方法，其中，對(duì)井口下包含微地震源的地下空間建立三維坐標(biāo)系統(tǒng)的步驟包括以井口或第一級(jí)三分量檢波器作為坐標(biāo)原點(diǎn)、正東為X方向、正北為Y方向、垂直向下為Z方向來建立三維坐標(biāo)系統(tǒng)。
5.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差的步驟包括使用射線追蹤方法來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差。
6.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，從實(shí)際采集的微地震資料中識(shí)別微地震事件的步驟包括對(duì)微地震資料進(jìn)行如下操作 對(duì)微地震資料設(shè)置預(yù)定門限值； 選擇微地震資料中的較強(qiáng)的三道數(shù)據(jù)以及合適的窗口長度； 沿選擇的三道數(shù)據(jù)在所述窗口內(nèi)進(jìn)行能量大小計(jì)算和偏振分析，以獲得并判斷窗口的前時(shí)窗中的數(shù)據(jù)的能量與窗口的后時(shí)窗中的數(shù)據(jù)的能量的比值是否達(dá)到預(yù)定門限值；如果達(dá)到預(yù)定門限值，則從微地震資料中的較強(qiáng)的三道數(shù)據(jù)中拾取出微地震事件，并對(duì)拾取出的微地震事件截取合適大小，如果未達(dá)到預(yù)定門限值，則進(jìn)行窗口滑動(dòng)并再次進(jìn)行能量大小計(jì)算和偏振分析，直到對(duì)三道數(shù)據(jù)完成能量大小計(jì)算和偏振分析以及微地震事件的拾取。
7.如權(quán)利要求3所述的基于方位角約束的微地震事件定位方法，其中，對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離的步驟包括當(dāng)微地震事件是雙軸掃描事件時(shí)，如果第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差大于第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)的三分量檢波器首尾走時(shí)差，則第一個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)S波，第二個(gè)被掃描的微地震事件的走時(shí)同相軸對(duì)應(yīng)P波，反之亦然。
8.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離的步驟包括 當(dāng)微地震事件是單軸掃描事件時(shí)，確定與被掃描的微地震事件相鄰的微地震事件的走時(shí)和位置，當(dāng)被掃描的微地震事件的走時(shí)同相軸的走時(shí)初至記錄與相鄰的微地震事件的走時(shí)同相軸的走時(shí)初至記錄的趨勢(shì)相同時(shí)，如果速度模型為P波的速度模型，則單軸事件為P波，如果速度模型為S波的速度模型，則單軸事件為S波。
9.如權(quán)利要求I所述的基于方位角約束的微地震事件定位方法，其中，對(duì)微地震事件進(jìn)行掃描疊加的步驟包括 沿所有假設(shè)的微地震事件的走時(shí)設(shè)置預(yù)定窗口，對(duì)微地震事件的能量進(jìn)行掃描疊加，以獲得至少一個(gè)非均一化互相關(guān)系數(shù)； 選擇所述至少一個(gè)非均一化互相關(guān)系數(shù)中的最大非均一化互相關(guān)系數(shù)所對(duì)應(yīng)的至少一個(gè)網(wǎng)格點(diǎn)作為微地震事件所對(duì)應(yīng)的網(wǎng)格點(diǎn)。
全文摘要
提供了一種基于方位角約束的微地震事件定位方法，包括根據(jù)測(cè)井資料及射孔資料對(duì)井口下包含微地震源的地下空間建立速度模型；對(duì)地下空間建立三維坐標(biāo)系統(tǒng)并進(jìn)行離散網(wǎng)格化，并通過速度模型來計(jì)算在每個(gè)網(wǎng)格點(diǎn)中假設(shè)的微地震事件的走時(shí)和走時(shí)差；從實(shí)際采集的微地震資料中識(shí)別微地震事件；對(duì)微地震事件進(jìn)行微地震震相識(shí)別分離，以獲得相應(yīng)震相的走時(shí)和走時(shí)差；通過使用相應(yīng)震相的走時(shí)和走時(shí)差以及假設(shè)的微地震事件的走時(shí)和走時(shí)差對(duì)微地震事件進(jìn)行掃描疊加，以獲得微地震事件所對(duì)應(yīng)的網(wǎng)格點(diǎn)；使用微地震事件所對(duì)應(yīng)的網(wǎng)格點(diǎn)和微地震事件三個(gè)分量中的兩個(gè)水平分量來計(jì)算微地震事件的方位角，并根據(jù)方位角以及網(wǎng)格點(diǎn)來確定微地震事件的位置。
文檔編號(hào)G01V1/28GK102841373SQ20121030134
公開日2012年12月26日 申請(qǐng)日期2012年8月23日 優(yōu)先權(quán)日2012年8月23日
發(fā)明者康亮, 尹陳, 劉鴻, 李亞林, 何光明, 巫芙蓉, 陳愛萍, 巫駿 申請(qǐng)人:中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司