新四季網

一種隨壓力變化的孔隙介質橫波速度預測方法

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日
【發明者】李琳, 馬勁風, 王香增, 高瑞民, 江紹靜, 黃春霞, 劉立 申請人:西北大學, 陝西延長石油(集團)有限責任公司研究院

同类文章

一種新型多功能組合攝影箱的製作方法

一種新型多功能組合攝影箱的製作方法【專利摘要】本實用新型公開了一種新型多功能組合攝影箱,包括敞開式箱體和前攝影蓋,在箱體頂部設有移動式光源盒,在箱體底部設有LED脫影板,LED脫影板放置在底板上;移動式光源盒包括上蓋,上蓋內設有光源,上蓋部設有磨沙透光片,磨沙透光片將光源封閉在上蓋內;所述LED脫影

壓縮模式圖樣重疊檢測方法與裝置與流程

本發明涉及通信領域,特別涉及一種壓縮模式圖樣重疊檢測方法與裝置。背景技術:在寬帶碼分多址(WCDMA,WidebandCodeDivisionMultipleAccess)系統頻分復用(FDD,FrequencyDivisionDuplex)模式下,為了進行異頻硬切換、FDD到時分復用(TDD,Ti

個性化檯曆的製作方法

專利名稱::個性化檯曆的製作方法技術領域::本實用新型涉及一種檯曆,尤其涉及一種既顯示月曆、又能插入照片的個性化檯曆,屬於生活文化藝術用品領域。背景技術::公知的立式檯曆每頁皆由月曆和畫面兩部分構成,這兩部分都是事先印刷好,固定而不能更換的。畫面或為風景,或為模特、明星。功能單一局限性較大。特別是畫

一種實現縮放的視頻解碼方法

專利名稱:一種實現縮放的視頻解碼方法技術領域:本發明涉及視頻信號處理領域,特別是一種實現縮放的視頻解碼方法。背景技術: Mpeg標準是由運動圖像專家組(Moving Picture Expert Group,MPEG)開發的用於視頻和音頻壓縮的一系列演進的標準。按照Mpeg標準,視頻圖像壓縮編碼後包

基於加熱模壓的纖維增強PBT複合材料成型工藝的製作方法

本發明涉及一種基於加熱模壓的纖維增強pbt複合材料成型工藝。背景技術:熱塑性複合材料與傳統熱固性複合材料相比其具有較好的韌性和抗衝擊性能,此外其還具有可回收利用等優點。熱塑性塑料在液態時流動能力差,使得其與纖維結合浸潤困難。環狀對苯二甲酸丁二醇酯(cbt)是一種環狀預聚物,該材料力學性能差不適合做纖

一種pe滾塑儲槽的製作方法

專利名稱:一種pe滾塑儲槽的製作方法技術領域:一種PE滾塑儲槽一、 技術領域 本實用新型涉及一種PE滾塑儲槽,主要用於化工、染料、醫藥、農藥、冶金、稀土、機械、電子、電力、環保、紡織、釀造、釀造、食品、給水、排水等行業儲存液體使用。二、 背景技術 目前,化工液體耐腐蝕貯運設備,普遍使用傳統的玻璃鋼容

釘的製作方法

專利名稱:釘的製作方法技術領域:本實用新型涉及一種釘,尤其涉及一種可提供方便拔除的鐵(鋼)釘。背景技術:考慮到廢木材回收後再加工利用作業的方便性與安全性,根據環保規定,廢木材的回收是必須將釘於廢木材上的鐵(鋼)釘拔除。如圖1、圖2所示,目前用以釘入木材的鐵(鋼)釘10主要是在一釘體11的一端形成一尖

直流氧噴裝置的製作方法

專利名稱:直流氧噴裝置的製作方法技術領域:本實用新型涉及ー種醫療器械,具體地說是ー種直流氧噴裝置。背景技術:臨床上的放療過程極易造成患者的局部皮膚損傷和炎症,被稱為「放射性皮炎」。目前對於放射性皮炎的主要治療措施是塗抹藥膏,而放射性皮炎患者多伴有局部疼痛,對於止痛,多是通過ロ服或靜脈注射進行止痛治療

新型熱網閥門操作手輪的製作方法

專利名稱:新型熱網閥門操作手輪的製作方法技術領域:新型熱網閥門操作手輪技術領域:本實用新型涉及一種新型熱網閥門操作手輪,屬於機械領域。背景技術::閥門作為流體控制裝置應用廣泛,手輪傳動的閥門使用比例佔90%以上。國家標準中提及手輪所起作用為傳動功能,不作為閥門的運輸、起吊裝置,不承受軸向力。現有閥門

用來自動讀取管狀容器所載識別碼的裝置的製作方法

專利名稱:用來自動讀取管狀容器所載識別碼的裝置的製作方法背景技術:1-本發明所屬領域本發明涉及一種用來自動讀取管狀容器所載識別碼的裝置,其中的管狀容器被放在循環於配送鏈上的文檔匣或託架裝置中。本發明特別適用於,然而並非僅僅專用於,對引入自動分析系統的血液樣本試管之類的自動識別。本發明還涉及專為實現讀