海洋学报

序号	参数名称	Wood 模型	Biot-Stoll 模型	Dvorkin等效介质模型
注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。
1	孔隙流体密度ρ_f/（kg·m⁻³）	1 025	1 025	1 025
2	基质颗粒密度ρ_g/（kg·m⁻³）	2 650	2 650	2 650
3	绝对黏度η/（kg·m⁻³）	−	0.001	−
4	密度ρ/（kg·m⁻³）	−	实际	实际
5	孔隙度ϕ/%	实际	实际	实际
6	临界孔隙度ϕ_c/%	−	−	38
7	基质颗粒体积模量K_r/Pa	1.47×10¹⁰	1.47×10¹⁰	1.47×10¹⁰
8	基质颗粒剪切模量µ_r/Pa	−	−	1.34×10¹⁰
9	孔隙流体体积模量K_f/Pa	2.18×10⁹	2.18×10⁹	2.18×10⁹
10	泊松比σ	−	0.15	0.15
11	曲折度c	−	1.35	−
12	渗透率κ/m²	−	公式（13）	−
13	体积衰减系数δ_b	−	0.1	−
14	剪切衰减系数δ_s	−	0.1	−
15	颗粒平均接触次数n	−	−	8.5

序号

参数名称

Wood
模型

Biot-Stoll
模型

Dvorkin等效
介质模型

注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。

孔隙流体密度ρ_f/（kg·m⁻³）

1 025

基质颗粒密度ρ_g/（kg·m⁻³）

2 650

绝对黏度η/（kg·m⁻³）

−

0.001

−

密度ρ/（kg·m⁻³）

−

实际

孔隙度ϕ/%

实际

临界孔隙度ϕ_c/%

−

基质颗粒体积模量K_r/Pa

1.47×10¹⁰

基质颗粒剪切模量µ_r/Pa

−

1.34×10¹⁰

孔隙流体体积模量K_f/Pa

2.18×10⁹

泊松比σ

−

0.15

曲折度c

−

1.35

−

渗透率κ/m²

−

公式（13）

−

体积衰减系数δ_b

−

0.1

−

剪切衰减系数δ_s

−

0.1

−

颗粒平均接触次数n

−

8.5

序号	参数名称	Wood 模型	Biot-Stoll 模型	Dvorkin等效介质模型
注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。
1	孔隙流体密度ρ_f/（kg·m⁻³）	1 025	1 025	1 025
2	基质颗粒密度ρ_g/（kg·m⁻³）	2 650	2 650	2 650
3	绝对黏度η/（kg·m⁻³）	−	0.001	−
4	密度ρ/（kg·m⁻³）	−	实际	实际
5	孔隙度ϕ/%	实际	实际	实际
6	临界孔隙度ϕ_c/%	−	−	38
7	基质颗粒体积模量K_r/Pa	1.47×10¹⁰	1.47×10¹⁰	1.47×10¹⁰
8	基质颗粒剪切模量µ_r/Pa	−	−	1.34×10¹⁰
9	孔隙流体体积模量K_f/Pa	2.18×10⁹	2.18×10⁹	2.18×10⁹
10	泊松比σ	−	0.15	0.15
11	曲折度c	−	1.35	−
12	渗透率κ/m²	−	公式（13）	−
13	体积衰减系数δ_b	−	0.1	−
14	剪切衰减系数δ_s	−	0.1	−
15	颗粒平均接触次数n	−	−	8.5

序号

参数名称

Wood
模型

Biot-Stoll
模型

Dvorkin等效
介质模型

注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。

孔隙流体密度ρ_f/（kg·m⁻³）

1 025

基质颗粒密度ρ_g/（kg·m⁻³）

2 650

绝对黏度η/（kg·m⁻³）

−

0.001

−

密度ρ/（kg·m⁻³）

−

实际

孔隙度ϕ/%

实际

临界孔隙度ϕ_c/%

−

基质颗粒体积模量K_r/Pa

1.47×10¹⁰

基质颗粒剪切模量µ_r/Pa

−

1.34×10¹⁰

孔隙流体体积模量K_f/Pa

2.18×10⁹

泊松比σ

−

0.15

曲折度c

−

1.35

−

渗透率κ/m²

−

公式（13）

−

体积衰减系数δ_b

−

0.1

−

剪切衰减系数δ_s

−

0.1

−

颗粒平均接触次数n

−

8.5

序号	密度/（g·cm⁻³）	孔隙度/%	孔隙比	水质量/g	砂土质量/g	黏土质量/g
1	1.54	51.5	1.06	3 061	11 130	7 420
2	1.57	50.5	1.02	3 122	11 350	7 567
3	1.61	49.2	0.97	3 201	11 638	7 759
4	1.64	48.5	0.94	3 251	11 818	7 879
5	1.67	47.4	0.9	3 320	12 067	8 045
6	1.71	46.2	0.86	3 391	12 327	8 218
7	1.75	44.8	0.81	3 485	12 667	8 445
8	1.8	43.2	0.76	3 583	13 027	8 685
9	1.84	42.5	0.74	3 625	13 177	8 785
10	1.88	40.8	0.69	3 732	13 567	9 045

序号

密度/（g·cm⁻³）

孔隙度/%

孔隙比

水质量/g

砂土质量/g

黏土质量/g

1.54

51.5

1.06

3 061

11 130

7 420

1.57

50.5

1.02

3 122

11 350

7 567

1.61

49.2

0.97

3 201

11 638

7 759

1.64

48.5

0.94

3 251

11 818

7 879

1.67

47.4

0.9

3 320

12 067

8 045

1.71

46.2

0.86

3 391

12 327

8 218

1.75

44.8

0.81

3 485

12 667

8 445

1.8

43.2

0.76

3 583

13 027

8 685

1.84

42.5

0.74

3 625

13 177

8 785

1.88

40.8

0.69

3 732

13 567

9 045

序号	密度/（g·cm⁻³）	孔隙度/%	孔隙比	水质量/g	砂土质量/g	黏土质量/g
1	1.54	51.5	1.06	3 061	11 130	7 420
2	1.57	50.5	1.02	3 122	11 350	7 567
3	1.61	49.2	0.97	3 201	11 638	7 759
4	1.64	48.5	0.94	3 251	11 818	7 879
5	1.67	47.4	0.9	3 320	12 067	8 045
6	1.71	46.2	0.86	3 391	12 327	8 218
7	1.75	44.8	0.81	3 485	12 667	8 445
8	1.8	43.2	0.76	3 583	13 027	8 685
9	1.84	42.5	0.74	3 625	13 177	8 785
10	1.88	40.8	0.69	3 732	13 567	9 045

序号

密度/（g·cm⁻³）

孔隙度/%

孔隙比

水质量/g

砂土质量/g

黏土质量/g

1.54

51.5

1.06

3 061

11 130

7 420

1.57

50.5

1.02

3 122

11 350

7 567

1.61

49.2

0.97

3 201

11 638

7 759

1.64

48.5

0.94

3 251

11 818

7 879

1.67

47.4

0.9

3 320

12 067

8 045

1.71

46.2

0.86

3 391

12 327

8 218

1.75

44.8

0.81

3 485

12 667

8 445

1.8

43.2

0.76

3 583

13 027

8 685

1.84

42.5

0.74

3 625

13 177

8 785

1.88

40.8

0.69

3 732

13 567

9 045

序号	参数	函数表达式	R²	RMSE
注：纵波速度单位为km/s，密度单位为g/cm³。
1	密度ρ	ρ=0.669V_p^2.699	0.943 1	0.029 19
2	孔隙度ϕ	ϕ=1.379V_p^−3.161	0.959 4	0.007 65

序号

参数

函数表达式

R²

RMSE

