用于確定儲層流體識別因子的敏感性的方法及系統(tǒng)的制作方法
【技術領域】
[0001] 本發(fā)明涉及石油地球物理勘探領域，尤其涉及一種用于確定儲層流體識別因子的 敏感性的方法及系統(tǒng)。
【背景技術】
[0002] 地震勘探技術是油氣勘探中應用最為廣泛的一種地球物理學方法，其原理是，利 用地震波在不同介質中傳播的速度、振幅、頻率、相位、波形等參數的變化來分析、預測油氣 儲層分布范圍W及儲層物性特征。目前油氣勘探面臨愈加復雜的情況，比如儲層通常表現(xiàn) 為具有厚度薄、物性高度非均質、有效儲層規(guī)模小、分布分散、巖石物理關系復雜及儲層巖 性差異小等特征。對于送些復雜情況，如果僅僅利用少量的或者單一的信息難W達到區(qū)分 識別儲層及流體的目的。
[0003] 為了從地震資料中獲得更多的反映儲層流體的信息，國內外很多學者已經提出了 多種不同的流體識別因子W及相應的流體區(qū)分方法。如化tti等定義的與縱橫波反射系數 有關的流體識別因子Δ F，Goodway等（1997)提出的識別流體異常的LMR法，Russell等提 出的Pf流體識別因子，Dillon等采用的直接姪類指示值HI)的波阻抗差分法、截距斜率 法（PG剖面）等。國內的比如有學者提出了具有波阻抗量綱的一次方和二次方的流體識別 因子 FIFI 和 FSFIF。
[0004] 在儲層預測及流體識別中，上述各種流體識別因子和流體識別方法均在某一方面 有較強的識別能力，但不能全面有效的識別不同類型的儲層及流體，不同的流體識別因子 對儲層流體的敏感性不同。如何在繁多的流體識別因子中選擇合適的、對儲層識別能力高 的流體識別因子，到目前為止仍缺乏理論指導，制約了油氣勘探工作的進一步開展。
[0005] 綜上，如何在眾多流體識別因子中優(yōu)選出對儲層流體最為敏感的流體識別因子， 成為亟待解決的問題。

【發(fā)明內容】

