專利名稱：利用微震事件生成的s尾波來識別q因子的制作方法
技術(shù)領域：
本發(fā)明涉及一種用于在油田數(shù)據(jù)采集中使用的方法和系統(tǒng)。具體地， 本發(fā)明提供一種用于利用從在儲液監(jiān)測期間記錄的微震事件獲得的S尾波 來估計Q因子的方法、儀器和系統(tǒng)。
背景技術(shù)：
當?shù)卣鹗录l(fā)生時，地震波從地震事件源傳播。在傳播期間，根據(jù)幾 種類型形式的波的傳播方向和特性，可以識別這些類型的波中的地震能量。 三個主要種類的地震能量一般在地震源已被觸發(fā)后傳播。這些類型的地震 源為橫初波(P波)和兩種類型的剪次波(s波)。剪s波包括縱向剪s波 和殘余s尾波。s尾波為隨著s波之后的能量鏈。這些殘余震動隨著時間 的衰減非常慢。s尾波信號取決于波通過的巖層的特性。
可從在水壓儲液刺激檢測期間記錄的微震事件獲得s尾波。由于微震 事件在地巖層內(nèi)和周圍以不同等級生成，可以通過用s尾波響應校準這樣 的差異，可推斷出關于感興趣的間隔的對比特性。由于s尾波對頻率吸收
敏感，所以可以建立這種關系。這種特性與傳播速度密切相關，其對于微 震事件的準確處理和解釋很重要。
由于s尾波傳播穿過周圍的巖石基質(zhì)，波的能量被周圍的巖石基質(zhì)吸
收，生成關于頻率的幅值和速度的變化，因此改變波形以及減弱高頻。這 些彌散和減弱現(xiàn)象主要是由地巖層內(nèi)的異構(gòu)性引起的。更好的理解這種地 震響應可導致更好的處理這種地巖層內(nèi)的異構(gòu)性。
因此，品質(zhì)因子，或Q因子，成為巖石的重要確定特性，因為Q因子 代表了與地巖層減弱特征的直接關系。Q因子為地表中的地巖層如何減弱 和彌散聲學地震能量的測量值。此外，Q因子為在地巖層中存在碳氫化合 物的指示。通過知道地巖層中的Q因子，可推斷水飽和度以及碳氫化合物 的濃度。Q因子將用于衰減振蕩物理系統(tǒng)的幅值的時間常數(shù)與用于物理系統(tǒng)的振蕩周期進行比較。等效地，Q因子對系統(tǒng)振蕩的頻率與系統(tǒng)散失能 量的比率進行比較。較高的Q因子表示相對于振蕩頻率較低的能量散失比 率。
由于孔隙中的液體特性和基質(zhì)的構(gòu)造影響吸收，Q因子可涉及幾個特 征，諸如孔隙率、滲透率和壓裂度及其他。

發(fā)明內(nèi)容
考慮到上述問題，本發(fā)明的目的是提供一種用于利用在儲液監(jiān)測期間 記錄的地震事件所獲得的S尾波來估計Q因子的方法、儀器和系統(tǒng)。描述 了一種用于建模井場的地巖層和巖石基質(zhì)的方法。記錄微震事件，識別微 震事件中的S尾波窗。S尾波窗包括一組頻率。隨后，識別頻率組中的每個 頻率的Q因子。
微震事件可生成在第一位置，其中，在不同于第一位置的第二位置記 錄微震事件。第一位置可以是井筒，而可從來自于井筒的水力壓裂生成微 震事件。
可用至少一個多分量地震檢波器記錄地震事件。 識別S尾波窗可包括記錄地震事件的幅值和地震事件的頻率。 用于建模井場中的地巖層和巖石基質(zhì)的方法可包括將頻率組分解為頻 域群，以及隨后從多個頻域群中選擇幾個代表性的頻率。識別頻率組的Q 因子可包括識別每個代表性頻率的相應Q因子。分析頻域群的速度的變化 以確定地巖層和巖石結(jié)構(gòu)的變化。Q因子可用于管理儲藏液地或在井筒處 的產(chǎn)量。
在識別了到達的P和S波后，從地震事件軌跡中析取S尾波窗，并將 其用作頻譜分解應用的輸入。考慮到最大可恢復頻率為奈奎斯特頻率，當S 尾波經(jīng)由頻譜分解法被分解成它的單頻時，感興趣的幾個頻率包含在所選 定的數(shù)據(jù)頻段中。利用^"皿^'—與時間關系的曲線的線性趨勢的斜率， 估計出每個單頻的Q因子的值。隨后，對于每個地震檢波器，重復同樣的 操作，并且Q因子與檢波器位置關系的曲線給出隨著深度而感測的衰減 (1/Q)。并且，從Q因子與頻率關系的曲線的斜率，可確定速度、Q因子 與頻率間的關系。利用地巖層內(nèi)的幾個微震事件，可以在一個區(qū)域內(nèi)研究Q因子的變化，以用對比衰減信號來區(qū)分不同地帶。并且，當計算頻譜分解
時，不同滑動窗可用于比較Q因子與傳播時間關系的變化。這些滑動窗在 觀察一個間隔內(nèi)的任何衰減變化中可能是有用的。
通過參考附圖、隨后的說明、以及權(quán)利要求，本發(fā)明的其它目的、特 征和優(yōu)點，對于那些本領域技術(shù)人員將變得顯而易見。