注：纵波速度单位为km/s，密度单位为g/cm³。

密度ρ

ρ=0.669V_p^2.699

0.943 1

0.029 19

孔隙度ϕ

ϕ=1.379V_p^−3.161

0.959 4

0.007 65

序号	参数	函数表达式	R²	RMSE
注：纵波速度单位为km/s，密度单位为g/cm³。
1	密度ρ	ρ=0.669V_p^2.699	0.943 1	0.029 19
2	孔隙度ϕ	ϕ=1.379V_p^−3.161	0.959 4	0.007 65

序号

参数

函数表达式

R²

RMSE

注：纵波速度单位为km/s，密度单位为g/cm³。

密度ρ

ρ=0.669V_p^2.699

0.943 1

0.029 19

孔隙度ϕ

ϕ=1.379V_p^−3.161

0.959 4

0.007 65

参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。
抗剪强度τ	τ=−61.22ρ²+259.2ρ−213.5	0.946 4	1.562
	τ=−555.3ϕ²+351.3ϕ+6.81	0.958 0	1.383
黏聚力c	c=−70.15ρ²+297.1ρ−256.2	0.924 1	2.161
	c=−637.2ϕ²+403.7ϕ−3.845	0.929 1	2.089
内摩擦角φ	φ=6.899ρ+7.547	0.565 1	0.740
	φ=−21.92ϕ+29.47	0.550 4	0.753

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。

抗剪强度τ

τ=−61.22ρ²+259.2ρ−213.5

0.946 4

1.562

τ=−555.3ϕ²+351.3ϕ+6.81

0.958 0

1.383

黏聚力c

c=−70.15ρ²+297.1ρ−256.2

0.924 1

2.161

c=−637.2ϕ²+403.7ϕ−3.845

0.929 1

2.089

内摩擦角φ

φ=6.899ρ+7.547

0.565 1

0.740

φ=−21.92ϕ+29.47

0.550 4

0.753

参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。
抗剪强度τ	τ=−61.22ρ²+259.2ρ−213.5	0.946 4	1.562
	τ=−555.3ϕ²+351.3ϕ+6.81	0.958 0	1.383
黏聚力c	c=−70.15ρ²+297.1ρ−256.2	0.924 1	2.161
	c=−637.2ϕ²+403.7ϕ−3.845	0.929 1	2.089
内摩擦角φ	φ=6.899ρ+7.547	0.565 1	0.740
	φ=−21.92ϕ+29.47	0.550 4	0.753

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。

抗剪强度τ

τ=−61.22ρ²+259.2ρ−213.5

0.946 4

1.562

τ=−555.3ϕ²+351.3ϕ+6.81

0.958 0

1.383

黏聚力c

c=−70.15ρ²+297.1ρ−256.2

0.924 1

2.161

c=−637.2ϕ²+403.7ϕ−3.845

0.929 1

2.089

内摩擦角φ

φ=6.899ρ+7.547

0.565 1

0.740

φ=−21.92ϕ+29.47

0.550 4

0.753

序号	参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。
1	抗剪强度τ	τ=27.4V_p^2.996+1.294V_p^1.498−29.89	0.927 7	1.961
2	黏聚力c	c=33.13V_p^3.022+4.858V_p^1.511−40.64	0.928 4	2.267
3	内摩擦角φ	φ=−4.594×10¹⁰V_p^−75.11+20	0.631 6	0.729

序号

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。

抗剪强度τ

τ=27.4V_p^2.996+1.294V_p^1.498−29.89

0.927 7

1.961

黏聚力c

c=33.13V_p^3.022+4.858V_p^1.511−40.64

0.928 4

2.267

内摩擦角φ

φ=−4.594×10¹⁰V_p^−75.11+20

0.631 6

0.729

序号	参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。
1	抗剪强度τ	τ=27.4V_p^2.996+1.294V_p^1.498−29.89	0.927 7	1.961
2	黏聚力c	c=33.13V_p^3.022+4.858V_p^1.511−40.64	0.928 4	2.267
3	内摩擦角φ	φ=−4.594×10¹⁰V_p^−75.11+20	0.631 6	0.729

序号

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。

抗剪强度τ

τ=27.4V_p^2.996+1.294V_p^1.498−29.89

0.927 7

1.961

黏聚力c

c=33.13V_p^3.022+4.858V_p^1.511−40.64

0.928 4

2.267

内摩擦角φ

φ=−4.594×10¹⁰V_p^−75.11+20

0.631 6

0.729

海底沉积物纵波速度与强度参数关系模型研究

PDF下载

张东昱甫 ¹ , 杨进 ¹^,^* , 王欢欢 ¹ , 李潇 ¹ , 徐飞 ¹

海洋学报 | 论文 2022,44(10): 90-99

收起

海洋学报 | 论文 2022, 44(10): 90-99

海底沉积物纵波速度与强度参数关系模型研究

全屏

张东昱甫¹, 杨进¹^,^*, 王欢欢¹, 李潇¹, 徐飞¹

作者信息

1.中国石油大学（北京）安全与海洋工程学院，北京102200

张东昱甫（1996-），男，北京市人，博士研究生，主要从事海洋地质与海洋地球物理研究。E-mail: zdyf.1013@163.com

通讯作者:

杨进（1966-），男，教授，主要从事石油工程与海洋工程研究。E-mail: cyjin1018@vip.163.com

Study on the relationship model between P-wave velocity and strength parameters of marine sediments

Dongyufu Zhang¹, Jin Yang¹^,^*, Huanhuan Wang¹, Xiao Li¹, Fei Xu¹

Affiliations

1. Safety and Ocean Engineering College, China University of Petroleum (Beijing), Beijing 102200, China

出版时间: 2022-10-01 doi: 10.12284/hyxb2022188

文章导航

摘要

收起

海底沉积物具有质地松软、强度低的特点，其强度参数与海上平台插桩就位、导管架的安装等工程关系紧密，与海上作业安全息息相关。常规钻孔取芯与静力触探等强度参数获取方法成本较高、取样点少，且对沉积物扰动较大，利用易获取的声学资料预测海底沉积物的强度参数具有重要研究意义。本文基于Wood方程、Biot-Stoll模型、Dvorkin等效介质模型等声波传播理论，计算不同物性参数（密度、孔隙度）梯度下的理论纵波速度，结合室内模拟地层声学实验，对比了模型计算声速与实测声速的变化特征，建立了声速与物性参数关系模型。基于室内模拟地层土力学试验，揭示了沉积物物性与抗剪强度、黏聚力等参数之间的关系，建立了物性与强度参数关系模型。以物性参数为桥梁，建立基于声学特征的海底沉积物强度参数预测模型。此模型既避免了原位取样沉积物失水扰动的问题，又弥补了经验公式局限性的缺点，具有普适性与准确性，有效提高无法取样地区沉积物强度参数的精度，提高了经济效益，对浅层的勘探与开发起到了理论指导的作用。

关键词

海底沉积物 / 声波传播理论 / 强度参数 / 物性参数 / 室内模拟实验

Abstract

收起

Marine sediments are characterized by soft texture and low strength, and their strength parameters are closely related to offshore platform pile placement and jacket installation, and are closely related to offshore operation safety. Conventional methods for obtaining strength parameters, such as drilling and cone penetration test, have high cost, few sampling points, and great disturbance to the soil. Therefore, it is of great significance to predict the strength parameters of marine sediments with easily accessible acoustic data. Based on the acoustic propagation theories such as Wood equation, Biot-Stoll model and Dvorkin equivalent medium model, the theoretical P-wave velocity under different physical parameter (density, porosity) gradients was calculated. Combined with the indoor simulated stratum acoustic experiment, the variation characteristics of the calculated and measured sound velocity were compared, and the relationship model between sound velocity and physical parameter was established. The relationship between soil physical parameter and shear strength, cohesive force and other parameters was revealed based on laboratory geotechnical test, and the relationship model between soil physical parameter and strength parameters was established. A prediction model for strength parameters of marine sediments based on acoustic characteristics was established by using physical parameters as bridge. This model not only avoids the problem of water loss disturbance of in-situ sampled soil, but also makes up for the limitation of empirical formula. It has universality and accuracy, and can effectively improve the precision of strength parameters of soil in unsampled areas, improve economic benefits, and play a theoretical guiding role in shallow exploration and development.

