波阻抗剖面分析研究成果摘要

波阻抗剖面技術是地震巖性勘探中的重要方法之一。它是將地震資料與測井資用線性正、反演技術求取地下地層的反射系數(shù)、波阻抗、擬速度等重要于波阻抗剖面可以提供眾多地下信息，故該方法已被地球物理學界所注自建的PDP－11－45機的波阻抗剖面系列，在合成記錄的研究中，提供的H一是一個成功的方法，達到了引進的CYBER一720所做的合成記錄水平；在常型范疇內實現(xiàn)了在時間域中求取反射系數(shù)、相對波阻和絕對波阻的方法H－T0轉換方法和最佳反褶積算子求取法的成功，保證了地層反褶積取其誤差約為1％。這些試驗結果，對實際處理是有意義的。為找油氣、和巖石性質，有較大的實用價值。 2100433B

地震波阻抗反演是儲層預測的有效手段之一，波阻抗與含油氣儲層有很好的對應性。

波阻抗反演實際上是從地震剖面上消除子波影響,留下反射系數(shù),再由反射系數(shù)計算出能反映地層物性變化的物理參數(shù)波阻抗。

波阻抗反演的方法很多，其中最常用的是廣義線性反演，該方法有較高的反演精度，但是在反演過程中，高頻噪聲對解病態(tài)方程組影響很大,即使加了正則化因子,也無法消除高頻噪聲的影響。

根據(jù)反射系數(shù)序列，并給定波阻抗初始值，利用遞推波阻抗的反演方法就可以算出相應的波阻抗序列。該方法的核心是反演反子波，但是反演出反子波局部化特征并不理想，因為兩翼都有振蕩，并不是標準的δ(t)函數(shù)，而是帶有剩余子波的δ(t)函數(shù)，其頻率過低，計算出的反射系數(shù)頻率也較低，相當于帶有剩余子波的反射系數(shù)。遞推波阻抗的反演方法對誤差及初值 Z0 均十分敏感，穩(wěn)定性也差?？梢杂眠f推波阻抗反演法反演出的波阻抗作為迭代反演波阻抗的初值。

因此，針對廣義線性反演，重點討論了迭代解方程的波阻抗算法、解方程的正則化方法、波阻抗反演中的分辨率、改善解的二次圓滑濾波、用于控制反演質量的統(tǒng)計量和均方差及波阻抗的約束反演等關鍵技術問題;針對遞推波阻抗的反演,重點討論了遞推波阻抗的反演方法、子波和反子波提取、誤差對提取反子波的影響等關鍵技術問題。同時做了大量的模型試算。2100433B

1 波阻抗計算

對有一定預應力的巖石進行循環(huán)沖擊試驗的裝置為自行研制的動靜組合加載裝置。對于有軸壓時的沖擊過程中，桿件與試樣仍然滿足裝置賴以存在的一維應力波傳輸理論。應力波從一種介質傳播到另一種波阻抗不同的介質時，入射波在2 種介質的界面處會產(chǎn)生反射波和透射波。由于其值比巖石的波阻抗大許多，認為在整個沖擊過程中其大小保持恒定；

在整個沖擊過程中并不是每一瞬時的反射波或透射波都夾雜多次反射和透射波的影響，當時，測試所得的反射波中沒有 2 次和更高次反射的影響。同理，在透射波采樣時間內的區(qū)段，試驗所得的透射波中也沒有夾雜多次透射波的影響。因此，可以利用SHPB 試驗系統(tǒng)采集的入射波和反射波或透射波在各自采樣區(qū)段內相對應的值計算巖石的波阻抗。而其他時間區(qū)段的波阻抗不可由試驗數(shù)據(jù)直接計算。

2 計算波阻抗時計算式的選擇

在實際試驗時，反射波和入射波是用入射桿上的同一個應變片測量的，特別是有軸壓的沖擊試驗時，反射波的起跳點受入射波的影響較大，其對應的時間點和大小不易確定，極易引起誤判。而透射波不在這種情況，起跳點相對容易判斷，因此計算巖石的波阻抗較為理想。

由于巖石與彈性桿間的波阻抗大小不同，每一入射波經(jīng)過透反射后必有一反射波和透射波與之相對應；沖擊過程中的不同時刻巖石所受的沖擊動載大小不同，巖石試件內部裂隙的閉合程度或擴展程度亦不相同，即巖石試件具有不同波阻抗。換言之，在一次沖擊過程中，入射波和反射波或透射波與巖石的波阻抗一一對應，不同時刻波阻抗的大小很可能不一樣，入射應力大小、透射應力大小和時間點3 個變量必須確定其中2 個不變。根據(jù)一維應力波理論可知，在異型沖頭形狀、長度和沖擊速度相同的前提下，在區(qū)間固定一時間點，則每次沖擊得到的理論上應為定值，對照時間確定出對應的透射波的值。具體選其中哪一個時間點視情況而定，它們之間的差別微乎其微。要強調的是對同一巖石試件的循環(huán)沖擊過程的損傷進行分析時，應選擇同一時間點。

另外，利用聲波定義循環(huán)沖擊過程巖石的累積損傷，也是在沖擊載荷作用后測量巖石的波速值。本文的方法與此相似，每次測量（除第1 次）的值都是前1 次沖擊作用后巖石在壓密后的波阻抗值，這個值表征巖石在前1 次沖擊后內部的累積損傷，因此，在測量程序上也與聲波定義損傷的一致 。