圖1是根據(jù)示例性實施例的油井數(shù)據(jù)采集系統(tǒng)的網(wǎng)絡的表示圖2是示出了根據(jù)示例性實施例的用于處理地震事件的系統(tǒng)的框圖3描述了根據(jù)示例性實施例的數(shù)據(jù)處理系統(tǒng)的圖表；
圖4描述了根據(jù)示例性實施例所示的典型的幅值與頻率關系的曲線圖5根據(jù)示例性實施例，描述了用于2—4Hz頻率的、己經(jīng)經(jīng)過頻譜分
解的尾波窗的典型的帶通濾波后的位移幅值；
圖6根據(jù)示例性實施例，描述了用于12—24Hz頻率的、已經(jīng)經(jīng)過頻譜
分解的尾波窗的典型的帶通濾波后的位移幅值；
圖7是根據(jù)示例性實施例的，用于2—4Hz頻率的、選定的尾波窗的典
型的幅值乘以間隔時間的均方根，及其最佳乘方擬合；
圖8是根據(jù)示例性實施例的，用于12—24Hz頻率的、選定的尾波窗的
典型的幅值乘以間隔時間的均方根，及其最佳乘方擬合；
圖9是根據(jù)示例性實施例的描述示出的過程步驟的高級流程圖10是根據(jù)示例性實施例的用于識別示出的至少一個地震事件的流程
圖11是根據(jù)示例性實施例的用于形成示出的地震事件的幅值與頻率關 系分析的流程圖12是根據(jù)示例性實施例的用于用頻率分析速度的變化以為油田確定 Q因子的變化的流程在下面的本發(fā)明的優(yōu)選實施例及其它實施例的詳細的描述中，對附圖 進行參考。要理解的是，本領域技術(shù)人員將容易理解其它不脫離本發(fā)明范 圍的實施例和改變。
具體實施例方式
在典型的鉆探操作中，幾個井筒戰(zhàn)略上圍繞油田分布以便于最大化油 田內(nèi)含有的油或氣的產(chǎn)量。在井筒鉆入含有油、天然氣和水的地巖層巖石 后，油和氣的產(chǎn)量最大化是一個首要的目標。可通過改善或最大化油井的 液體流量來增加油和氣的產(chǎn)量。通過使許多先前存在的裂縫和儲液巖石中 的流動通道與較大的裂縫連接，較大的裂縫充當管道，更加快且有效地從 地巖層巖石引導液體到井筒。這個較大的、人造的裂縫開始于油井并向外 延伸入儲液巖石多達數(shù)百英尺。
通過記錄從油井鉆空中的水力壓裂所生成的微震事件，可以識別微震
事件的S尾波窗。隨后識別出S尾波窗內(nèi)的一組頻率的Q因子。通過知道 Q因子，可以推斷出地巖層中的水飽和度和碳氫化合物濃度。因此，可以 對井場和/或井場所位于的儲油層的鉆探和碳氫化合物的產(chǎn)量，作出更精明 的決定。
現(xiàn)在參考圖1，描述了可以實現(xiàn)本發(fā)明的示例性實施例的網(wǎng)絡數(shù)據(jù)處理 系統(tǒng)的表示圖。在這個實例中，網(wǎng)絡數(shù)據(jù)處理系統(tǒng)100為可以在其中實現(xiàn) 本發(fā)明的不同實施例的計算設備網(wǎng)絡。網(wǎng)絡數(shù)據(jù)處理系統(tǒng)100包括網(wǎng)絡102， 其為用于在網(wǎng)絡數(shù)據(jù)處理系統(tǒng)100內(nèi)互相通信的各種設備和計算機之間提 供通信鏈路的介質(zhì)。網(wǎng)絡102可包括連接，諸如導線、無線通信鏈路或光 纜。通過將數(shù)據(jù)存儲于諸如硬盤驅(qū)動、DVD或閃存器的存儲設備上，數(shù)據(jù) 甚至可手動傳遞。
在該所描述的實例中，井場104、 106、 108和110具有用于生成關于 位于這些井場的油井的數(shù)據(jù)的計算機或其它計算設備。在這些實例中，井 場104、 106、 108和IIO位于地理區(qū)域112。正如在這些實例中所描述的， 此地震檢波器為單個儲油層。當然，取決于具體的實施，這些井場可跨不 同的地理區(qū)域和/或跨多個儲油層。井場104和106具有到網(wǎng)絡102的有線 通信鏈路114和116。井場108和110具有到網(wǎng)絡102的無線通信鏈路118 和120。
分析中心122為一位置，在該位置設置有數(shù)據(jù)處理系統(tǒng)以處理從井場 104、 106、 108和110收集的地震數(shù)據(jù)和其它數(shù)據(jù)。取決于特定的實施，可 呈現(xiàn)多個分析中心。例如，取決于具體實施，這些分析中心可在辦公室或地理位置112現(xiàn)場。在這些示例性實施例中，分析中心122利用用于本發(fā) 明的不同實施例的過程，分析來自于井場104、 106、 108和110的數(shù)據(jù)。
在該所描述的實例中，網(wǎng)絡數(shù)據(jù)處理系統(tǒng)100為具有網(wǎng)絡102和網(wǎng)關 的互聯(lián)網(wǎng)，網(wǎng)絡102代表世界范圍網(wǎng)絡集合，網(wǎng)關使用傳輸控制協(xié)議/互聯(lián) 網(wǎng)協(xié)議(TCP/IP)的協(xié)議組來互相通信。位于互聯(lián)網(wǎng)中心的是主節(jié)點或主 機之間的高速數(shù)據(jù)通信線的主干，包括數(shù)千個用于發(fā)送數(shù)據(jù)和信息的商業(yè)、 政府、教育及其它計算機系統(tǒng)。當然，網(wǎng)絡數(shù)據(jù)處理系統(tǒng)100也可以實施 為許多不同類型的網(wǎng)絡，例如，內(nèi)聯(lián)網(wǎng)、局域網(wǎng)(LAN)或廣域網(wǎng)(WAN)。 圖1是用來作為實例，而不是作為用于不同實施例的建筑上的限制。該信 息可用于井場的管理或潛在的鉆探基地的確定。例如，信息可用于指導井 場的鉆探操作。
現(xiàn)在參考圖2，根據(jù)示例性實施例示出了闡明用于處理微震事件的系統(tǒng) 的框圖。通過將壓裂液體灌注到諸如圖1的井場104、 106、 108和110的 井場的井筒212中，形成水力壓裂210。以一個速率將壓裂液體灌注到井筒 212，該速率足以增加向下鉆進壓力，使其超過地巖層巖石214的斷裂梯度 的值。隨后增加的壓力使得地巖層巖石214破裂，這允許壓裂液體進入并 使裂縫進一步延伸到地巖層。地巖層巖石214的破裂及由此的任何擴大或 張馳，生成微震事件216。
微震事件216是由釋放所儲存的應變能量的巖石的突然破裂引起的振 動或一系列振動。微震事件216生成地震波218，其包括有壓縮性的初波(P 波)、剪次波(S波)和殘余S尾波。這些波遠離微震事件216傳播并被接 收器220檢測到。接收器220可以是VSI通用地震成像儀，其可從斯倫貝 謝有限公司(Schlumberger Ltd)獲得。在接收器120內(nèi)的地震檢波器可以 是地震檢波器加速度計(GAC-D)，也可從斯倫貝謝有限公司(Schlumberger Ltd)獲得。接收器220是能夠檢測微震事件216的靈敏的地震儀器。
微震水力壓裂監(jiān)測為用于檢測由在巖石基質(zhì)的最初應力位置的變化引 起的微聲能量的鉆探地震技術(shù)。通過在生成儲油層內(nèi)或附近放置接收器 220，通常在第二鄰近的井筒，儲油層中的微震活動的連續(xù)的測量表現(xiàn)了正 在那個儲油層內(nèi)變化的物理特性?？梢栽谠搩τ蛯觾?nèi)變化的物理特性包括， 但不限于，應力破裂、流量、流動性和熱特性。接收器220可以是放置于監(jiān)測油井222中的一個或多個多分量地震檢 波器228-236，用于確定微震事件216的位置。在這些實例中，地震檢波器 228-236為高靈敏度地面移動傳感器，其通過將地面運動轉(zhuǎn)換為電壓來測量 地面內(nèi)的振動，可在記錄站記錄該電壓。此測量電壓從基線的偏離被稱為 地震響應，并且其用于分析地球結(jié)構(gòu)。
地震檢波器228-236 —般由懸掛在彈簧上的質(zhì)點組成。導線纏繞在質(zhì)點 周圍。相對于周圍地巖層靜止的磁鐵圍繞在被導線纏繞的質(zhì)點周圍。當?shù)?面移動時，諸如當微震事件發(fā)生時，質(zhì)點趨向于保持不動，而磁鐵與周圍 的地巖層一起圍繞質(zhì)點移動。移動磁鐵的磁場在導線中產(chǎn)生電壓。該電壓 通過簡單的電壓計放大和記錄。由電壓計記錄的電壓與地面移動處的震級 成比例。
可通過測量壓力或初級(P-)波與剪力或次級(S-)波間的波至時間的 差值來計算接收器220到微震事件216的距離。用于研究P波的粒子運動 的Hodogmm分析可同樣用于確定事件的方位角。通過利用在監(jiān)測油井222 所觀察的接收器間的P-和S-波的波至延遲來限定微震事件216的深度。該 定位技術(shù)需要用于計算事件位置的準確的速度模型以、用于記錄微震事件 216的低噪聲環(huán)境高靈敏度地震檢波器、和接收器220的確切位置及方向的 知識。
接收器220被配置在井筒212的2000英尺內(nèi)的監(jiān)測油井222中，在大 約與地巖層巖石214同樣的深度處。監(jiān)測油井222 —般為油田內(nèi)的生產(chǎn)井 筒，其中接收器220用于測量諸如微震事件216的微震事件。監(jiān)測油井222 內(nèi)的接收器220的最佳放置和幾何形狀主要取決于周圍的巖石基質(zhì)結(jié)構(gòu)。 然而，由于其它井場的變化，接收器220和井筒212之間的理想空間配置 很難達到。
接收器220發(fā)送從地震事件216收集到的數(shù)據(jù)給數(shù)據(jù)處理系統(tǒng)224以 用于分析。數(shù)據(jù)處理系統(tǒng)224可以是圖1的分析中心122。運行在數(shù)據(jù)處理 系統(tǒng)224上的軟件組件226從微震事件116中分解出S尾波。軟件組件226 隨后利用分解出的S尾波確定地巖層巖石214的Q因子。該信息可用于管 理井場或確定潛在的鉆探基地。例如，該信息可用于指導井場的鉆探操作。
現(xiàn)在轉(zhuǎn)到圖3，描述了根據(jù)示例性實施例的數(shù)據(jù)處理系統(tǒng)的圖表。數(shù)據(jù)
ii處理系統(tǒng)300可以是圖1的數(shù)據(jù)處理系統(tǒng)124。在這個示例性實例中，數(shù)據(jù) 處理系統(tǒng)300包括通信結(jié)構(gòu)302，其提供在處理器單元304、存儲器306、 永久性存儲器308、通信單元310、輸入/輸入(I/O)單元312和顯示器314
之間的通信。
處理器單元304用于執(zhí)行用于軟件的指令，其可加載到存儲器306。取 決于具體的實施，處理器單元304可以是一個或'多個處理器的組，或可以 是多個處理器中心。此外，可利用一個或多個不同的處理系統(tǒng)實施處理器 單元304，其中，主處理器與次級處理器一起存在于單個芯片上。在另一個 示例性實例中，處理器單元304可以是包含同樣類型的多個處理器的對稱 多處理器系統(tǒng)。
在這些實例中，存儲器306，例如可以是隨機存取存儲器或任何其它適 當?shù)囊资曰蚍且资源鎯υO備。取決于具體實施，永久性存儲器308可 使用多種形式。例如，永久性存儲器308可包含一個或多個組件或設備。 例如，永久性存儲器308可以是硬盤驅(qū)動、閃存存儲器、可重寫光盤、可 重寫磁帶或上述的一些組合。永久性存儲器308所用的介質(zhì)也可以是移動 的。例如，移動硬盤驅(qū)動可用于永久性存儲器308。
在這些實例中，通信單元310提供與其它數(shù)據(jù)處理系統(tǒng)或設備的通信。 在這些實例中，通信單元310是網(wǎng)絡接口卡。通信單元310可通過使用物 理和無線通信鏈路其中之一或兩者來提供通信。
輸入/輸出單元312允許其它可連接到數(shù)據(jù)處理系統(tǒng)300的設備輸入及 輸出數(shù)據(jù)。例如，輸入/輸出單元312可通過鍵盤和鼠標為用戶輸入提供連 接。此外，輸入/輸出單元312可將輸出發(fā)送到打印機。顯示器314向使用 者提供用于顯示信息的機構(gòu)。
用于操作系統(tǒng)的指令和應用或程序設置在永久性存儲器308上。這些 指令可加載到存儲器306，以用于由處理單元304執(zhí)行。利用執(zhí)行設置于諸 如存儲器306的存儲器中的指令的計算機，處理器304可完成不同的實施 例的程序。這些可由處理器單元304中的處理器讀取和執(zhí)行的指令可被稱 作程序代碼、計算機可用程序代碼或計算機可讀程序代碼。不同實施例中 的程序代碼可包含在不同的物理的或?qū)嶓w的計算機可讀介質(zhì)上，諸如存儲 器306或永久性存儲器308。
12將程序代碼316以函數(shù)的形式設置在計算機可讀介質(zhì)318上且加載到 或傳遞到數(shù)據(jù)處理系統(tǒng)300，以用于通過處理單元304執(zhí)行。在這些實例中， 程序代碼316和計算機可讀介質(zhì)318形成計算機程序產(chǎn)品320。在一個實例 中，計算機可讀介質(zhì)318可以是實體形式，例如，被插入或放入到作為永 久性存儲器308 —部分的驅(qū)動或其它設備的光或磁盤，以用于轉(zhuǎn)移到存儲 設備，諸如為永久性存儲器308 —部分的硬盤驅(qū)動。在實體形式中，計算 機可讀介質(zhì)318也可采取永久性存儲器的形式，諸如連接到數(shù)據(jù)處理系統(tǒng) 300的硬盤驅(qū)動或閃存存儲器等。計算機可讀介質(zhì)318的實體形式也被稱作 計算機可錄存儲器介質(zhì)。
可選擇地，通過到通信單元310的通信鏈路和/或通過到輸入/輸出單元 312的連接，程序代碼316可從計算機可讀介質(zhì)318傳遞到數(shù)據(jù)處理系統(tǒng) 300。通信鏈路和/或連接可以是示例性實例中的物理的或無線的形式。計算 機可讀介質(zhì)也可采取非實體介質(zhì)的形式，諸如含有程序代碼的通信鏈路或 無線傳輸。
所闡述的用于數(shù)據(jù)處理系統(tǒng)的不同元件不意味著限制不同的實施例的 構(gòu)造實施方式?？稍诎ǔ嘶蛱鎿Q那些用于闡述數(shù)據(jù)處理系統(tǒng)300的元 件的數(shù)據(jù)處理系統(tǒng)中實施不同的示例性實施例。圖3中示出的其它元件可 不同于示出的示例性實施例。
例如，總線系統(tǒng)可用于實現(xiàn)通信結(jié)構(gòu)302，并且總線系統(tǒng)可由諸如系統(tǒng) 總線或輸入/輸出總線的一個或多個總線組成。當然，可利用用于提供連接 到總線系統(tǒng)的不同元件或設備間的數(shù)據(jù)傳遞的任何適當類型的構(gòu)造來實現(xiàn) 總線系統(tǒng)。另外，通信單元可包括一個或多個用于傳輸和接收數(shù)據(jù)的設備， 諸如調(diào)制解調(diào)器或網(wǎng)絡適配器。此外，存儲器可以是存儲器306或高速緩 存器，其處于可出現(xiàn)在通信結(jié)構(gòu)302中的接口或內(nèi)存控制中心。
現(xiàn)在參考圖4，根據(jù)示例性實施例示出了典型的幅值與頻率關系曲線。 圖4為諸如圖2中的微震事件216的微震事件的記錄，該微震事件已經(jīng)被 諸如圖2中的接收器220的接收器記錄。
微震事件發(fā)生在初始時刻410。在時刻p412，較快移動的P橫波到達 諸如圖1中的接收器120的微震接收器。在時刻s414，剪S波到達微震接 收器。時刻c416指示S尾波窗418的開始。時刻c416發(fā)生在等于時刻s414
13和初始時刻410之間的差的兩倍的時刻。S尾波窗418的末端通常在一個點 上，在該點上S尾波信號衰減了，以至于不能從來自于接收器的固有噪聲 中區(qū)分該點。S尾波窗418包含與微震事件相關的S尾波數(shù)據(jù)。
現(xiàn)在參考圖5，根據(jù)示例性實施例描述了濾波后的S尾波數(shù)據(jù)的圖表。 圖5為S尾波窗(諸如圖4的S尾波窗418)的帶通濾波后的位移幅值的典 型的幅值與頻率的關系分析，該S尾波窗已經(jīng)經(jīng)過頻譜分解。圖5為在2 一4Hz頻率處的S尾波窗的帶通濾波后的位移幅值。
諸如圖4的S尾波窗418的尾波窗經(jīng)過2—4Hz的窄帶濾波器，以將圖 4的S尾波窗418分解為具有窄頻段的濾波后的S尾波窗510。 一般，越高 的地震頻率反映越薄的地層，而越低的地震頻率反映越厚的地層。頻譜分 解允許利用地震頻帶的離散分量。通過利用離散的頻率分量，可以更準確 地觀察和辨別出地巖層中的不同地層。
窄帶濾波器可以是任何能將S尾波信號分解成離散的波長段的濾波器。 窄帶濾波器可以是，例如但不限于，多層帶通濾波器、前端濾波器、介質(zhì) 濾波器、高介電常數(shù)陶瓷濾波器、微帶濾波器、光學帶通濾波器、交指型 帶通濾波器、有限沖激響應(FIR)濾波器或其它本領域已知的頻帶濾波器。
現(xiàn)在參考圖6，根據(jù)示例性實施例描述了濾波后的S尾波數(shù)據(jù)。圖6 為已經(jīng)經(jīng)過頻譜分解的S尾波窗(諸如圖4的S尾波窗418)的帶通濾波后 的位移幅值的典型的幅值與頻率的關系分析。圖6為在12—24Hz頻率處的 S尾波窗的帶通濾波后的位移幅值。
尾波窗(諸如圖4的S尾波窗418)經(jīng)過12—24Hz的窄帶濾波器以將 圖4的S尾波窗418分解為具有窄頻段的濾波后的S尾波窗610。 一般，越 高的地震頻率反映越薄的地層，而越低的地震頻率反映越厚的地層。頻譜 分解允許利用地震頻帶的離散分量。通過利用離散的頻率分量，可以更準 確地觀察和識別出地巖層中的不同地層。
窄帶濾波器可以是任何能將S尾波信號分解成離散的波長段的濾波器。 窄帶濾波器可以是例如但不限于，多層帶通濾波器、前端濾波器、介質(zhì)濾 波器、高介電常數(shù)陶瓷濾波器、微帶濾波器、光學帶通濾波器、交指型帶 通濾波器、有限沖激響應(FIR)濾波器或其它本領域已知的頻帶濾波器。
基于從圖5和圖6的濾波后的S尾波窗510和610中的頻譜分解獲得
14的結(jié)果，從瞬時幅值計算均方根幅值?？衫孟旅娴墓綇乃矔r幅值計算 均方根(RMS)幅值
其中，P為瞬時幅值，并且N為時域波形中的幅值測量值。
接下來確定S尾波的幅值與它們被記錄器記錄的時刻間的關系?？梢?通過繪制RMS幅值與時間的乘積取自然對數(shù)與時間的關系的曲線，將該關 系線性表征。就是說，所述關系可被線性表征為
Z〃(^MS ""乃'麼)與時間。
現(xiàn)在參考圖7，其根據(jù)示例性實施例，描述了闡明S尾波的幅值與記錄
S尾波的時刻間的關系的圖表。圖7描述了圖5的濾波后的S尾波窗510
的典型的幅值乘以間隔時間的均方根，及其最佳乘方擬合。就是說，圖6 為^"皿^'^)與時間的圖示。
圖7用于2-4Hz頻率。使用中心頻率3.0Hz的均方根幅值，中心頻率 3.0Hz是圖5中識別的帶通濾波后的2—4Hz的位移頻段的中心頻率。
基于S尾波的幅值與它們被記錄的時刻間的關系，可以確定圖4中識 別的某些濾波后的頻率的Q因子。可利用下面的關系計算Q因子
其中，f為來自于圖5中識別的2—4Hz的帶通濾波后的位移頻段的中心頻
率3.0Hz;且b為，例如在圖7中獲得的"，^'^)與時間的關系的曲線
的線性趨勢的斜率。建立了用于2—4Hz頻率Q因子407。
現(xiàn)在參考圖8,根據(jù)示例性實施例，描述了闡明S尾波的幅值與S尾波
被記錄的時刻間的關系的圖表。圖8描述了圖6的濾波后的S尾波窗610
的典型的幅值乘以間隔時間的均方根，及其最佳乘方擬合。就是說，圖6 為丄A^4羅4!'附e)與時間的圖示。
圖8用于2-4Hz頻率。所使用中心頻率18.0Hz的均方根幅值，中心頻
率18.0Hz是圖6中識別的帶通濾波后的12—24Hz的位移頻段的中心頻率。
基于S尾波的幅值與它們被記錄的時刻間的關系，可以確定圖4中識
別的某些濾波后的頻率的Q因子。可利用下面的關系計算Q因子 2 = n*//z>
15其中，f為來自于圖6中識別的12—24Hz的帶通濾波后的位移頻段的中心 頻率18.0Hz;且b為，例如在圖7中獲得的^"麗^z'歴)與時間的關系的 曲線的線性趨勢的斜率。建立了用于12—24Hz頻率Q因子2618。
隨后在每個選定的帶通濾波后的頻率處，確定每個接收器的剩余地震 檢波器的Q因子。通過分析在每個接收器(諸如圖2的接收器220)位置 處的Q因子，可確定在整個地巖層的各個深度處的Q因子?？梢栽诮邮掌?的位置上畫出為每個地震檢波器獲得的各種Q因子的值，以分析Q根據(jù)深 度的變化。
對于每個選定的頻率，在記錄器的每個地震檢波器處確定Q因子。可 以在地震檢波器的位置上畫出為每個地震檢波器獲得的各種Q因子的值， 以分析Q根據(jù)深度的變化。
現(xiàn)在參考圖9，根據(jù)示例性實施例示出了過程步驟的高級流程圖。過程 900為運行在諸如圖2的軟件組件226的軟件組件上的軟件過程，該軟件組 件在諸如圖2的數(shù)據(jù)處理系統(tǒng)224的數(shù)據(jù)處理系統(tǒng)上執(zhí)行。
過程800通過識別一組微震事件開始(步驟910)。 一組微震事件是至 少一個諸如圖2的微震事件216的微震事件。在諸如圖2的監(jiān)測油井222 的監(jiān)測油井處，由諸如圖2的接收器220的接收器記錄一組微震事件。識 別出的一組微震事件優(yōu)選地具良好的信噪比，帶有清楚限定的P波和S波 波至以及限定的S尾波窗。
響應于識別一組微震事件，對一組微震事件執(zhí)行幅值與頻率關系的分 析(步驟920)。幅值與頻率關系的分析可以是，如圖6中所描述的，由諸 如圖2的接收器220的接收器記錄的， 一組微震事件的幅值和頻率的曲線 圖。
利用向前傅里葉變換，可從由響應于微震事件的接收器所生成的電壓 中確定幅值和頻率。執(zhí)行幅值和頻率的分析，以識別S尾波窗的開始。S 尾波窗的開端通常位于如圖4的時刻c416所指示的2ts時刻。尾波窗的末 端較模糊尾波窗的末端通常取S尾波信號已經(jīng)衰減到不再能夠從來自于 接收器的固有噪聲中識別出的點。
位于特定頻率處的尾波幅值隨著間隔時間的減小僅僅是由于能量從巖 石基質(zhì)的減弱和彌散導致的，而與地震源、傳播路徑和基地放大無關。地
16震波的衰減為固有的和分散的衰減的總和，其中在第一種情形中，能量轉(zhuǎn) 換成熱量，而在第二種情形中，其通過在整個儲油層中隨機的折射、反射 和散射來重新分配能量。
響應于對一組微震事件執(zhí)行幅值與頻率關系的分析，過程900分析速 度關于頻率的變化及最終儲油層Q因子中的變化(步驟930)。儲油層內(nèi)Q 因子中的變化可用于更好的表征地巖層內(nèi)的各種不同的地層。可從這個儲 油層的Q因子中的變化中建模潛在油田、產(chǎn)量估計及井油田巖層的一般特 性。
響應于分析速度關于頻率的變化和儲油層的Q因子中的變化，過程900 識別地巖層的總Q因子(步驟940)。隨后過程900利用Q值管理井場(步 驟950),之后過程結(jié)束。通過知道Q因子，可推斷地巖層的水飽和度和碳 氫化合物濃度。因此能對井場處的鉆探和碳氫化合物的產(chǎn)量做出更精明的 決定。例如，信該息可用于指導井場的鉆探操作。
現(xiàn)在參考圖10，根據(jù)示例性實施例示出了用于識別至少一組微震事件 的流程圖。過程1000為圖9的步驟910的處理步驟的詳圖。過程1000發(fā) 生在諸如圖2的接收器220的接收器內(nèi)。
過程1000通過監(jiān)測一組微震事件開始(步驟IOIO)。 一組微震事件為 由釋放儲存的壓力能量的巖石突然破裂而引起的振動或一系列振動。 一組 微震事件生成地震波，其包括可壓縮的初波(P波)、剪次波(S波)和殘 余S尾波。上述波從一組微震事件向外傳播并被至少一個諸如圖2的接收 器220的接收器檢測。
在這些實例中，接收器被配置在受到水壓或其它刺激的油井的2000英 尺內(nèi)的監(jiān)測井場，位于與處理間隔大約相同的深度。監(jiān)測油井內(nèi)的接收器 的最佳放置和幾何形狀有助于減小信號中的噪聲和衰減，當在接收器處檢 測到P波、S波和S尾波時，將它們的到達時刻分隔開。也就是，在諸如 圖4中所描述的典型的幅值與頻率關系的曲線中，最佳放置允許P波、S 波和S尾波互相區(qū)分。監(jiān)測油井內(nèi).的接收器的最佳放置和幾何形狀嚴重依 賴于周圍的速度結(jié)構(gòu)，因此準確的大地模型有助最佳化監(jiān)測布局。
過程IOOO記錄由接收器生成的電壓信號，并將所記錄的數(shù)據(jù)發(fā)送到數(shù) 據(jù)處理系統(tǒng)以用于分析(步驟1020)，之后過程結(jié)束。數(shù)據(jù)處理系統(tǒng)可以是
17圖2的數(shù)據(jù)處理系統(tǒng)224。所識別的地震事件組最好具有好的信噪比，帶有 清楚限定的P波和S波波至以及限定的S尾波窗。
現(xiàn)在參考圖11，根據(jù)示例性實施例示出了用于形成一組微震事件的幅 值與頻率關系的分析的流程圖。過程1100為運行在諸如圖2的軟件組件226 的軟件組件上的軟件過程，在諸如圖2的數(shù)據(jù)處理系統(tǒng)224的數(shù)據(jù)處理系 統(tǒng)上執(zhí)行該軟件組件。過程1100為圖9的過程步驟920。
過程1100通過接收來自于諸如圖2的接收器220的接收器的電壓信號 開始(步驟1110)。接收器由一個或多個地震檢波器組成。地震檢波器一般 由懸掛于彈簧上的質(zhì)點組成。導線纏繞在質(zhì)點上。相對于周圍的地巖層靜 止的磁鐵，圍繞在纏繞有導線的質(zhì)點上。當大地移動時，諸如當微震事件 發(fā)生時，質(zhì)點趨向于保持不動，而磁鐵在質(zhì)點周圍關于周圍的地巖層移動。 移動磁鐵的磁場在導線中生成電壓。該電壓被簡單的電壓計放大和記錄。 由電壓計記錄的電壓與大地正在移動處的震級成比例。
隨后，過程1100使用向前傅里葉變換來將從地震檢波器生成的電壓分 解成其分量的幅值和頻率(步驟1120)，隨后過程結(jié)束。幅值與頻率關系的 分析可以是諸如圖6中所描述的地震事件的幅值和頻率的曲線圖。
從生成的曲線圖中，容易看出微震事件內(nèi)的相關事件。微震事件在初 始時刻發(fā)生。在時刻p，諸如圖4中的時刻p412，較快移動的P橫波到達 接收器，諸如圖2的接收器1220。在時刻s，諸如圖4的時刻s414，剪S 波到達微震接收器。S尾波窗，諸如圖43的S尾波窗418，在時刻c，諸如 圖4中的時刻c416，開始。時刻c發(fā)生在等于時刻s和微震事件發(fā)生的初 始時刻間的差值的兩倍的時刻。S尾波窗的末端較模糊S尾波窗的末端通 常取S尾波信號已經(jīng)衰減到不再可從來自于接收器的固有噪聲中識別出其 的點。
位于特定頻率處的尾波幅值隨著間隔時間的減小不但是由能量從巖石 基質(zhì)的減弱和彌散導致的，而且與地震源、傳播路徑和基地放大無關。地 震波的衰減是固有的衰減和分散的衰減的總和，其中在固有的衰減中將能 量轉(zhuǎn)換成熱量，在分散的衰減中在整個儲油層中通過隨機的折射、反射和 散射來重新分配能量。
現(xiàn)在參考圖12,根據(jù)示例性實施例示出了用于估計Q因子的流程圖。
18在這些實施例中，過程1100用于分析速度根據(jù)頻率的變化，以確定油田Q 因子中的變化。過程1200是在諸如圖2的軟件組件226的軟件組件上執(zhí)行 的軟件過程，并且過程1200為圖9的步驟930的更詳細的說明。
過程1200通過接收為微震事件所執(zhí)行的幅值與頻率關系的分析開始 (步驟1210)。如圖6中所描述的，幅值與頻率關系的分析可以是接收器(被 諸如圖2的接收器220)記錄的微震事件的幅值與頻率關系的曲線。幅值與 頻率關系的分析可以圖9的步驟920所描述的分析。
S尾波窗的開端通常取如圖4的在時刻c416指示的時刻2ts。 S尾波窗 的末端較模糊S尾波窗的末端通常取S尾波信號已經(jīng)衰減到S尾波信號 不再可從來自于接收器的固有噪聲中識別出其的點。
選擇包含S尾波窗的時間窗(步驟1215)。 S尾波窗的開端通常取如圖 4的在時刻c416指示的時刻2ts。也就是，S尾波窗的開端通常取時刻2ts， 或
其中，②da^為s尾波窗的開端，S',為S波窗的開端，而A^為微震事 件初始時刻。
S尾波窗的末端較模糊S尾波窗的末端通常取S尾波信號已經(jīng)衰減到 不再可從接收器中的固有噪聲中識別出其的點。
基于S尾波窗的長度，將滑動窗應用于S尾波窗(步驟1220)?；瑒哟?是用于尾波窗處置的統(tǒng)計平滑技術(shù)，用于使S尾波時間序列的長期走勢更 清楚。盡管實施例描述為使用滑動窗，也可使用其它統(tǒng)計處理，諸如移動 平均。同樣地，可完全不采用數(shù)據(jù)的統(tǒng)計平滑。
基于S尾波窗的大小選擇初始滑動窗。隨后選擇步長，該步長用于在 接下來的運行中改變窗大小。例如，如果識別出S尾波窗具有150毫秒的 長度，應用于統(tǒng)計平滑的初始窗可以是30毫秒(ms)。如果確定有4個運 行，識別出隨后的滑動窗的大小可為60ms、 90ms和120ms。
隨后對S尾波窗的每個識別出的滑動窗執(zhí)行頻譜分解(步驟1225)。可 用窄帶濾波器執(zhí)行頻譜分解。來自于地巖層的反射具有頻域中的特征表達 式，其指示地巖層的厚度。越高的分解頻率一般提供更準確的越薄的地巖 層的映像，而越低的分解頻率提供更準確的越厚的地巖層的映像。頻譜分
19解將時間或深度地震數(shù)據(jù)轉(zhuǎn)換到頻域。可利用諸如離散傅里葉變換或最大 熵方法等的適當?shù)某绦蚧蛩惴▽崿F(xiàn)該頻譜分解。頻譜分解允許工程師或其 它使用者獨立地分析用于揭示各種地巖層的先前模糊的特征的每個頻率。
窄帶濾波器可以是任何能將s尾波信號分解成為離散波長段的濾波器。
窄帶濾波器可以是例如但不限于，多層帶通濾波器、前端濾波器、介質(zhì)濾 波器、高介電常數(shù)陶瓷濾波器、微帶濾波器、光學帶通濾波器、交指型帶
通濾波器、有限沖激響應(FIR)濾波器或其它本領域已知的頻帶濾波器。 當S尾波被分解成它的單頻時，過程1200隨后選擇幾個包含在記錄的 數(shù)據(jù)頻帶內(nèi)的感興趣的代表性的頻率(步驟1230)。可以通過協(xié)作數(shù)據(jù)分析 過程，通過獨立的分析確定這些代表性的頻率。相反地，可以基于典型觀 察的代表性的頻率，預先確定該代表性頻率。代表性頻率還可用于定位特 定的使用者輸入。
基于來自于幅值與頻率關系的分析的結(jié)果，過程1200從瞬時幅值計算 均方根幅值(步驟1235)。可利用下面的公式從瞬時幅值計算均方根(RMS) 幅值
其中，P為瞬時幅值，并且N為時域波形中的幅值測量值。
此外，過程1200確定S尾波的幅值與它們被記錄器記錄的時刻之間的
關系(步驟1240)。通過繪制RMS幅值與時間的乘積取自然對數(shù)與時間的
關系的曲線，可以線性表征該關系。就是說，過程1200繪制曲線 丄W(^咖s *乃'附e)與時間。
圖7和8中示出了典型的^"廳"^e)與時間的曲線。這里所使用的 均方根幅值是從分別在圖5和圖6中識別出的位移頻率范圍2—4Hz和12 一24Hz帶通濾波出的中心頻率3.0和18.0的均方根幅值。
基于在步驟1240中建立的關系，為幾個代表性的頻率計算Q因子(步 驟1245)。過程1200可計算Q因子，其利用下面的關系
其中，b為^"麗^^e)與時間的關系的曲線的線性趨勢的斜率；且f為 頻率。
如圖7和8中所表明的，建立了用于2—4Hz頻率的Q因子407，同時
20建立了用于12—24Hz頻率的Q因子2618。諸如圖2的數(shù)據(jù)處理系統(tǒng)224 的數(shù)據(jù)處理系統(tǒng)確定結(jié)果曲線的線性趨勢的斜率。對于每個地震檢波器的 每個選定的頻率，確定^^"皿#77^)與時間的關系的曲線。
一旦針對幾個代表性頻率計算出Q因子，對于每個地震檢波器獲得的 各種Q因子的值可根據(jù)地震檢波器的位置繪制曲線，以分析Q因子隨著深 度的變化(步驟1250)。為每個選定的頻率獲得單獨的Q因子。
過程1200隨后根據(jù)選定的頻率，為每個選定的頻率(步驟1255)繪制所 獲得的Q因子的曲線，以分析它們的習性并獲得總的Q因子。根據(jù)圖表的 斜率，考慮到下面的關系式來確定速度根據(jù)頻率的變化-
其中，b為，例如圖7和8中獲得的，^"鵬^'^)與時間的關系的曲 線的線性趨勢的斜率；f為頻率；且V為選定的頻率的速度。
通過返回到步驟1210，可以為每個選定的微震事件重復過程1200，直 到存在Q值的相關統(tǒng)計總體(步驟1260)為止。 一旦Q值的相關的統(tǒng)計總 體被確定，Q因子信息可用于管理井場(步驟1265)或確定潛在的鉆探基 地，之后過程結(jié)束。例如，所述信息可用于指導井場的鉆探操作。統(tǒng)計總 體的均值為被處理的地巖層提供總的Q因子。所述值或值的范圍可用作用 于進一步建模的輸入，或用于將來的感興趣的地巖層的綜合知識。
因此，不同的示例性實施例提供了一種通過使用從在儲油層監(jiān)測期間 記錄的微震數(shù)據(jù)中獲得的S尾波來估計Q因子，以管理井場的方法、儀器 或計算機程序產(chǎn)品。在識別出P和S波的波至之后，從微震事件的軌跡中 提取S尾波窗，并將其用作頻譜分解應用的輸入。當尾波被分解成它的單 頻時，考慮到最大恢復頻率為奈奎斯特頻率，選擇幾個包含在數(shù)據(jù)頻帶內(nèi) 的感興趣的頻率。利用^"鵬^'^與時間關系的曲線的線性趨勢的斜率，
估計出每個單頻的Q因子的值。隨后，對于每個地震檢波器，重復同樣的 操作，并且Q因子與檢波器位置關系的曲線給出隨著深度而感測的衰減 (1/Q)。并且，從Q與頻率關系的曲線的斜率，可確定速度、Q與頻率間 的關系。利用地巖層內(nèi)的幾個微震事件，可以研究Q因子在一個區(qū)域內(nèi)的 變化，以便用對比衰減特性來區(qū)分不同的地帶。并且，當計算頻譜分解時， 不同滑動窗可用于比較Q因子與傳播時間關系的變化，其在觀察間隔內(nèi)的
21衰減變化中可能是有用的。
在井場生產(chǎn)期間，從Q因子確定中獲得的數(shù)據(jù)可用于規(guī)劃地巖層油田 中的鉆探基地?？梢詫χ車牡貛r層和巖石基質(zhì)進行更好的建模，以將關 于地巖層的可能含量和產(chǎn)量值的信息更準確地提供給工程師。通過知道Q 因子，可推斷地巖層中的水飽和度和碳氫化合物濃度。因此，可作出關于 井場處的鉆探和碳氫化合物的產(chǎn)量的更精明的決定。
盡管為了在具體細節(jié)上闡述、解釋和描述某些本發(fā)明的實施例的目的 提供前述內(nèi)容，但是顯然本領域技術(shù)人員可以在不脫離本方面的范圍和精 神的前提下，對所描述的方法、系統(tǒng)和其它實施例作出更改和修改。
2權(quán)利要求
1、一種用于在井場中識別品質(zhì)因子的方法，所述方法包括檢測微震事件，從所述微震事件中識別次級尾波數(shù)據(jù)；以及利用所述次級尾波數(shù)據(jù)識別所述品質(zhì)因子。
2、 如權(quán)利要求l所述的方法，其中，所述方法進一步包括-在第一位置生成所述微震事件，其中，在不同于所述第一位置的第二位置處記錄所述微震事件。
3、 如權(quán)利要求2所述的方法，其中，所述第一位置為井筒，并且所述 微震事件是從來自于所述井筒的水力壓裂生成的。
4、 如權(quán)利要求l所述的方法，其中，用至少一個多分量地震檢波器記 錄所述微震事件。
5、 如權(quán)利要求l所述的方法，其中，所述識別所述次級尾波數(shù)據(jù)的步 驟進一步包括記錄所述微震事件的幅值和所述微震事件的頻率。
6、 如權(quán)利要求l所述的方法，其中，所述次級尾波數(shù)據(jù)包括一個頻率 組，所述方法進一步包括將所述頻率組分解為多個單頻；以及 從所述多個單頻中選擇多個代表性的頻率。
7、 如權(quán)利要求6所述的方法，其中，利用所述地震尾波數(shù)據(jù)來識別所 述品質(zhì)因子的步驟包括為所述多個代表性的頻率中的每一個識別相應的品質(zhì)因子。
8、 如權(quán)利要求1所述的方法，進一步包括從水力壓裂生成所述微震事件。
9、如權(quán)利要求1所述的方法，進一步包括利用所述品質(zhì)因子管理一 組儲油層。
10、 一種用于管理井場的方法，所述方法包括: 檢測微震事件，從所述微震事件中識別次級尾波數(shù)據(jù)； 利用所述次級尾波數(shù)據(jù)識別品質(zhì)因子； 利用所述品質(zhì)因子管理所述井場。
11、 如權(quán)利要求10所述的方法，其中，所述方法進一步包括 在第一位置生成所述微震事件，其中，在不同于所述第一位置的第二位置處記錄所述微震事件。
12、 如權(quán)利要求ll所述的方法，其中，所述微震事件是從來自于井筒 的水力壓裂生成的。
13、 如權(quán)利要求10所述的方法，其中，用至少一個多分量地震檢波器 記錄所述微震事件。
14、 如權(quán)利要求10所述的方法，其中，所述識別所述次級尾波數(shù)據(jù)的 步驟進一步包括記錄所述微震事件的幅值和所述微震事件的頻率。
15、 如權(quán)利要求10所述的方法，其中，所述次級尾波數(shù)據(jù)包括一個頻 率組，所述方法進一步包括將所述頻率組分解為多個頻率段；以及 從所述多個頻率段中選擇多個代表性的頻率。
16、 如權(quán)利要求15所述的方法，其中，所述利用所述次級尾波數(shù)據(jù)識別所述品質(zhì)因子的步驟包括為所述多個代表性的頻率中的每一個識別相 應的品質(zhì)因子。
17、 如權(quán)利要求16所述的方法，進一步包括 從水力壓裂生成所述微震事件。
18、 一種數(shù)據(jù)處理系統(tǒng)包括 總線；連接到所述總線的通信單元；連接到所述總線的存儲設備，其中，所述存儲設備包括計算機可用程 序代碼；以及連接到所述總線的處理器單元，其中，所述處理器單元執(zhí)行所述計算 機可用程序代碼以檢測微震事件，從所述微震事件中識別次級尾波數(shù)據(jù)， 并且利用所述次級尾波數(shù)據(jù)識別品質(zhì)因子。
19、 如權(quán)利要求18所述的數(shù)據(jù)處理系統(tǒng)，其中，所述識別所述次級尾 波數(shù)據(jù)的步驟進一步包括所述處理器單元執(zhí)行所述計算機可用程序代碼， 以記錄所述微震事件的幅值和所述微震事件的頻率。
20、 如權(quán)利要求18所述的數(shù)據(jù)處理系統(tǒng)，其中，所述次級尾波數(shù)據(jù)包 括一個頻率組，并且其中，所述處理器單元進一步執(zhí)行所述計算機可用程 序代碼，以將所述頻率組分解為多個頻率段，并且從所述多個頻率段中選 擇多個代表性的頻率。
21、 如權(quán)利要求20所述的數(shù)據(jù)處理系統(tǒng)，其中，為所述頻率組識別所 述品質(zhì)因子的步驟包括所述處理器單元執(zhí)行所述計算機可用程序代碼， 以為所述多個代表性的頻率中的每一個識別相應的品質(zhì)因子，并且其中， 所述處理器單元進一步執(zhí)行所述計算機可用程序代碼，以分析所述多個單 頻的速度中的變化以確定地巖層和巖石分量中的變化。
22、 一種用于建模井場中的地巖層和巖石基質(zhì)的儀器，所述儀器包括: 用于記錄微震事件的接收器；以及數(shù)據(jù)處理系統(tǒng)，其包括 總線；連接到所述總線的通信單元；連接到所述總線的存儲設備，其中，所述存儲設備包括計算機可 用程序代碼；以及連接到所述總線的處理器單元，其中，所述處理器單元執(zhí)行所述 計算機可用程序代碼以記錄微震事件，識別所述微震事件中的次級尾 波數(shù)據(jù)，并且利用所述次級尾波數(shù)據(jù)識別品質(zhì)因子。
23、 如權(quán)利要求22所述的儀器，其中，所述次級尾波數(shù)據(jù)包括一個頻 率組，并且其中，所述處理器單元進一步執(zhí)行所述計算機可用程序代碼， 以將所述頻率組分解為多個頻率段，以及從所述多個頻率段中選擇多個代 表性的頻率。
24、 如權(quán)利要求23所述的儀器，其中，為所述頻率組識別所述品質(zhì)因 子的步驟包括所述處理器單元執(zhí)行所述計算機可用程序代碼，以為所述 多個代表性的頻率中的每一個識別相應的品質(zhì)因子，并且其中，所述處理 器單元進一步執(zhí)行所述計算機可用程序代碼，以分析所述多個單頻的速度 中的變化以確定地巖層和巖石分量中的變化。
25、 一種計算機程序產(chǎn)品，包括具有計算機可用程序代碼的計算機可讀介質(zhì)，所述計算機可用程序代 碼用于識別品質(zhì)因子，所述計算機程序產(chǎn)品包括-用于記錄微震事件的計算機可用程序代碼；用于識別所述微震事件中的次級尾波數(shù)據(jù)的計算機可用程序代 碼；以及用于利用所述次級尾波數(shù)據(jù)來識別品質(zhì)因子的計算機可用程序代碼。
全文摘要
本發(fā)明利用微震事件生成的S尾波來識別Q因子涉及一種井場中地巖層和巖石基質(zhì)的方法和系統(tǒng)模型。在監(jiān)測井場處記錄井筒中的水力壓裂導致的微震事件。識別微震事件的S尾波窗。隨后識別S尾波窗內(nèi)的一組頻率的Q因子。
文檔編號G01V1/28GK101470211SQ20081019084
公開日2009年7月1日 申請日期2008年12月26日 優(yōu)先權(quán)日2007年12月27日
發(fā)明者J·勒卡爾韋, Y·蘇亞雷斯 申請人:普拉德研究及開發(fā)股份有限公司