Key words

marine sediments / acoustic propagation theory / strength parameters / physical parameters / indoor simulated experiment

引用本文

张东昱甫, 杨进, 王欢欢, 李潇, 徐飞. 海底沉积物纵波速度与强度参数关系模型研究. 海洋学报, 2022 , 44 (10) : 90 -99 . DOI: 10.12284/hyxb2022188

Dongyufu Zhang, Jin Yang, Huanhuan Wang, Xiao Li, Fei Xu. Study on the relationship model between P-wave velocity and strength parameters of marine sediments[J]. Haiyang Xuebao, 2022 , 44 (10) : 90 -99 . DOI: 10.12284/hyxb2022188

正文

收起

1 引言

收起

海底沉积物是海床以下区域沉积地层土的总称，与底部海水相连接，其类型多样且复杂，与区域环境、沉积作用过程等因素相关。随着石油的勘探与开发逐渐转向深海，熟悉与确定海底沉积物特性变得尤为重要。

海底沉积物的强度参数与海上平台插桩就位、导管架的安装等工程关系密切，与海上作业安全息息相关。常规强度参数的获取主要通过海底钻孔取芯或静力触探技术（Cone Penetration Test）来完成。钻孔取芯法发展较早，根据沉积物的种类与取样器类型等因素划分出不同的取样方式^[1-2]。虽然此方式钻入深度范围广，但钻孔数量较少，而且钻孔对于沉积物的扰动较大，导致强度参数精度较差。相较于钻孔取样方法，静力触探法能提供高分辨率的原位测试数据^[3]，具有快速性、可靠性等特点，但其成本较高，应用局限性较大。所以亟需一种准确、便捷、低成本的海底沉积物强度参数获取方法。

在海洋石油钻井中，地球物理勘探是一种非常成熟的研究手段。依据地震资料可以获得海底沉积物的声学信息，如声波速度、声衰减等。其声学特性的分布是海域声场环境的重要组成部分，也是进行海洋多波束或浅地层精细解释的基础资料^[4]。目前，获取海底沉积物声学参数的方法包括声波传播理论与经验公式的构建。其中，声波传播理论将海底沉积物视为固、液两相介质，探究介质之间运动与相互作用对声学特征的影响。Gassmann^[5]考虑到沉积物颗粒之间相互摩擦，建立了饱和沉积物的声学传播公式，但没有考虑沉积物中孔隙流体与骨架的相对运动。Biot^[6–8]考虑到各向异性和黏弹性的问题，建立了流体饱和多孔介质弹性波波动方程。Buckingham^[9-10]提出了一种颗粒间耗散形式，建立了新型声波传播理论。经验公式法是基于海域原位测量结果，建立符合本区域特征的沉积物物性与声学特征的关系。Hamilton^[11]根据北太平洋沉积物的资料，建立了声速与孔隙度、密度、黏土含量等参数的经验公式。侯正瑜等^[12]依据南沙群岛海域沉积物资料，建立声速与孔隙度等参数的经验公式。Zheng等^[13]采用不同频率的声学测量手段，建立了声速与含砂量等参数的经验公式。Kim等^[14]建立了韩国海区声速与平均粒径等参数的经验关系。

考虑到海域地震资料的声学响应特征，探究土体声学与强度特性之间的关系，是在缺乏原位数据情况下，获取沉积物强度参数的有效方式。Duan等^[15]建立了江苏地区土体横波速度与抗剪强度、有效应力间的关系，黄雅虹等^[16]建立了渤海黏土的横波速度与抗剪强度的关系。但原位测量具有地域局限性，建立的相互关系仅仅符合本区域的特点，难以将关系模型进行外部拓展。

本文基于南海沉积物的实际资料，在室内制备不同物性参数（密度、孔隙度）梯度的小尺寸模拟地层，并用声波仪器测取不同模拟地层的纵波速度。依托于声波理论传播模型，带入预设参数求解理论波速，并与实际配置模拟的纵波速度进行对比与验证。采用室内土力学试验测量强度参数，探究海底沉积物声学特性与强度参数的关系，建立基于声学特征的海底沉积物强度参数预测模型，有效提高无法取样地区地层强度参数的精度，为浅层钻井提供理论技术支持。

2 沉积物声学传播理论

收起

因为声波波长远大于沉积物颗粒尺寸大小，因此将海底沉积物视为连续介质。研究初期，多数学者用基质颗粒与孔隙流体两相组分体积模量和密度的加权平均数表示纵波速度，代表理论为Wood方程^[9]。此后Stoll^[17]在Biot的饱和多孔介质声波传播模型基础上，加入了两相组分位移活动的影响。Dvorkin等^[18]和Helgerud等^[19]考虑到骨架和基质颗粒间的关系，对弹性模量进行了进一步修正，建立了等效介质理论。

2.1 Wood理论模型

海底沉积物是一种两相介质，由基质颗粒和孔隙流体组成，在粒间没有摩擦耗散的情况下，纵波速度与两相介质的体积模量和密度相关，计算公式为

(1)

$ {V_{\rm{p}}} = \sqrt {\frac{{{K_0}}}{{{\rho _0}}}} \text{，} $

式中，${V_{\rm{p}}}$为纵波速度；$ {K_0} $为两相介质体积模量；$ {\rho _0} $为两相介质密度。

两相介质的体积模量与密度计算公式为

(2)

$ \frac{1}{{{K_0}}} = \frac{\phi }{{{K_{\rm{f}}}}} + \frac{{1 - \phi }}{{{K_{\rm{r}}}}} \text{，} $

(3)

$ {\rho _0} = {\rho _{\rm{f}}}\phi + {\rho _{\rm{g}}}(1 - \phi ) \text{，} $

式中，${K_{\rm{r}}}$为基质颗粒体积模量；${K_{\rm{f}}}$为孔隙流体体积模量；${\rho _{\rm{f}}}$为孔隙流体密度；${\rho _{\rm{g}}}$为基质颗粒密度；$ \phi $为孔隙度。带入式（1）可求得纵波速度为

(4)

$ {V_{\rm{p}}} = \sqrt {\frac{{{K_{\rm{r}}}{K_{\rm{f}}}}}{{[\phi {\rho _{\rm{f}}} + (1 - \phi ){\rho _{\rm{g}}}][\phi {K_{\rm{r}}} + (1 - \phi ){K_{\rm{f}}}]}}} . $

2.2 Biot-Stoll理论模型

Stoll在Biot模型基础上提出了多孔介质模型，考虑到了介质的孔隙性与弹性。对于海底沉积物，多孔弹性理论模型可以很好地模拟实际情况。Stoll根据可测量的沉积物特性推导出Biot模量的表达式为

(5)

$ D = {K_{\rm{r}}}\left[1 + \phi \left(\frac{{{K_{\rm{r}}}}}{{{K_{\rm{f}}}}} - 1\right)\right] \text{，} $

(6)

$ H = {\frac{{({K_{\rm{r}}} - {K_{\rm{b}}})}}{{D - {K_{\rm{b}}}}}^2} + {K_{\rm{b}}} + \frac{4}{3}{\mu _{\rm{b}}} \text{，} $