[0006] 本發(fā)明所要解決的技術問題之一是需要提供一種在眾多流體識別因子中優(yōu)選出 對儲層流體最為敏感的流體識別因子的方法。
[0007] 為了解決上述技術問題，本發(fā)明提供了一種用于確定儲層流體識別因子的敏感性 的方法，包括W下步驟：利用Biot-Gassmann方程進行流體替換來獲得儲層段完全飽水狀 態(tài)和完全飽氣狀態(tài)時的物性參數；基于儲層段完全飽水狀態(tài)和完全飽氣狀態(tài)時的物性參數 得到兩種狀態(tài)的多種流體識別因子，并基于所述多種流體識別因子構建多個敏感性定量因 子；基于各個敏感性定量因子的數值的大小確定流體識別因子對儲層流體的敏感性的大 小。
[0008] 在一個實施例中，所述流體替換包括；根據儲層測井數據中的飽含流體巖 石的縱波速度、橫波速度W及密度得到飽含流體巖石的體積模量和剪切模量；利用 Biot-Gassmann方程，根據飽含流體巖石的體積模量和剪切模量W及孔隙所飽含的流體的 體積模量得到干巖石的體積模量和剪切模量；利用Biot-Gassmann方程，根據干巖石的體 積模量和剪切模量、W及巖石完全飽水和完全飽氣兩種狀態(tài)下的孔隙所飽含的流體的體積 模量，分別得到巖石在完全飽水和完全飽氣時的體積模量和剪切模量；分別利用巖石在完 全飽水和完全飽氣時的體積模量、剪切模量W及密度得到完全飽水巖石和完全飽氣巖石的 縱波速度和橫波速度。
[0009] 在一個實施例中，通過W下表達式來構建所述敏感性定量因子t :
[0010]
[0011] 式中，和Fg分別為儲層巖石完全飽含水和完全飽含氣時的流體識別因子； max化，F(xiàn)g)表示取和Fg中的較大值，min化，F(xiàn)g)表示取和Fg中的較小值。
[0012] 在一個實施例中，在確定流體識別因子對儲層流體的敏感性的大小的步驟中，包 括：對于數值大的敏感性定量因子，確定與該敏感性定量因子對應的流體識別因子對儲層 流體的敏感性大；對于數值小的敏感性定量因子，確定與該敏感性定量因子對應的流體識 別因子對儲層流體的敏感性小。
[0013] 在一個實施例中，當儲層流體為油水混合態(tài)時，通過W下表達式來構建所述敏感 性定量因子t:
[0014]
[0015] 式中，F(xiàn)。為儲層巖石完全飽含油時的流體識別因子。
[0016] 另一方面，還提供了一種用于確定儲層流體識別因子的敏感性的系統(tǒng)，包括W下 模塊：流體替換模塊，其利用Biot-Gassmann方程進行流體替換來獲得儲層段完全飽水狀 態(tài)和完全飽氣狀態(tài)時的物性參數；敏感性定量因子的建立模塊，其基于儲層段完全飽水狀 態(tài)和完全飽氣狀態(tài)時的物性參數得到兩種狀態(tài)的多種流體識別因子，并基于所述多種流體 識別因子構建多個敏感性定量因子；敏感性分析模塊，其基于各個敏感性定量因子的數值 的大小確定流體識別因子對儲層流體的敏感性的大小。
[0017] 在一個實施例中，流體替換模塊包括：第一彈性模量計算單元，其根據儲層測井數 據中的飽含流體巖石的縱波速度、橫波速度W及密度得到飽含流體巖石的體積模量和剪切 模量；第二彈性模量計算單元，其利用Biot-Gassmann方程，根據飽含流體巖石的體積模量 和剪切模量W及孔隙所飽含的流體的體積模量得到干巖石的體積模量和剪切模量；第Η彈 性模量計算單元，其利用Biot-Gassmann方程，根據干巖石的體積模量和剪切模量、W及巖 石完全飽水和完全飽氣兩種狀態(tài)下的孔隙所飽含的流體的體積模量，分別得到巖石在完全 飽水和完全飽氣時的體積模量和剪切模量；物性參數計算單元，其分別利用巖石在完全飽 水和完全飽氣時的體積模量、剪切模量W及密度得到完全飽水巖石和完全飽氣巖石的縱波 速度和橫波速度。
[0018] 在一個實施例中，敏感性定量因子的建立模塊通過W下表達式來構建所述敏感性 定量因子t:
[0019]
[0020] 式中，和Fg分別為儲層巖石完全飽含水和完全飽含氣時的流體識別因子； max化，F(xiàn)g)表示取和Fg中的較大值，min化，F(xiàn)g)表示取和Fg中的較小值。
[0021] 在一個實施例中，敏感性分析模塊在確定流體識別因子對儲層流體的敏感性的大 小時，對于數值大的敏感性定量因子，確定與該敏感性定量因子對應的流體識別因子對儲 層流體的敏感性大；對于數值小的敏感性定量因子，確定與該敏感性定量因子對應的流體 識別因子對儲層流體的敏感性小。
[0022] 在一個實施例中，當儲層流體為油水混合態(tài)時，所述敏感性定量因子的建立模塊 通過W下表達式來構建所述敏感性定量因子t:
[0023]
[0024] 式中，F(xiàn)。為儲層巖石完全飽含油時的流體識別因子。
[00巧]與現(xiàn)有技術相比，上述方案中的一個或多個實施例可W具有如下優(yōu)點或有益效 果：
[0026] 通過建立敏感性定量因子，實現(xiàn)了對各種流體識別因子對儲層流體的敏感性的定 量描述，進一步根據敏感性定量因子對流體識別因子進行優(yōu)選，不但能為測井資料解釋提 供重要的參考，還能有效地指導地震屬性參數的提取，進而可為儲層預測與流體識別服務。 本方法具有計算速度快、穩(wěn)定性好的優(yōu)點，可直接用于地震資料的解釋工作。
[0027] 本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且部分地從說明書中變得 顯而易見，或者通過實施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點可通過在說明書、權利要 求書W及附圖中所特別指出的結構來實現(xiàn)和獲得。
【附圖說明】
[0028] 附圖用來提供對本發(fā)明的進一步理解，并且構成說明書的一部分，與本發(fā)明的實 施例共同用于解釋本發(fā)明，并不構成對本發(fā)明的限制。在附圖中：
[0029] 圖1是根據本發(fā)明第一實施例的用于確定儲層流體識別因子的敏感性的方法的 流程圖；
[0030] 圖2是第一類砂巖的流體識別因子的敏感性定量因子的柱狀圖；
[0031] 圖3是第二類砂巖的流體識別因子的敏感性定量因子的柱狀圖；
[0032] 圖4是第Η類砂巖的流體識別因子的敏感性定量因子的柱狀圖；
[0033] 圖5是某區(qū)域的實際測井曲線圖；
[0034] 圖6是某區(qū)域實際測井資料某段的流體識別因子的敏感性定量因子的柱狀圖；
[0035] 圖7是根據表1的11號流體識別因子的縱向剖面圖；
[0036] 圖8是根據表1的11號流體識別因子的沿層切片圖；
[0037] 圖9是根據本發(fā)明第二實施例的用于確定儲層流體識別因子的敏感性的系統(tǒng)的 結構示意圖。
【具體實施方式】
[0038] W下將結合附圖及實施例來詳細說明本發(fā)明的實施方式，借此對本發(fā)明如何應用 技術手段來解決技術問題，并達成技術效果的實現(xiàn)過程能充分理解并據W實施。需要說明 的是，只要不構成沖突，本發(fā)明中的各個實施例w及各實施例中的各個特征可w相互結合， 所形成的技術方案均在本發(fā)明的保護范圍之內。
[0039] 第一連施例
[0040] 圖1是根據本發(fā)明第一實施例的用于確定儲層流體識別因子的敏感性的方法的 流程圖。下面參考圖1來詳細說明本方法的各個步驟。
[0041] 步驟S110、利用Biot-Gassmann方程進行流體替換來獲得儲層段完全飽水