[發(fā)明專利]一種高含泥質砂巖聲速的計算方法

有效

• 申請?zhí)枺?/td>: 202110289819.4
• 申請日: 2021-03-15
• 公開（公告）號: CN112987096B
• 公開（公告）日: 2021-11-26
• 發(fā)明（設計）人: 韓學輝;江佳洋;韓振鵬;劉智楨;羅振
• 申請（專利權）人: 中國石油大學（華東）;東營市永昇能源科技有限責任公司;克拉瑪依市昂科能源科技有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京精金石知識產權代理有限公司 11470
• 代理人: 尉月麗

1.一種高含泥質砂巖聲速的計算方法，包括以下步驟：

S1、利用樣品薄片鑒定結果，確定高含泥質砂巖中泥質分布方式；

S2、利用染色標記-像素拾取的方法，定量區(qū)分膠結泥質含量V_cementclay；

S3、利用物質平衡方法，計算結構泥質含量V_{structureclay}和分散泥質含量V_{dispersedclay}；

S4、基于修正CCT模型，計算石英和膠結泥質的等效體積模量K_cct和等效剪切模量G_cct；

S5、利用DEM模型將結構泥質替換部分石英作為骨架，計算結構泥質替換骨架后的等效體積模量K_cctdem和等效剪切模量G_cctdem；

S6、基于未固結砂巖模型將分散泥質填充到孔隙中，計算得到結構泥質、膠結泥質和分散泥質約束下的高含泥質砂巖的等效體積彈性模量K_eff和等效剪切彈性模量G_eff；

S7、步驟S6得到的高含泥質砂巖的模量，計算縱波速度V_p1和橫波速度V_s1。

2.根據(jù)權利要求1所述的計算方法，其特征在于，步驟S1中所述高含泥質砂巖的泥質含量為30-50％。

3.根據(jù)權利要求1所述的計算方法，其特征在于，步驟S1中所述泥質分布方式為膠結泥質、分散泥質和結構泥質。

4.根據(jù)權利要求1所述的計算方法，其特征在于，所述步驟S1中還包括測量所述高含泥質砂巖的孔隙度密度ρ、縱波速度V_p、橫波速度V_s和泥質含量V_sh。

5.根據(jù)權利要求1所述的計算方法，其特征在于，步驟S3中所述結構泥質含量V_{structureclay}的計算公式為：

V_g+V_{structureclay}+Φ₀＝1 (1)

其中，V_g為石英骨架含量，V_{structureclay}為結構泥質含量，Φ₀為臨界孔隙度。

6.根據(jù)權利要求1所述的計算方法，其特征在于，步驟S3中所述分散泥質V_{dispersedclay}的計算公式為：

V_sh＝V_{dispersedclay}+V_{structureclay}+V_cementclay (2)

其中，V_sh為總泥質含量。

7.根據(jù)權利要求1所述的計算方法，其特征在于，步驟S4中所述石英和膠結泥質的等效體積模量K_cct和等效剪切模量G_cct的計算公式為：

其中，

膠結物分布模式1：

膠結物分布模式2：

式中，

α₁為膠結平面半徑與砂巖顆粒半徑的比值；

ε是膠結物中心厚度與顆粒半徑的比值；

Φ₀為疏松砂巖的臨界孔隙度；

Φ為疏松砂巖的真實孔隙度；

G_c和G分別為泥質膠結物和石英顆粒的剪切模量；

v_c和v分別為泥質膠結物和石英顆粒的泊松比；

膠結物分布模式1為膠結物僅分布在顆粒之間的情況；

膠結物分布模式2是膠結物均勻分布在顆粒表面時的情況。