(7)

$ C = \frac{{{K_{\rm{r}}}({K_{\rm{r}}} - {K_{\rm{b}}})}}{{D - {K_{\rm{b}}}}} \text{，} $

(8)

$ M = \frac{{{K_{\rm{r}}}^2}}{{D - {K_{\rm{b}}}}} \text{，} $

(9)

$ {K_{\rm{b}}} = {K_{{\rm{b}}0}}\left(1 + \frac{{{\delta _{\rm{b}}}}}{\text{π} }{\rm{i}}\right) \text{，} $

(10)

$ {\mu _{\rm{b}}} = {\mu _{{\rm{b}}0}}\left(1 + \frac{{{\delta _{\rm{s}}}}}{\text{π} }{\rm{i}}\right) \text{，} $

式中，D为中间变量；$H$为弹性模量；$ C $为附加弹性模量；$ M $为复弹性模量；${K_{\rm{b}}}$为骨架体积模量（复数）；i为虚数单位；${\mu _{\rm{b}}}$为骨架剪切模量（复数）； ${\delta _{\rm{b}}}$为体积衰减系数；${\delta _{\rm{s}}}$为剪切衰减系数；骨架体积模量实部$ {K_{{\rm{b}}0}} $，骨架剪切模量实部$ {\mu _{{\rm{b}}0}} $，渗透率$ \kappa $是根据土体本身性质得出^[20]，推导公式为

(11)

$ {K_{{\rm{b}}0}} = \frac{{2{\mu _{{\rm{b}}0}}(1 + \sigma )}}{{3(1 - 2\sigma )}} \text{，} $

(12)

$ {\mu _{{\rm{b}}0}} = 1.835 \times {\text{1}}{{\text{0}}^5}{\left(\frac{\phi }{{1 - \phi }}\right)^{ - 1.12}}\sqrt {(1 - \phi )({\rho _{\rm{g}}} - {\rho _{\rm{f}}})gd} \text{，} $

(13)

$ \kappa = \frac{{{d_{\rm{a}}}^2{\phi ^3}}}{{180{{(1 - \phi )}^2}}}\frac{1}{{\sqrt {10} }}, $

多孔介质中的简谐平面波传播的表达式为

(14)

$ \left| {\begin{array}{*{20}{c}} {H{k^2} - \rho {\omega ^2}}&{{\rho _{\rm{f}}}{\omega ^2} - c{k^2}} \\ {c{k^2} - {\rho _{\rm{f}}}{\omega ^2}}&{m{\omega ^2} - M{k^2} - \dfrac{{\omega F\eta }}{\kappa }{\rm{i}}} \end{array}} \right| = 0 \text{，} $

式中，$ \sigma $为泊松比；$ g $为重力加速度；$ d $为模拟深度；$ {d_{\rm{a}}} $为平均粒径；$ k $为复波数；$\ \rho $为实际密度；$ \omega $为角频率；$ F $为黏性校正系数；$\eta $为绝对黏度；$ m $为考虑了流体流动相对于宏观压力梯度的相位；$ c $为曲折度，其中，

(15)

$ m = \frac{{c{\rho _{\rm{f}}}}}{\phi } . $

求解上述方程，得到以$ {k^2} $为基底二次型方程，利用求根公式解出复波数$ k $并取其实部，带入下式即得纵波速度：

(16)

$ {k_{\rm{r}}} = {{\rm{Re}}} \{ k\} \text{，} $

(17)

$ {V_{\rm{p}}} = \frac{\omega }{{{k_{\rm{r}}}}} . $

2.3 Dvorkin等效介质理论模型

Dvorkin提出等效介质理论，建立了骨架与基质颗粒之间的关系，是饱和多孔介质理论的补充^[18-19]。沉积物骨架的弹性模量与孔隙度、基质颗粒弹性模量、有效应力等因素相关。以临界孔隙度$ {\phi _{\rm{c}}} $为基准分别建立不同情况下干骨架弹性模量的计算公式。

当实际孔隙度大于临界孔隙度时：

(18)

$ {K_{{\rm{Dry}}}} = {\left[\dfrac{{(1 - \phi )/(1 - {\phi _{\rm{c}}})}}{{{K_{\rm{e}}} + \dfrac{4}{3}{\mu _{\rm{e}}}}} + \dfrac{{(\phi - {\phi _{\rm{c}}})/(1 - {\phi _{\rm{c}}})}}{{\dfrac{4}{3}{\mu _{\rm{e}}}}}\right]^{ - 1}} - \dfrac{4}{3}{\mu _{\rm{e}}} \text{，} $

(19)

$ {\mu _{{\rm{Dry}}}} = {\left[\frac{{(1 - \phi )/(1 - {\phi _{\rm{c}}})}}{{{\mu _{\rm{e}}} + Z}} + \frac{{(\phi - {\phi _{\rm{c}}})/(1 - {\phi _{\rm{c}}})}}{Z}\right]^{ - 1}} - Z, $

式中，Z为中间系数，当实际孔隙度小于临界孔隙度时：

(20)

$ {K_{{\rm{Dry}}}} = {\left(\dfrac{{\phi /{\phi _{\rm{c}}}}}{{{K_{\rm{e}}} + \dfrac{4}{3}{\mu _{\rm{e}}}}} + \dfrac{{1 - \phi /{\phi _{\rm{c}}}}}{{{K_{\rm{r}}} + \dfrac{4}{3}{\mu _{\rm{e}}}}}\right)^{ - 1}} - \dfrac{4}{3}{\mu _{\rm{e}}} \text{，} $

(21)

$ {\mu _{{\rm{Dry}}}} = {\left(\dfrac{{\phi /{\phi _{\rm{c}}}}}{{{\mu _{\rm{e}}} + Z}} + \dfrac{{1 - \phi /{\phi _{\rm{c}}}}}{{{\mu _{\rm{r}}} + Z}}\right)^{ - 1}} - Z \text{，} $

(22)

$ Z = \frac{{{\mu _{\rm{e}}}}}{6}\left(\frac{{9{K_{\rm{e}}} + 8{\mu _{\rm{e}}}}}{{{K_{\rm{e}}} + 2{\mu _{\rm{e}}}}}\right), $

有效骨架体积模量$ {K_{\rm{e}}} $，有效骨架剪切模量$ {\mu _{\rm{e}}} $，有效应力$ P $表达式如下：

(23)

$ {K_{\rm{e}}} = {\left[\frac{{{n^2}{{(1 - {\phi _{\rm{c}}})}^2}{\mu _{\rm{r}}}}}{{18{\text{π} ^2}{{(1 - \sigma )}^2}}}P\right]^{\frac{1}{3}}} \text{，} $

(24)

$ {\mu _{\rm{e}}} = \frac{{5 - 4\sigma }}{{5(2 - \sigma )}}{\left[\frac{{3{n^2}{{(1 - {\phi _{\rm{c}}})}^2}{\mu _{\rm{r}}}^2}}{{2{\text{π} ^2}{{(1 - \sigma )}^2}}}P\right]^{\frac{1}{3}}} \text{，} $

(25)

$ P = (\rho - {\rho _{\rm{f}}})gd \text{，} $

式中，${K_{{\rm{Dry}}}}$为干骨架体积模量；${\ \mu _{{\rm{Dry}}}}$为干骨架剪切模量；${\ \mu _{\rm{r}}}$为基质颗粒剪切模量；$ n $为每粒颗粒的平均接触次数。

由于海底沉积物是饱和多孔介质，Gassmann^[5]建立了流体饱和体积模量$ {K_{{\rm{sat}}}} $与孔隙流体体积模量、干骨架体积模量等参数之间的关系。流体饱和体积模量$ {K_{{\rm{sat}}}} $会随着流体的性质而改变，但流体饱和剪切模量$ {\mu _{{\rm{sat}}}} $是不变的，公式为

