一種隨壓力變化的孔隙介質橫波速度預測方法
2023-11-02 03:20:57
一種隨壓力變化的孔隙介質橫波速度預測方法
【專利摘要】一種隨壓力變化的孔隙介質橫波速度預測方法,對測量配位數進行擬合併且加權,得到配位數Cp,將配位數Cp代入考慮壓力變化的Digby公式得到Kdry,將配位數Cp代入考慮壓力變化的Mindlin公式中得到μdry,將Kdry與μdry代入Gassmann方程的變形公式中計算預測縱波速度由此得到加權係數W;將加權係數W代入公式(2),得到配位數Cp,得到Cp之後,把Cp值代入公式(8)中得到幹巖石的切變模量μdry,進而將μdry代入Gassmann方程的變形公式中計算得到預測橫波速度依據隨壓力變化的橫波速度可以建立包括四維AVO模型及彈性阻抗模型等,進而預測開發階段儲層壓力變化的油氣藏屬性;這種預測得到的隨壓力變化的橫波速度更加符合實際儲層變化情況。
【專利說明】一種隨壓力變化的孔隙介質橫波速度預測方法
【技術領域】
[0001]本發明屬於地震勘探中巖石物理領域,具體涉及一種隨壓力變化的橫波速度預測方法。
【背景技術】
[0002]通過AVO分析,地球物理學家可以更好地評估油氣藏巖石屬性,包括孔隙度、密度、巖性與流體屬性,而橫波速度是建立AVO模型、轉換波分析過程中的一個不可缺少的彈性參數。多數情況下研究工區是沒有橫波速度測井資料的;人們通常用縱橫波速度比為常數來代替橫波速度,但是對於不同的介質縱橫波速度比應該是不同的,這樣的假設不合理。因此很多的地球物理學者都在研究橫波速度的預測方法。有經驗公式,也有基於巖石物理理論的。
[0003]Pickett於1963年給出了灰巖的縱橫波關係Vs=Vp/l.9,對於白雲巖他則給出了 Vs=Vp/l.8。Castagna等人於1993年將這個公式進行了修改,灰巖為Vs=-0.055Vp2+l.017Vp-l.031,而白雲巖則為Vs=0.583Vp-0.078,他同時還提出了碎屑巖的縱橫波速度關係式為Vs=0.804Vp-0.856。
[0004]比較著名的經驗公式包括Castagna等人在1985年提出的著名的泥巖線為Vs=0.862Vp-l.172。Gardner於1974年給出了不同的巖性之間的速度與密度的關係,其中它的平均變換式為P =0.23V°_25,這個平均變換式是對所有巖性的速度與密度關係的最佳擬合,它適合於所有巖性,不僅僅適用用某種巖性。而Castagna於1993年又對Gardner的
公式進行了擴充,得到了不同巖性的速度與密度之間的關係:對於砂巖有P = MwO2'對於頁巖有Z7 = ^K2b5,對於石灰巖有P = a243C225,對於白雲巖有P = 0.226V;243,
對硬石膏則是廣= (U)fK)C'。Wyllie等人在1958年和1963年陸續提出了充滿鹽水的孔`隙介質的孔隙度與速度之間的經驗關係:l/V=(l-(j5)/VM+(j5/Vfl,其中V為巖石的整體速度,Vma為巖石骨架的速度,Vfl為孔隙流體的速度,Φ為孔隙度。這個公式還通常可以寫成層間旅行時的表達式:Λ?=(1_Φ) Λ--+Φ Atfl,其中At代表整個巖層的旅行時,Atma為骨架的旅行時,而Atfl為孔隙流體的旅行時,而Wyllie的這個時間平均公式還包括許多假設和限制,如:這個方程要用於孔隙流體是鹽水的情形,用於深度小於2700米的巖石,而且這個巖石的膠結程度和固結程度要很好,並且孔隙度為中等。當某些測井曲線缺失的時候,或者地震振幅異常都可以應用這些期限進行質量監控,但是這些公式對於巖性的依賴很強,而且依賴於局部條件,而Mavko等人多次在他們的書中提到:「這些關係公式都是經驗公式,因而嚴格來說它們只能用在當時研究的巖石上」,因此這些經驗公式不具有普遍性。
[0005]隨著巖石物理理論的完善,基於巖石物理理論的橫波速度預測方法逐漸成為研究的主流。如 Greenberg 和 Castagna 於 1992 年利用 Biot-Gassmann Theory (BGT)進行橫波速度的預測,也就是假設在縱橫波速度之間存在一個穩固的關係同時假設固體巖石成分之間的混合定律是線性的。J0rstad等人於1999年根據基於內含物的有效介質理論進行
橫波速度的預測,並且得出結論,即使有效介質理論比回歸統計方法更加複雜,但是它仍然有優勢因為它可以將泥巖的影響以及孔隙形狀的影響體現在公式中。許多地球物理學家都喜歡用Gassmann方程進程橫波速度的預測,這是因為Gassmann方程的大部分參數如顆粒的體變模量Kma與切變模量μ ma等都是很好獲得的,因此很多方法的給出無論是對於砂巖還是對碳酸鹽都是基於Gassmann方程的。不過在Gassmann方程中幹巖石的體變模量與切變模量是個很難解決的問題,因此很多的地球物理學家給出了幹巖石體變模量與切變模量的計算方法,知道了幹巖石的體變模量與切變模量,巖石的縱橫波速度就很好獲得。Xu和White將Kuster和Toks5z於1974年建立的理論與差分有效介質理論結合,進行巖石彈性模量的計算,具體表現為利用孔隙縱橫比來表徵砂泥成分之間的關係。Nolen-Hoeksema與Wang Zhijing於1996年根據實驗室測得的幹巖石的橫波速度,利用Gassmann方程計算出幹巖石的彈性模量,進而用到流體飽和巖石的橫波速度預測中。在2006年的時候Lee提出的用固結參數聯繫基質彈性模量與骨架彈性模量之間的關係,通過實測縱波速度與預測縱波速度的比較得出固結參數,然後利用固結參數計算橫波速度。2008年孫福利等人利用實際的數據對Lee的方法進行了驗證,並且提出了固結參數的取值範圍。
[0006]但是這些方法均未考慮到壓力對介質屬性的影響。比如在實際生產中,隨著油田開發的進行,無論是注水還是注氣,儲層壓力是會發生變化的。而在CCS(碳捕捉與封存)技術中,我們知道在CO2注入地下及C02-E0R的過程中,注入井點的壓力大而生產井點的壓力較小。隨著CO2的不斷注入,包括在CO2地質封存的不同階段(注入過程中、注入完畢和注入完成後相當長的時間內),儲層內孔隙壓力會發生很大變化,會使得差異壓力發生變化,從而幹巖石的體變模量和切變模量也發生了變化,那麼縱橫波速度也會隨著發生變化。利用四維地震監測CO2在地下封存的狀態過程中,無論是四維地震解釋還是四維AVO反演、彈性阻抗反演,都需要利用隨壓力變化的橫波速度才能進行。
【發明內容】
·
[0007]本發明的目的在於提供一種隨壓力變化的孔隙介質橫波速度預測的方法,利用採集的數據進行儲層橫波速度的預測,預測的橫波速度更加符合實際情況。
[0008]為實現上述目的,本發明採用如下的技術方案:
[0009]本發明包括以下步驟:
[0010]I)採集數據:採集巖石的孔隙度Φ,巖石的體積密度P,流體的體變模量Kf,實際縱波速度Vp_SUMd,巖石骨架的體變模量Kma,巖石骨架的切變模量μ ma,差異壓力P,測量配位數Cp』,巖石顆粒變形之前接觸區域的半徑a與巖石顆粒的半徑R ;
[0011]對測量配位數Cp』進行加權,得到配位數Cp的公式(2),式中W為加權係數:
[0012]Cp=W*Cp』(2)
[0013]2)利用配位數Cp以及Digby公式得到幹巖石的體變模量Kty,利用配位數Cp以及Mindlin公式得到幹巖石的切變模量μ ,然後根據得到的幹巖石的體變模量Ktey、幹巖石的切變模量μ dry以及Gassmann方程的變形公式,得到含有加權係數W的預測縱波速度Vpsm,根據預測的縱波速度等於測量的實際縱波速度Vpnreasural,得到加權係數W ;[0014]3)根據加權係數W以及Mindlin公式,得到幹巖石的切變模量μ ,將幹巖石的切變模量μ dry代入Gassmann方程的變形公式中得到預測橫波速度^
[0015]4)依據預測橫波速應K 建立四維AVO模型及彈性阻抗模型,預測開發階段儲層壓力變化的油氣藏屬性。
[0016]所述步驟I)中測量配位數Cp』是通過以下過程得到的:對C;與eH進行線性擬合得到測量配位數與孔隙度的關係,如公式(I)所示:
[0017]C; =11.759θ1_φ-12.748(I)。
[0018]所述步驟2)中加權係數W是通過以下過程得到的:
[0019]利用Gassmann方程的變形公式進行縱橫波速度的預測,公式(3)-(5)為Gassmann方程的變形公式:
【權利要求】
1.一種隨壓力變化的孔隙介質橫波速度預測方法,其特徵在於,包括以下步驟: 1)採集數據:採集巖石的孔隙度Φ,巖石的體積密度P,流體的體變模量Kf,實際縱波速度Vp_SUMd,巖石骨架的體變模量Kma,巖石骨架的切變模量μ ma,差異壓力P,測量配位數cp』,巖石顆粒變形之前接觸區域的半徑a與巖石顆粒的半徑R ; 對測量配位數Cp』進行加權,得到配位數Cp的公式(2),式中W為加權係數: Cp=ff*C;(2) 2)利用配位數Cp以及Digby公式得到幹巖石的體變模量Kty,利用配位數Cp以及Mindlin公式得到幹巖石的切變模量μ dl7,然後根據得到的幹巖石的體變模量Kdl7、幹巖石的切變模量μ 以及Gassmann方程的變形公式,得到含有加權係數W的預測縱波速度,根據預測的縱波速度等於測量的實際縱波速度Vpnreasural,得到加權係數W ; 3)根據加權係數W以及Mindlin公式,得到幹巖石的切變模量μ,將幹巖石的切變模量μ dry代入Gassmann方程的變形公式中得到預測橫波速度K、; 4)依據預測橫波速度,建立四維AVO模型及彈性阻抗模型,預測開發階段儲層壓力變化的油氣藏屬性。
2.根據權利要求1所述的一種隨壓力變化的孔隙介質橫波速度預測方法,其特徵在於,所述步驟I)中測量配位數Cp』是通過以下過程得到的:對Cp』與e1—1*.進行線性擬合得到測量配位數與孔隙度的關係,如公式(I)所示: Cp' =11.759θ1_φ-12.748(I)。
3.根據權利要求1所述的一種隨`壓力變化的孔隙介質橫波速度預測方法,其特徵在於,所述步驟2)中加權係數W是通過以下過程得到的: 利用Gassmann方程的變形公式進行縱橫波速度的預測,公式(3)- (5)為Gassmann方程的變形公式:
(1-^)2
K / 4
K H--"Ui__μ η
1^drv 丁 / λ I V-> ^ sat
Φ , '-Φ3VJ^α?__π; Ilρ y sat= y dry⑷ V=Vrl = I— = 1-^z- 一 \丨 P ? PCS; 其中與分別為預測的縱波速度、橫波速度,Usat為孔隙介質的切變模量,Uty為幹巖石的切變模量;Kdl7為幹巖石的體變模量,Kma為巖石骨架的體變模量;Φ為巖石的孔隙度;Kf為流體的體變模量,P為巖石的體積密度; 將配位數Cp代入Digby公式中得到:舊(7) … 3kR(\-v)將配位數Cp代入Mindlin公式中得到:
4.根據權利要求1或3所述的一種隨壓力變化的孔隙介質橫波速度預測方法,其特徵在於,所述Mindlin公式為:
5.根據權利要求4所述的一種隨壓力變化的孔隙介質橫波速度預測方法,其特徵在於,所述Mindlin公式(9)中
6.根據權利要求1或3所述的一種隨壓力變化的孔隙介質橫波速度預測方法,其特徵在於,所述Digby公式為 :
【文檔編號】G01V1/30GK103576196SQ201310521170
【公開日】2014年2月12日 申請日期:2013年10月28日 優先權日:2013年10月28日
【發明者】李琳, 馬勁風, 王香增, 高瑞民, 江紹靜, 黃春霞, 劉立 申請人:西北大學, 陝西延長石油(集團)有限責任公司研究院