(26)

$ {K_{{\rm{sat}}}} = {K_{{\rm{Dry}}}} + \dfrac{{{{\left(1 - \dfrac{{{K_{{\rm{Dry}}}}}}{{K_{\rm{r}}}}\right)}^2}}}{{\dfrac{\phi }{{{K_{\rm{f}}}}} + \dfrac{{1 - \phi }}{{K_{\rm{r}}}} - \dfrac{{{K_{{\rm{Dry}}}}}}{{K_{\rm{r}}^2}}}}, $

(27)

$ {\mu _{{\rm{sat}}}} = {\mu _{{\rm{Dry}}}} . $

根据弹性介质理论，纵波速度与流体饱和体积模量、剪切模量和密度相关，公式为

(28)

$ {V_{\rm{p}}} = \sqrt {\dfrac{{{K_{{\rm{sat}}}} + \dfrac{4}{3}{\mu _{{\rm{sat}}}}}}{\rho }} . $

2.4 理论模型参数选取

本文的参数选取中，密度和孔隙度根据实际土样配置情况得到，实验的孔隙流体为水，则流体密度和体积模量参考常温下水的数值。基质颗粒体积模量、基质颗粒剪切模量、泊松比、渗透率、临界孔隙度等其余数值参考经验公式与几位学者预设数值^{[4, 9, 18, 20]}，并结合实际情况与实验特性进行微调（表1）。

3 室内模拟实验

收起

3.1 实验背景

南海蕴藏着许多油藏，是我国油气开采的重点区域，也是国家能源战略的重要地区。南海的地貌大致分为大陆架、大陆坡和深海平原3个部分。南海北部大陆架是华南大陆向南海的延伸部分，现有油气田的开发大多位于本区域内，实验参数范围的选取也参考此处区域土质参数信息。

从横向上看，南海北部大陆架近岸以粉砂质黏土和黏土质粉砂为主，向外延伸至混合沉积区，以粉砂、黏土质砂和砂质黏土为主，至陆架外缘沉积物以细砂分布为主。在纵向上，通过钻孔数据表明，浅部地层上段主要为黏土质粉砂^[21]。对于南沙群岛附近沿海大陆架海域，其粉砂含量范围为38.16%～70.73%，均值为55.37%^[22]，粉砂和黏土是南海沉积物的重要组成部分。在物性参数方面，孔隙度为海底沉积物孔隙体积与总体积的比值，南海北部大陆架沉积物的密度范围为1.48～1.92 g/cm³，均值为1.70 g/cm³^[23]，孔隙度为27.2%～69.7%，均值为53.3%。海南岛东南外海域沉积物密度范围为1.414～1.648 g/cm³，均值为1.565 g/cm³，孔隙度为44.77%～77.35%^[24]。南海南沙群岛海域沉积物密度为1.46～1.83 g/cm³，孔隙度为50.36%～58.56%^[25]。海上原位取样测量因受土体失水等因素影响，造成土质样品受扰动，数据与实际情况偏差较大，并且成本较高。而室内模拟实验可以避免上述问题，按梯度模拟沉积物物性参数范围。通过南海实际沉积物物性参数调查，为实验开展提供数据支持。

本实验采用EDJ-1电动应变直剪仪与微型十字剪切板测量强度参数。声波测量装置采用北京中拓科仪公司ZT802型声波测试仪，声波换能器频率为50 kHz，通过内置计算机与声波测量装置相结合，可进行声波数据的记录与处理。声学与强度参数测量在一个规格为35 cm×20 cm×30 cm的测试箱中进行，箱侧边均标有刻度线，内壁贴有泡沫板作为吸声装置，可防止声波在壁面反射，影响实验效果。将声波换能器分别置于箱底部两侧，将配置好的土样填入制样箱中，进行纵波速度的测量。装置示意图如图1。

3.2 实验步骤

选取100目的石英砂作为砂土部分，选取1 250目的膨润土作为黏土成分，为保持土质的一致性与减少误差，一次性备置足量土样。其中，砂土占总质量的60%，黏土占总质量的40%，通过改变孔隙比与密度的大小，根据公式在测试箱中分别添加不同质量的基础土样和水，保持含水率固定在20%^[26]。为充分模拟海底沉积物特性，加土的同时使用洒水壶均匀浸润土样，并击实土样以排出孔隙中的空气。土水添加完成后使用压实装置对土样充分压实到同一高度（20 cm），使实验土样的体积保持一致。填充完成后立即用保鲜膜紧密盖住土样，将土样静置6 h，使其水分充分均匀渗透进孔隙中。

配置完成后根据密度公式计算土样的密度，配置10组不同密度的土样（表2），相关公式为

(29)

$ {m_{\rm{s}}} = \frac{{{\rho _{\rm{s}}}V}}{{1 + e}}(1 + {\omega _0}) \text{，} $

(30)

$ {m_\text{ω}} = \frac{{{\rho _{\rm{s}}}V}}{{1 + e}}(\omega - {\omega _0}) \text{，} $

式中，$ {m_{\rm{s}}} $为添加土的质量；$ {m_\text{ω} } $为添加水的质量；$ {\rho _{\rm{s}}} $为颗粒密度；$ e $为孔隙比；$ {\omega _0} $为风干土含水率；$ \omega $为设计含水率；$ V $为土样体积。

为了防止土样受剪切发生扰动，先进行声学测量。将声波换能器分别置于箱底部两侧，间距保持10 cm，填土前在换能器上涂抹凡士林使其与土样充分耦合。在箱中对每一组样品进行声波速度的测量，通过测量发射端与接收端信号接收时间，计算出土样的声速。为减小误差，每组土样测量3次取其平均值。其声波测量原理为

(31)

$ {V_{\rm{p}}} = {d_1}/\Delta t \text{，} $

式中，$ {d_1} $为两个声波换能器的间距；$ \Delta t $为发射端与接收端的时间差。

参照土工试验方法标准（GB/T 50123−2019）测量强度参数。首先使用微型十字剪切板测量不同密度土样的抗剪强度，继而使用环刀法对土样进行取样，在不同垂向压力（50 kPa、100 kPa、150 kPa、200 kPa）条件下进行直接剪切试验，采用快剪法并使用较快速度的水平剪切速度（0.8 mm/min）剪切，当剪应力数值达到峰值时，记录土样在不同压力条件下抗剪强度。建立垂向压力与抗剪强度坐标轴，依据库伦公式计算出每组土样的内摩擦角和黏聚力数值。

4 结果分析

收起

4.1 声学特性分析

将3种声学理论传播模型带入预设参数求解出纵波速度，为控制误差大小，密度孔隙度换算关系参考Liu等^[27]对于海底沉积物的原位测量规律，公式为

(32)

$ \rho = - 1.734\phi + 2.748 . $

3种声学传播模型声波速度随着密度的变化趋势相似，密度为1.7 g/cm³区域为临界点，小于此密度时声速保持不变或增加缓慢，超过此密度后声速增加趋势较剧烈。原因在于，密度较小区域，孔隙度较大，孔隙内流体含量较多，密度增加带来的孔隙变化对于土体模量的影响甚微，声速变化不大。随着密度达到一定数值，声波传播主要通过土体的固相框架进行，而密度增加带来孔隙度的变化，会急剧影响土体的模量变化，声速会随着密度的增加急剧升高。相较于其他两种模型，Wood模型求解出的声速较高，密度为1.9 g/cm³时较Biot-Stoll模型高15.9 m/s，较Dvorkin模型高5.8 m/s。虽然Wood模型没有考虑到颗粒间的能量耗散，但与其余两种模型的最大差值仅为1.12%和0.47%，可见物性参数中的密度对于声学传播模型的敏感度较低。

基于室内模拟实验，将室内实测纵波速度与理论模型计算速度进行对比，验证室内实验的准确性（图2）。结果表明，两者曲线趋势相似，当密度小于1.75 g/cm³时，实测值与理论值契合度较高，大于1.75 g/cm³时，实测值略高于理论值。原因在于，当密度大于一定数值时，声速的变化主要依赖于土体的模量变化，但理论模型对于模量变化的估计存在局限性，使其相较于实测数值较小。但考虑到整个密度区间，实测声速与Dvorkin模型求解声速的最大差值为2.91%，与Biot-Stoll模型的最大差值为3.52%，与Wood模型最大差值仅为2.47%。总体来看，室内模拟实验得出的声波速度符合声波理论传播模型。

因本文的3种理论模型较为复杂，参数众多且相互影响，无法根据模型建立物性参数与纵波速度的单一参数模型。基于Gardner等^[28]和Castagna等^[29]提出的地层密度与纵波速度的经验模型，带入实测数值进行公式拟合，得到室内模拟沉积物的纵波速度与物性参数关系（表3），因实测声速契合声学理论传播规律，此关系模型具有准确性。

4.2 强度特性分析

基于室内土力学试验结果，建立不同密度土体直接剪切试验三维图（图3）。图像整体呈斜面趋势，远端高近端低。密度保持一致的条件下，抗剪强度与垂向压力呈正相关。密度为1.8 g/cm³时，垂向压力自50 kPa增加至200 kPa，抗剪强度分级增长16.14 kPa、20.59 kPa、16.68 kPa。同样，当垂向压力保持一致的情况下，随着土质密度的增加，抗剪强度虽有小幅波动，但整体也逐步增加。对于不同孔隙度土体的试验结果（图4），当垂向压力不变时，孔隙度与抗剪强度呈反比的关系。垂向压力为100 kPa时，不同孔隙度土体中抗剪强度的最大差值为17.34 kPa。垂向压力为200 kPa时，最大差值为27.76 kPa。原因在于，模拟土体的密度与孔隙度呈负相关关系，随着土质密度的增加或孔隙度的减小，或是垂向压力的增加，单位体积内的土质变得更加密实，土质颗粒关系愈加紧凑，咬合作用越强，抵抗外部的剪切作用能力越强。

建立不同密度、孔隙度与强度参数的关系（图5，图6）。其中，抗剪强度根据十字剪切板法获取，黏聚力和内摩擦角基于直剪试验结果，使用库伦公式计算得出。整体来看，随着密度的增加或是孔隙度的减小，3种强度参数均呈现增加的趋势，相较于黏聚力与抗剪强度，内摩擦角波动性较大，但基本稳定在一定数值区间。密度为1.8 g/cm³时为黏聚力的极值，孔隙度40.8%时为抗剪强度的极值。

基于室内土力学试验，带入实测数值进行公式拟合，得到沉积物的强度特性与物性参数关系（表4）。其中内摩擦角与物性参数呈线性关系，抗剪强度、黏聚力与物性参数呈二次型关系。除内摩擦角以外，其余两种强度参数模型拟合程度较高。原因在于，随着密度增长，孔隙空间逐渐减小，颗粒之间愈发紧密，形成了公共结合水膜，颗粒之间的胶结作用就更加明显，颗粒接触点的化合键作用增加，使得黏聚力与抗剪强度逐步增加。土体内摩擦角受咬合摩擦与滑动摩擦共同影响，受剪切作用时，密度增加导致土体颗粒移动阻力增大，使得咬合摩擦增大。密度增加导致颗粒间接触点增多，总荷载分配到更多的接触点中，每个接触点上所承受的力就相应减小，使得滑动摩擦增加，但增加的幅度并不一致。综合两种因素可知，随着密度的增长，内摩擦角的涨幅并不均匀，导致拟合程度不高。

4.3 沉积物强度参数预测模型

根据室内实验结果，建立强度参数与纵波速度关系（图5）。整体看强度参数、物性参数（密度）与纵波速度3种变量呈现正相关的关系，即随着其中一种变量的增加，其余两项也随之增长。随着抗剪强度增加，纵波速度呈现初期缓慢增长，中后期快速增加的趋势。在实际海域沉积物的原位测量中，随着地层深度的增加，土质抗剪强度增加，地震资料获取的层位声波时差减小，对应的纵波速度增长，这与实验结果相吻合。黏聚力与声波速度的关系趋势与抗剪强度类似，原因为土质密度增大导致密实程度较高。而内摩擦角的大小波动较为明显，与纵波速度联系不太紧密。

根据前文研究结果，土体物性参数与纵波速度的拟合型式为

(33)

$ \rho = a{V_{\rm{p}}}^b . $

强度参数与物性参数的拟合型式为

(34)

$ \varphi = a'\rho + b'， $

(35)

$ \tau = a'{\rho ^2} + b\rho + c'， $

(36)

$ c = a'{\rho ^2} + b\rho + c'， $

式中，$a、b、a'、b'、 c'$均为待定系数。将式（33）带入式（34）至式（36）中，得到两种型式的强度与声学参数关系模型。其中式（37）适用于内摩擦角与纵波速度的关系式，式（38）适用于抗剪强度、黏聚力与纵波速度的关系式。

(37)

$ y = A{x^B} + C \text{，} $

(38)

$ y = A'{x^{2D}} + B'{x^D} + C' \text{，} $

式中，$ A $、$ B $、$ C $、$ D $、$ A' $、$ B' $、$ C' $均为待定系数。

将实测强度参数与纵波速度带入公式进行拟合，得到基于声学特性的强度参数预测模型（表5）。此模型将物性参数作为中间变量，直接建立强度参数与纵波速度的关系。在海洋石油钻井中，地球物理资料较易获取，所以沉积物的声学资料较为完备，且涉及范围较大。但强度参数往往是通过钻孔取芯或是原位静力触探获取，受生产成本、工作效率等影响，强度参数较为缺乏，且局限于几个点。在现场应用中，基于地球物理方法获取的声学资料，直接代入预测模型可直接求取此区域工程强度性质，节约了成本，提高了经济效益。

因原位测量存在土体失水扰动和应用范围较窄等问题，室内测量可以充分控制土体的物性参数，建立不同梯度参数条件下声学特性与强度参数的联系。

此方法基于土体本身物性参数性质，且通过声波理论传播模型的验证，具有普适性与准确性，对于浅层的勘探与开发起到了理论指导的作用。

此外，相关模型还可用于海底浅层地质灾害预测。如已知区块的声学与强度参数资料符合上述预测模型，则代表此处区域没有浅层地质灾害。反之，若实测资料与模型偏差较大，则表明此处海底沉积物区域可能存在浅层气或浅水流等地质灾害。

5 结论

收起

（1）基于3种声波传播模型，带入预设参数求解纵波速度，比较了不同理论模型之间声速的差异。结果表明，3种模型变化趋势相似，均与土体密度呈正相关的关系，原因在于密度变化引起的土体模量变化，导致声速改变。相较于其他两种理论，Wood模型理论声速偏高，但模型之间声速差值较小，物性参数中的密度对于声波传播模型的敏感度较低。

（2）室内实验测量值与理论模型计算值的对比分析表明，实测声速与理论模型中Biot-Stoll模型的最大误差为3.52%，与Wood模型的最大误差仅为2.47%，室内测量结果符合理论模型趋势，契合度较高。基于Gardner方程与实测数据，建立模拟海底沉积物地层纵波速度与物性参数的关系模型。

（3）室内土力学试验结果表明，随着密度增加或孔隙度减小，抗剪强度和黏聚力的数值逐渐增大，而内摩擦角趋势不太稳定。原因在于，密度增加使得单位体积内的土质更加密实，颗粒关系愈加紧凑，胶结作用明显，导致抗剪强度和黏聚力增加。而内摩擦角受土体咬合摩擦与滑动摩擦共同影响，密度增加导致土体颗粒移动阻力增大，颗粒间接触点增多，总荷载分配到更多的接触点中，每个接触点上所承受的力就相应减小，使得内摩擦角数值增加，但增加的幅度并不一致，导致拟合程度不高。基于实测数据，建立了模拟海底沉积物地层强度特性与物性参数的关系模型。

（4）将物性参数作为中间变量，建立基于声学特性的强度参数预测模型。因实际海域声学资料易获取且较为完备，将实测声速直接带入模型中，可直接求解目标区块的强度参数。此预测模型分梯度控制了土体的物性参数，既避免了原位测量土体失水扰动的问题，又弥补了经验公式地域局限性的缺点，具有普适性与准确性，有效提高了无法取样海域土质强度参数的精度，节约了成本，提高了经济效益，对于浅层的勘探与开发起到了理论指导的作用。

基金

收起

国家自然科学基金(51434009，51774301)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

耿雪樵, 徐行, 刘方兰, 等. 我国海底取样设备的现状与发展趋势[J]. 地质装备, 2009, 10(4): 11−16.

Geng Xueqiao, Xu Xing, Liu Fanglan, et al. The current status and development trends of marine sampling equipment[J]. Equipment for Geotechnical Engineering, 2009, 10(4): 11−16.

王淑云, 鲁晓兵. 深水土工调查技术和分析方法新进展[J]. 海洋工程, 2007, 25(2): 126−130.

Wang Shuyun, Lu Xiaobing. Recent advances in techniques and analytic methods of deepwater geotechnical investigation[J]. The Ocean Engineering, 2007, 25(2): 126−130.

何旭涛, 张秀峰, 舒琪, 等. 海底麻坑内外土体物理力学特性差异研究[J]. 海洋科学, 2020, 44(2): 131−137.

He Xutao, Zhang Xiufeng, Shu Qi, et al. Study on the differences in physical and mechanical properties of soil inside and outside a seabed pockmark[J]. Marine Sciences, 2020, 44(2): 131−137.

王景强, 郭常升, 刘保华, 等. 基于Buckingham模型和Biot-Stoll模型的南沙海域沉积物声速分布特征[J]. 地球学报, 2016, 37(3): 359−367.

Wang Jingqiang, Guo Changsheng, Liu Baohua, et al. Sound speed distribution of seafloor sediments in Nansha Islands sea based on Buckingham model and Biot-Stoll model[J]. Acta Geoscientica Sinica, 2016, 37(3): 359−367.

Gassmann F. Elastic waves through a packing of spheres[J]. Geophysics, 1951, 16(4): 673−685.

Biot M A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range[J]. The Journal of the Acoustical Society of America, 1956, 28(2): 179−191.

Biot M A. Generalized theory of acoustic propagation in porous dissipative media[J]. The Journal of the Acoustical Society of America, 1962, 34(9A): 1254−1264.

Biot M A. Generalized boundary condition for multiple scatter in acoustic reflection[J]. The Journal of the Acoustical Society of America, 1968, 44(6): 1616−1622.

Buckingham M J. Theory of acoustic attenuation, dispersion, and pulse propagation in unconsolidated granular materials including marine sediments[J]. The Journal of the Acoustical Society of America, 1997, 102(5): 2579−2596.

Buckingham M J. Theory of compressional and shear waves in fluidlike marine sediments[J]. The Journal of the Acoustical Society of America, 1998, 103(1): 288−299.

Hamilton E L. Sound velocity and related properties of marine sediments, North Pacific[J]. Journal of Geophysical Research, 1970, 75(23): 4423−4446.

侯正瑜, 郭常升, 王景强. 南沙海域深水区表层沉积物声速与孔隙度相关关系[J]. 海洋科学, 2013, 37(7): 77−82.

Hou Zhengyu, Guo Changsheng, Wang Jingqiang. Surface sediments acoustic velocity and porosity correlation in Nansha sea area abyssal region[J]. Marine Sciences, 2013, 37(7): 77−82.

Zheng Jiewen, Liu Baohua, Kan Guangming, et al. The sound velocity and bulk properties of sediments in the Bohai Sea and the Yellow Sea of China[J]. Acta Oceanologica Sinica, 2016, 35(7): 76−86.

Kim D C, Sung J Y, Park S C, et al. Physical and acoustic properties of shelf sediments, the South Sea of Korea[J]. Marine Geology, 2001, 179(1/2): 39−50.

Duan Wei, Cai Guojun, Liu Songyu, et al. Correlations between shear wave velocity and geotechnical parameters for Jiangsu clays of China[J]. Pure and Applied Geophysics, 2019, 176(2): 669−684.

黄雅虹, 吕悦军, 彭艳菊, 等. 渤海海域黏性土剪切波速与抗剪强度统计关系的初步研究[J]. 中国地震, 2020, 36(3): 527−538.

Huang Yahong, Lü Yuejun, Peng Yanju, et al. A preliminary study on statistical relationship between shear wave velocity and shear strength of clayey soil in Bohai Sea area[J]. Earthquake Research in China, 2020, 36(3): 527−538.

Stoll R D. Acoustic waves in ocean sediments[J]. Geophysics, 1977, 42(4): 715−725.

Dvorkin J, Prasad M, Sakai A, et al. Elasticity of marine sediments: rock physics modeling[J]. Geophysical Research Letters, 1999, 26(12): 1781−1784.

Helgerud M B, Dvorkin J, Nur A, et al. Elastic-wave velocity in marine sediments with gas hydrates: effective medium modeling[J]. Geophysical Research Letters, 1999, 26(13): 2021−2024.

Schock S G. A method for estimating the physical and acoustic properties of the sea bed using chirp sonar data[J]. IEEE Journal of Oceanic Engineering, 2004, 29(4): 1200−1217.

潘国富. 南海北部海底浅部沉积物声学特性研究[D]. 上海: 同济大学, 2003.

Pan Guofu. Research on the acoustic characteristics of seabed sediments in the northern South China Sea[D]. Shanghai: Tongji University, 2003.

Hou Zhengyu, Guo Changsheng, Wang Jingqiang, et al. Seafloor sediment study from South China Sea: acoustic & physical property relationship[J]. Remote Sensing, 2015, 7(9): 11570−11585.

卢博, 李赶先, 黄韶健, 等. 中国黄海、东海和南海北部海底浅层沉积物声学物理性质之比较[J]. 海洋技术, 2005, 24(2): 28−33.

Lu Bo, Li Ganxian, Huang Shaojian, et al. The comparing of seabed sediment acoustic-physical properties in the Yellow Sea, the East China Sea and northern the South China Sea[J]. Ocean Technology, 2005, 24(2): 28−33.

卢博, 李赶先, 刘强, 等. 海南岛东南外海海底沉积物特征及其声学物理性质研究[J]. 海洋学报, 2007, 29(4): 34−42.

Lu Bo, Li Ganxian, Liu Qiang, et al. A study on seafloor sediment and its acouso-physical properties in the southeast offshore sea area of Hainan Island in China[J]. Haiyang Xuebao, 2007, 29(4): 34−42.

Wang Jingqiang, Guo Changsheng, Hou Zhengyu, et al. Distributions and vertical variation patterns of sound speed of surface sediments in South China Sea[J]. Journal of Asian Earth Sciences, 2014, 89: 46−53.

王星华, 黄长溪, 隆威. 直接快剪条件下黏土抗剪强度影响因素探讨[J]. 铁道科学与工程学报, 2012, 9(5): 46−49.

Wang Xinghua, Huang Changxi, Long Wei. Shear strength factors of clayey soils under the condition of direct fast shear[J]. Journal of Railway Science and Engineering, 2012, 9(5): 46−49.

Liu Baohua, Han Tongcheng, Kan Guangming, et al. Correlations between the in situ acoustic properties and geotechnical parameters of sediments in the Yellow Sea, China[J]. Journal of Asian Earth Sciences, 2013, 77: 83−90.

Gardner G H F, Gardner L W, Gregory A R. Formation velocity and density—The diagnostic basics for stratigraphic traps[J]. Geophysics, 1974, 39(6): 770−780.

Castagna J P, Batzle M L, Eastwood R L. Relationships between compressional-wave and shear-wave velocities in clastic silicate rocks[J]. Geophysics, 1985, 50(4): 571−581.

2022年第44卷第10期

PDF下载

195

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022188

接收时间：2022-03-24
首发时间：2026-02-01
出版时间：2022-10-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2022-03-24
修回日期：2022-06-22

基金

国家自然科学基金(51434009，51774301)

作者信息

1.中国石油大学（北京）安全与海洋工程学院，北京102200

通讯作者:

杨进（1966-），男，教授，主要从事石油工程与海洋工程研究。E-mail: cyjin1018@vip.163.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022188

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

序号	参数名称	Wood 模型	Biot-Stoll 模型	Dvorkin等效介质模型
注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。
1	孔隙流体密度ρ_f/（kg·m⁻³）	1 025	1 025	1 025
2	基质颗粒密度ρ_g/（kg·m⁻³）	2 650	2 650	2 650
3	绝对黏度η/（kg·m⁻³）	−	0.001	−
4	密度ρ/（kg·m⁻³）	−	实际	实际
5	孔隙度ϕ/%	实际	实际	实际
6	临界孔隙度ϕ_c/%	−	−	38
7	基质颗粒体积模量K_r/Pa	1.47×10¹⁰	1.47×10¹⁰	1.47×10¹⁰
8	基质颗粒剪切模量µ_r/Pa	−	−	1.34×10¹⁰
9	孔隙流体体积模量K_f/Pa	2.18×10⁹	2.18×10⁹	2.18×10⁹
10	泊松比σ	−	0.15	0.15
11	曲折度c	−	1.35	−
12	渗透率κ/m²	−	公式（13）	−
13	体积衰减系数δ_b	−	0.1	−
14	剪切衰减系数δ_s	−	0.1	−
15	颗粒平均接触次数n	−	−	8.5

序号

参数名称

Wood
模型

Biot-Stoll
模型

Dvorkin等效
介质模型

注：“实际”代表各参数根据实际配置情况得到；“−”代表没有取值。

孔隙流体密度ρ_f/（kg·m⁻³）

1 025

基质颗粒密度ρ_g/（kg·m⁻³）

2 650

绝对黏度η/（kg·m⁻³）

−

0.001

−

密度ρ/（kg·m⁻³）

−

实际

孔隙度ϕ/%

实际

临界孔隙度ϕ_c/%

−

基质颗粒体积模量K_r/Pa

1.47×10¹⁰

基质颗粒剪切模量µ_r/Pa

−

1.34×10¹⁰

孔隙流体体积模量K_f/Pa

2.18×10⁹

泊松比σ

−

0.15

曲折度c

−

1.35

−

渗透率κ/m²

−

公式（13）

−

体积衰减系数δ_b

−

0.1

−

剪切衰减系数δ_s

−

0.1

−

颗粒平均接触次数n

−

8.5

序号	密度/（g·cm⁻³）	孔隙度/%	孔隙比	水质量/g	砂土质量/g	黏土质量/g
1	1.54	51.5	1.06	3 061	11 130	7 420
2	1.57	50.5	1.02	3 122	11 350	7 567
3	1.61	49.2	0.97	3 201	11 638	7 759
4	1.64	48.5	0.94	3 251	11 818	7 879
5	1.67	47.4	0.9	3 320	12 067	8 045
6	1.71	46.2	0.86	3 391	12 327	8 218
7	1.75	44.8	0.81	3 485	12 667	8 445
8	1.8	43.2	0.76	3 583	13 027	8 685
9	1.84	42.5	0.74	3 625	13 177	8 785
10	1.88	40.8	0.69	3 732	13 567	9 045

序号

密度/（g·cm⁻³）

孔隙度/%

孔隙比

水质量/g

砂土质量/g

黏土质量/g

1.54

51.5

1.06

3 061

11 130

7 420

1.57

50.5

1.02

3 122

11 350

7 567

1.61

49.2

0.97

3 201

11 638

7 759

1.64

48.5

0.94

3 251

11 818

7 879

1.67

47.4

0.9

3 320

12 067

8 045

1.71

46.2

0.86

3 391

12 327

8 218

1.75

44.8

0.81

3 485

12 667

8 445

1.8

43.2

0.76

3 583

13 027

8 685

1.84

42.5

0.74

3 625

13 177

8 785

1.88

40.8

0.69

3 732

13 567

9 045

序号	参数	函数表达式	R²	RMSE
注：纵波速度单位为km/s，密度单位为g/cm³。
1	密度ρ	ρ=0.669V_p^2.699	0.943 1	0.029 19
2	孔隙度ϕ	ϕ=1.379V_p^−3.161	0.959 4	0.007 65

序号

参数

函数表达式

R²

RMSE

注：纵波速度单位为km/s，密度单位为g/cm³。

密度ρ

ρ=0.669V_p^2.699

0.943 1

0.029 19

孔隙度ϕ

ϕ=1.379V_p^−3.161

0.959 4

0.007 65

参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。
抗剪强度τ	τ=−61.22ρ²+259.2ρ−213.5	0.946 4	1.562
	τ=−555.3ϕ²+351.3ϕ+6.81	0.958 0	1.383
黏聚力c	c=−70.15ρ²+297.1ρ−256.2	0.924 1	2.161
	c=−637.2ϕ²+403.7ϕ−3.845	0.929 1	2.089
内摩擦角φ	φ=6.899ρ+7.547	0.565 1	0.740
	φ=−21.92ϕ+29.47	0.550 4	0.753

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°）。

抗剪强度τ

τ=−61.22ρ²+259.2ρ−213.5

0.946 4

1.562

τ=−555.3ϕ²+351.3ϕ+6.81

0.958 0

1.383

黏聚力c

c=−70.15ρ²+297.1ρ−256.2

0.924 1

2.161

c=−637.2ϕ²+403.7ϕ−3.845

0.929 1

2.089

内摩擦角φ

φ=6.899ρ+7.547

0.565 1

0.740

φ=−21.92ϕ+29.47

0.550 4

0.753

序号	参数	函数表达式	R²	RMSE
注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。
1	抗剪强度τ	τ=27.4V_p^2.996+1.294V_p^1.498−29.89	0.927 7	1.961
2	黏聚力c	c=33.13V_p^3.022+4.858V_p^1.511−40.64	0.928 4	2.267
3	内摩擦角φ	φ=−4.594×10¹⁰V_p^−75.11+20	0.631 6	0.729

序号

参数

函数表达式

R²

RMSE

注：抗剪强度、黏聚力单位为kPa，内摩擦角单位为（°），纵波速度单位为km/s。

抗剪强度τ

τ=27.4V_p^2.996+1.294V_p^1.498−29.89

0.927 7

1.961

黏聚力c

c=33.13V_p^3.022+4.858V_p^1.511−40.64

0.928 4

2.267

内摩擦角φ

φ=−4.594×10¹⁰V_p^−75.11+20

0.631 6

0.729