海洋学报

分类	海床以下深度/m	含水量/%	干密度/（g·cm⁻³）	比重	孔隙比	塑限/%	液限/%	塑性指数
未液化粉土	0.5	23.3	1.62	2.71	0.70	18.7	26.8	8.1
	1.0	23.4	1.62	2.70	0.69	19.1	27.0	7.9
	1.5	23.0	1.63	2.71	0.69	17.6	26.8	9.2
	2.5	26.7	1.58	2.70	0.71	21.3	28.5	7.2
	3.0	25.9	1.59	2.71	0.71	21.9	30.5	8.6
	3.5	25.3	1.60	2.70	0.69	21.7	29.9	8.2
	4.3	27.9	1.57	2.70	0.72	22.6	31.6	9.0
	4.7	25.6	1.59	2.70	0.69	21.6	29.7	8.1
已液化粉土	0.5	24.5	1.61	2.70	0.70	20.6	27.6	7.0
	1.0	24.2	1.61	2.70	0.70	20.6	26.3	5.7
	1.5	24.9	1.61	2.70	0.69	18.6	26.1	7.5
	2.5	23.4	1.62	2.70	0.66	20.6	28.1	7.5
	3.0	23.9	1.62	2.70	0.67	20.6	27.9	7.3
	3.5	23.6	1.62	2.71	0.67	19.7	27.7	8.0
	4.3	21.9	1.64	2.70	0.64	15.9	25.2	9.3
	4.7	22.1	1.64	2.71	0.65	16.7	24.5	7.8

分类	海床以下深度/m	含水量/%	干密度/（g·cm⁻³）	比重	孔隙比	塑限/%	液限/%	塑性指数
未液化粉土	0.5	23.3	1.62	2.71	0.70	18.7	26.8	8.1
	1.0	23.4	1.62	2.70	0.69	19.1	27.0	7.9
	1.5	23.0	1.63	2.71	0.69	17.6	26.8	9.2
	2.5	26.7	1.58	2.70	0.71	21.3	28.5	7.2
	3.0	25.9	1.59	2.71	0.71	21.9	30.5	8.6
	3.5	25.3	1.60	2.70	0.69	21.7	29.9	8.2
	4.3	27.9	1.57	2.70	0.72	22.6	31.6	9.0
	4.7	25.6	1.59	2.70	0.69	21.6	29.7	8.1
已液化粉土	0.5	24.5	1.61	2.70	0.70	20.6	27.6	7.0
	1.0	24.2	1.61	2.70	0.70	20.6	26.3	5.7
	1.5	24.9	1.61	2.70	0.69	18.6	26.1	7.5
	2.5	23.4	1.62	2.70	0.66	20.6	28.1	7.5
	3.0	23.9	1.62	2.70	0.67	20.6	27.9	7.3
	3.5	23.6	1.62	2.71	0.67	19.7	27.7	8.0
	4.3	21.9	1.64	2.70	0.64	15.9	25.2	9.3
	4.7	22.1	1.64	2.71	0.65	16.7	24.5	7.8

分类	钻孔编号	海床以下深度/m	固结围压/kPa	循环应力比	初始模拟波浪状态
未液化粉土	ZK1	0.75～1.60	20	0.35, 0.40, 0.50	6.0	60.0	10.0
2.55～3.60	30	0.23, 0.33, 0.37, 0.42, 0.47	4.7	60.0	10.0
3.75～4.70	40	0.23, 0.25, 0.30, 0.44, 0.48	5.2	60.0	10.0
已液化粉土	ZK3	0.95～1.85	20	0.37, 0.45, 0.50	6.4	60.0	10.0
2.48～3.50	30	0.37, 0.42, 0.50	7.5	60.0	10.0
3.70～4.70	40	0.43, 0.45, 0.50	9.0	60.0	10.0

分类	钻孔编号	海床以下深度/m	固结围压/kPa	循环应力比	初始模拟波浪状态
未液化粉土	ZK1	0.75～1.60	20	0.35, 0.40, 0.50	6.0	60.0	10.0
2.55～3.60	30	0.23, 0.33, 0.37, 0.42, 0.47	4.7	60.0	10.0
3.75～4.70	40	0.23, 0.25, 0.30, 0.44, 0.48	5.2	60.0	10.0
已液化粉土	ZK3	0.95～1.85	20	0.37, 0.45, 0.50	6.4	60.0	10.0
2.48～3.50	30	0.37, 0.42, 0.50	7.5	60.0	10.0
3.70～4.70	40	0.43, 0.45, 0.50	9.0	60.0	10.0

分类	围压/kPa	模型	CSR	参数a	参数b	相关系数R²
未液化粉土	20	双曲线	0.35	0.512	0.918	0.99
0.40	0.818	0.761	0.99
0.50	0.375	1.494	0.98
30	双曲线	0.23	0.041	1.025	0.99
0.33	0.042	0.971	0.99
0.37	0.297	0.890	0.99
0.42	0.951	0.845	0.99
0.47	0.662	0.567	0.99
40	双曲线	0.23	0.133	0.973	0.99
0.25	0.154	0.955	0.99
0.30	0.394	1.149	0.99
0.44	0.275	0.890	0.99
0.48	0.365	1.090	0.99
已液化粉土	20	双曲线	0.37	0.146	0.906	0.99
0.45	0.128	1.072	0.99
0.50	0.022	0.969	0.99
30	双曲线	0.37	0.410	1.072	0.98
0.42	0.277	1.148	0.99
0.50	0.519	0.712	0.99
40	双曲线	0.43	0.577	0.427	0.99
0.45	0.292	0.913	0.99
0.50	0.410	0.952	0.98

分类	围压/kPa	模型	CSR	参数a	参数b	相关系数R²
未液化粉土	20	双曲线	0.35	0.512	0.918	0.99
0.40	0.818	0.761	0.99
0.50	0.375	1.494	0.98
30	双曲线	0.23	0.041	1.025	0.99
0.33	0.042	0.971	0.99
0.37	0.297	0.890	0.99
0.42	0.951	0.845	0.99
0.47	0.662	0.567	0.99
40	双曲线	0.23	0.133	0.973	0.99
0.25	0.154	0.955	0.99
0.30	0.394	1.149	0.99
0.44	0.275	0.890	0.99
0.48	0.365	1.090	0.99
已液化粉土	20	双曲线	0.37	0.146	0.906	0.99
0.45	0.128	1.072	0.99
0.50	0.022	0.969	0.99
30	双曲线	0.37	0.410	1.072	0.98
0.42	0.277	1.148	0.99
0.50	0.519	0.712	0.99
40	双曲线	0.43	0.577	0.427	0.99
0.45	0.292	0.913	0.99
0.50	0.410	0.952	0.98

波浪作用下黄河三角洲海底粉土液化特征对比研究

PDF下载

宋玉鹏 ¹^,³ , 孙永福 ²^,³^,^* , 宋丙辉 ¹ , 董立峰 ¹ , 杜星 ¹

海洋学报 | 论文 2021,43(6): 129-138

收起

海洋学报 | 论文 2021, 43(6): 129-138

波浪作用下黄河三角洲海底粉土液化特征对比研究

全屏

宋玉鹏¹^,³, 孙永福²^,³^,^*, 宋丙辉¹, 董立峰¹, 杜星¹

作者信息

¹自然资源部第一海洋研究所海洋工程环境研究中心，山东青岛 266061

²国家深海基地管理中心，山东青岛 266237

³青岛海洋科学与技术试点国家实验室海洋地质过程与环境功能实验室，山东青岛 266237

宋玉鹏（1979-），男，黑龙江省鸡西市人，研究员，博士，主要从事海洋工程地质方面工作。.E-mail：songyupeng@fio.org.cn

通讯作者:

孙永福（1964-），男，山东省潍坊市人，研究员，博士，主要从事海洋工程地质方面研究。E-mail：sunyongfu@fio.org.cn

Comparative study on the liquefaction properties of seabed silt under wave loading in the Huanghe River Delta

Yupeng Song¹^,³, Yongfu Sun²^,³^,^*, Binghui Song¹, Lifeng Dong¹, Xing Du¹

Affiliations

¹Marine Engineering Environment Research Center, First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China

²National Deep Sea Center, Qingdao 266237, China

³Laboratory for Marine Geology and Environment, Pilot National Oceanography Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

出版时间: 2021-06-25 doi: 10.12284/hyxb2021100

文章导航

摘要

收起

作为一种常见的近海海底灾害地质现象，波致海床液化严重威胁着黄河三角洲地区海底工程设施的安全。粉质海床液化后，海底粉土的结构、物理和力学性质均发生了改变，研究该变化规律尤其是评估液化后海底粉土再次发生液化的可能性具有重要的理论意义和应用价值。本文利用室内动三轴仪对取自黄河三角洲已液化和未液化海底粉土开展了液化试验对比研究，讨论了已液化和未液化海底粉土在孔压增长模式和轴向动应变发展趋势方面的异同，对比分析了二者的液化势。研究结果表明：应变标准比孔压标准更适用于评估黄河三角洲地区海底粉土的液化势；孔压和动应变发展模式均表明与未液化粉土相比，已液化海底粉土再次发生液化的抗力有所提高；已液化和未液化海底粉土归一化孔压比u_d/σ₃与循环加载次数比N/N_f间相关关系可采用双曲线或指数函数模型进行定量化描述；未液化海底粉土的波致液化临界循环应力比约为0.20，已液化海底粉土的临界循环应力比约为0.35。研究成果有助于加深对海底粉土波致液化特性的认识，亦可为循环应力历史影响下的土体力学性质研究提供参考。

关键词

黄河三角洲 / 海底粉土 / 再次液化 / 动三轴试验

Abstract

收起

As a common submarine geological disaster, wave-induced seabed liquefaction seriously threatens the safety of subsea engineering facilities in the Huanghe River Delta. The structure, physical and mechanical properties of seabed soil after wave-induced liquefaction all have changed, so it has important theoretical significance and practical value to study on the evaluation of potential possibility of re-liquefaction of seabed soil after previous liquefaction. In this paper, a series of cyclic triaxial liquefaction tests were conducted on core samples collected from submarine non-liquefied and liquefied zone in the Huanghe River Delta, respectively. The differences between non-liquefied and liquefied seabed soil in the developing trends of pore pressure and axial dynamic strain with cycles were analyzed and discussed, and the corresponding liquefaction potentials were also comparatively evaluated. The test results show that compared to pore pressure, the strain standard is more suitable to evaluate the liquefaction potential of the seabed silt in the Huanghe River Delta. The pore pressure and dynamic axial strain development characteristics indicate that the re-liquefaction resistance of the liquefied seabed silt is improved to some extent compared with the non-liquefied silt. Furthermore, the correlations between the normalized pore pressure ratio u_d/σ₃ and the normalized cycle ratio N/N_f could be described quantitatively by the hyperbolic or exponential functions for liquefied and non-liquefied seabed silts. Finally, the critical cyclic stress ratio for the non-liquefied seabed silt is around 0.20 compared to 0.35 for the liquefied one in the Huanghe River Delta. The research findings will contribute to deepening the understanding of the wave-induced liquefaction mechanism of seabed silt, and also provide an example reference for the study of the mechanical properties of soil subjected to previous cyclic stress history.

Key words

Huanghe River Delta / seabed silt / re-liquefaction / cyclic triaxial test

引用本文

宋玉鹏, 孙永福, 宋丙辉, 董立峰, 杜星. 波浪作用下黄河三角洲海底粉土液化特征对比研究. 海洋学报, 2021 , 43 (6) : 129 -138 . DOI: 10.12284/hyxb2021100

Yupeng Song, Yongfu Sun, Binghui Song, Lifeng Dong, Xing Du. Comparative study on the liquefaction properties of seabed silt under wave loading in the Huanghe River Delta[J]. Haiyang Xuebao, 2021 , 43 (6) : 129 -138 . DOI: 10.12284/hyxb2021100

正文

收起

1 引言

收起

黄河是世界上含沙量最高的河流^[1]，每年巨量的泥沙在三角洲地区快速沉积造就了黄河三角洲海底粉土高含水量、高孔隙比、欠固结和低强度的工程特性^[2−3]，并表现出较高的潜在液化势^[4−5]。同时，黄河三角洲地区又是我国重要的滨海油气资源开发区，伴随着海上油气资源的开发利用，该区域海底粉土的工程特性受到了研究人员的持续关注。

前人通过声学测深、浅地层剖面、侧扫声呐、地质钻探及原位测试等多种手段对黄河三角洲海底沉积地层开展了多期次调查研究，发现海底发育有塌陷凹坑、液化扰动土层、滑塌、粉砂流冲沟等多种地质灾害^[6-9]。其中，液化粉土土层由于其分布面积较广、扰动深度大、对海上油气资源开发活动威胁程度高^[10]而备受关注。孙永福等^[11]、宋玉鹏等^[12]通过静力触探和室内物理力学试验，对比研究了该区域液化前后海底粉土的工程地质特征，发现与未液化粉土相比，已液化海底粉土的密度和抗剪强度有所增大而含水率和孔隙比却有所减小；许国辉等^[13]通过现场工程地质调查也发现黄河口液化扰动土层的工程特性好于周边原始土层；众多学者还通过理论分析或物理模型试验探讨了黄河三角洲海底粉土液化扰动土层的形成机理，认为波浪或者风暴潮引起的海床液化起主导作用^{[4-5,8,14-16]}。

鉴于波致液化对海底沉积物性质的显著改造作用，评估已液化海床再次发生液化的潜在可能性是十分必要的，然而目前关于波致液化前后黄河三角洲海底粉土液化势的对比性研究极少。本文依据声学浅地层剖面探测资料对黄河三角洲海底液化海床进行了定位识别，从已液化和临近未液化海床分别采集原位土样，利用室内动三轴仪开展已液化和未液化海底天然粉土波致液化势对比性研究，重点探讨液化前后海底粉土在动孔压和轴向动应变发展演变模式方面的异同，研究成果有助于加深对海底粉土波致液化特性的理解与认识，亦可为循环应力历史影响下的土体力学性质演变研究提供参考。

2 研究区概况

收起

2.1 研究区位置

研究区位于黄河三角洲埕岛海域（图1）。该区域为黄河走河刁口流路及神仙沟流路期间形成的亚三角洲，海底表层多为粉土沉积，形成时代新、固结时间短，在波浪、风暴潮等动力作用下容易发生液化。

2.2 海底液化土层特征

基于海底原位观测、海上工程物探和钻探、室内试验以及数值分析等多种研究手段^[17]，研究区内发现多处海底液化土层，典型液化扰动粉土层浅地层剖面影像如图2所示，未液化粉土层浅地层剖面影像如图3所示。从图中可以看出，与未液化土层相比，已液化海底土层在浅地层声学记录上主要表现为杂乱反射，土层内部无明显层理，原始地层结构遭到显著破坏。

依据高精度浅地层剖面探测结果对研究区内海底液化土层范围进行了划分^[12]，见图1中红线标示区域。在已液化和未液化海床分别进行海上原位静力液压取芯，严格遵循取芯操作流程以获取高质量原位土样，钻孔站位如图1所示，海上GPS定位精度为亚米级。液化土层厚度约5 m，钻孔贯穿整个液化土层。

2.3 海底粉土的物理性质

依据室内土工试验标准对海底粉土的基本物理性质进行了测试，试验结果如表1所示。

由表1可知，海床以下埋深较浅时（深度小于1.5 m），未液化和已液化海底粉土的物理性质彼此比较接近，但随着埋深增大，二者的物理性质表现出较大的差异，即与未液化海底土相比，已液化海底粉土的含水量和孔隙比有所减小，而干密度有所增大，揭示液化后重固结作用有助于提高黄河三角洲海底粉土的排水固结度，使其变的愈加密实^[16]。此外，液塑限测试结果显示，随着埋深的增大，已液化海底粉土的液限和塑限比未液化粉土有所减小，尤其是对于埋深较深（深度为4.3 m和4.7 m）的情况，而二者相应的土体塑性指数却比较接近，显示液化后海底粉土中黏粒（粒径小于0.005 mm）含量有所降低，但控制土体塑性大小的关键胶粒（粒径小于0.002 mm）含量保持相对稳定^[18]。

图4描绘了已液化和未液化海底粉土颗粒组成特征。从图4中也可以看出，相比起未液化粉土，已液化海底粉土中黏粒（粒径小于0.005 mm）含量有所减少，表现出一定的颗粒粗化现象^[19]。

3 海底粉土的液化特性

收起

3.1 试验设备和方案

采用美国GCTS公司SSH-100型动三轴循环剪切试验系统开展海底粉土波致液化特性试验研究。该系统可施加最大垂直荷载达25 kN，循环加载频率范围为0～20 Hz，围压为0～2 MPa，标准试样尺寸：直径为50 mm，高度为100 mm。

为评估已液化海底粉土再次发生液化的可能性，分别选取已液化和未液化原位土样经仔细切取符合上述尺寸要求的试样若干组以便开展室内动三轴液化试验。试样先经过真空抽气饱和而后施加反压进一步饱和，直至试样饱和度超过0.95。依据土样埋深处实际上覆土压力（20～40 kPa）对试样进行各向同性固结，固结稳定标准为30 min内试样轴向变形不超过0.01 mm。固结结束后，对试样施加应力控制等幅正弦波激振荷载，依据波浪特征加载频率设定为0.2 Hz，从工程应用角度出发，同时参考前人有关研究经验^[20]，以单幅轴向动应变达到5%作为液化试验终止条件，循环动应力加载方案如表2所示，其中依据线性波浪理论^[11]，每组循环加载试验起始循环应力比（Cyclic Stress Ratio, CSR）对应的模拟波浪状态亦列于表2中。

3.2 液化试验结果分析

3.2.1 孔压增长模式

孔压增长模式是分析土体液化产生机理及建立液化破坏本构模型的重要参考^[21]。图5绘制了已液化和未液化海底粉土的超孔压随循环加载次数增大的发展变化曲线。注意，图5b中CSR为0.23时，实际循环加载次数N约为2 100次，N大于1 200后孔压发展基本稳定，为了更清楚地展示循环加载前期孔压的发展演变特征，图5b采用了与图5a和图5c相同的横坐标尺度（N=1 200）。

由图5可知，考虑到天然海底土赋存状态的复杂性，超孔压试验结果难免表现出一定的不规律性，但整体来看，黄河三角洲海底粉土在循环加荷初期，土体内部孔压增长迅速，随着循环加载次数N的不断增大，孔压增长速率逐渐减小，某些情况下孔压会逐渐趋于稳定（如图5a已液化土CSR为0.50时）。相同围压σ′下，除个别情况外CSR越大，海底粉土超孔压表现出快速增长的态势。值得注意的是，除图5a中CSR为0.50（已液化粉土）和图5b中CSR为0.33（未液化粉土）时的土样在5%动应变范围内孔压增长达到了有效围压外，其余情况下海底粉土的超孔压均小于其相应的有效围压，说明以孔压标准来评判黄河三角洲海底粉土的液化势存在一定局限性。因此，从工程角度出发，本文后续分析将在孔压标准的基础上主要以单幅轴向动应变达到5%作为海底粉土液化破坏的判别依据^[22]。

为了更合理地对比已液化和未液化海底粉土的孔压增长模式，绘制相同或相近CSR作用下已液化和未液化海底粉土归一化孔压比（u_d/σ₃）−循环加载次数N关系曲线如图6所示。

由图6可知，总的来说，相比起未液化粉土，相同条件下已液化海底粉土的孔压比随循环加载次数的增大而增长的更加缓慢，虽然个别试样的孔压比随CSR的增大出现了不规律发展（如图6a中CSR为0.50和图6b中CSR为0.42时），说明从孔压标准来看，已液化海底粉土再次发生液化破坏的抗力有所增大。另外，从等效循环加载次数^[23]的角度分析可知，CSR越小，即波浪作用越弱，已液化海底粉土抗再次液化的能力提高的越明显。

为定量预测循环荷载作用下土中超静孔隙水压力发展演变趋势，前人根据室内动力试验结果提出了多种归一化孔压比（u_d/σ₃）和归一化循环加载次数比（N/N_f）之间的经验关系式^[24−25]。

幂函数模型：

(1) $\frac{{{u_d}}}{{{\sigma _3}}} = a{\left(\frac{N}{{{N_f}}}\right)^b},$

双曲线模型：

(2) $\frac{{{u_d}}}{{{\sigma _3}}} = \frac{{N/{N_f}}}{{a + b(N/{N_f})}},$

式中，u_d为超静孔隙水压力；σ₃为初始围压；N_f为液化破坏时的循环加载次数；a和b为模型拟合参数。

针对上述模型在求取极限孔压比（u_d/σ₃）时的不足，曾长女等^[26−27]根据室内饱和粉土动三轴试验结果提出了一种孔压发展指数模型：

(3) $\frac{{{u_d}}}{{{\sigma _3}}} = a(1 - {{\rm{e}}^{ - b\frac{N}{{{N_f}}}}}),$

式中模型参数的代表意义同上。

本文尝试采用上述3种经验模型对已液化和未液化海底粉土的归一化孔压比（u_d/σ₃）和归一化循环加载次数比（N/N_f）间相关关系分别进行了曲线拟合，结果表明，幂函数模型无法很好地描述未液化海底粉土在较小CSR情况下（CSR为0.23或0.33）的孔压发展特征，而双曲线和指数模型均可以用来描述未液化海底粉土的孔压增长规律，但双曲线模型的拟合效果相对更好。对已液化海底粉土而言，幂函数模型同样无法很好地描述其在较小围压情况下（σ₃=20 kPa）的孔压发展特征，而双曲线和指数模型均可以用来描述已液化海底粉土的孔压增长规律，但整体而言，仍是双曲线模型的拟合效果相对更优一些，代表性拟合结果如图7和表3所示。

3.2.2 轴向动应变发展模式

孔压增长规律研究表明，在一定动应变范围内（ε_d<5%），海底粉土的动孔压很难增长达到围压，因此探讨循环荷载作用下海底粉土轴向动应变发展特征更具应用价值。

图8描绘了已液化和未液化海底粉土轴向动应变随循环加载次数N的增大而发展演变的趋势。从图8中可以看出，在初始循环加载阶段，海底粉土轴向动应变ε_d先快速增长，随着循环加载次数N的进一步增大，动应变ε_d增长速率呈逐渐减小趋势，这与前述孔压发展模式比较相似。整体而言，海底粉土的轴向动应变随着CSR的增大也表现出增长越来越快的趋势。

郑刚等^[28]依据循环加载作用下黏土轴向动应变发展特征总结提出了4类动应变发展模式：破坏型、直线型、发展型和渐稳型。在一定的动应变范围内（ε_d<5%），对照该分类标准发现，黄河三角洲海底粉土轴向动应变发展模式可归类为发展型和直线型两类，其中发展型表示加载初始阶段轴向动应变快速增大，之后随循环加载次数增幅逐渐放缓（如图8a中已液化土CSR为0.37时）；直线型代表整个循环加载过程中轴向动应变与循环加载次数N之间呈近似线性相关关系（如图8b中未液化土CSR为0.47时）。另外值得注意的是，随着CSR的增大，海底粉土动应变发展模式逐渐由发展型向直线型过渡。

为了更好地对比已液化和未液化海底粉土的轴向动应变发展模式，绘制了相同或相近CSR下已液化和未液化海底粉土轴向动应变ε_d−循环加载次数N关系曲线，如图9所示。总的来说，相比起未液化土，图9中相同条件下已液化海底粉土的动应变随循环加载次数的增大而增长的更加缓慢，虽然个别试样的动应变随CSR的增大出现了不规律发展（如图9b中CSR为0.42时），表明以应变标准来看，已液化海底粉土再次发生液化破坏的抗力也有所增强^[29]。

3.2.3 循环应力比 CSR

为了更直观地对比已液化和未液化海底粉土抗液化能力的强弱，以ε_d=5%作为液化判别标准，依据动三轴试验结果绘制黄河三角洲海底粉土CSR与循环加载破坏次数N_f之间相关关系如图10所示。

由图10可知，整体而言，在同一应变破坏标准下，随着循环加载破坏次数N_f的增大，黄河三角洲海底粉土液化CSR呈非线性减小趋势，其中循环加载破坏次数N_f较小时，CSR衰减速度较快，随着破坏次数N_f增大，CSR衰减速度逐渐减小。对比已液化和未液化海底粉土CSR与N_f相关关系可知，未液化海底粉土的波致液化临界CSR约为0.20，而已液化海底粉土的临界CSR约为0.35，相同循环加载破坏次数N_f下，已液化海底粉土的CSR整体上大于未液化土样的CSR，直观地表明已液化海底粉土再次发生液化破坏的抗力得到了增强。另外值得注意的是，在探讨的试验围压范围内（20～40 kPa），黄河三角洲海底天然粉土液化CSR与N_f关系受围压的影响不甚明显。

4 结论

收起

以往针对黄河三角洲海底粉土开展的波致液化特性研究对液化扰动土层复杂影响的考虑不够，容易误判海底粉土的液化势，造成不利后果。本文通过室内动三轴试验对比研究了已液化和未液化黄河三角洲海底粉土的波致液化特性，得出了如下研究结论：

（1）已液化和未液化海底粉土归一化孔压比u_d/σ₃与循环加载次数比N/N_f间相关关系可采用双曲线或指数函数模型进行描述。

（2）黄河三角洲海底粉土动应变发展模式可分为发展型和直线型两类，随着CSR增大，动应变发展模式逐渐由发展型向直线型过渡。

（3）无论基于孔压标准还是应变标准，黄河三角洲已液化海底粉土再次发生液化的抗力均有所增大。

（4）黄河三角洲未液化海底粉土的波致液化临界CSR约为0.20，已液化海底粉土的临界CSR约为0.35。

基金

收起

国家自然科学基金青年基金（41806079，42007274）；青岛海洋科学与技术试点国家实验室鳌山科技创新计划项目（2015ASKJ03）

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Milliman J D, Meade R H. World-wide delivery of river sediment to the oceans[J]. The Journal of Geology, 1983, 91(1): 1−21.

冯秀丽, 沈渭铨, 杨荣民, 等. 现代黄河口区沉积环境与沉积物工程性质的关系[J]. 青岛海洋大学学报, 1994(S1): 21−28.

Feng Xiuli, Shen Weiquan, Yang Rongmin, et al. Relation between the geotechnical character of sediment and the sedimentary environment of the modern Huanghe Estuary area[J]. Journal of Ocean University of Qingdao, 1994(S1): 21−28.

冯秀丽, 林霖, 庄振业, 等. 现代黄河水下三角洲全新世以来土层岩土工程参数与沉积环境之间的关系[J]. 海岸工程, 1999, 18(4): 1−7.

Feng Xiuli, Lin Lin, Zhuang Zhenye, et al. The relationship between geotechnical parameters and sedimentary environment of soil layers since Holocene in modern Huanghe subaqueous delta[J]. Coastal Engineering, 1999, 18(4): 1−7.

常方强, 贾永刚. 黄河口粉质土海床液化过程的现场试验研究[J]. 土木工程学报, 2012, 45(1): 121−126.

Chang Fangqiang, Jia Yonggang. In-situ test to study silt liquefaction at the subaqueous delta of Yellow River[J]. China Civil Engineering Journal, 2012, 45(1): 121−126.

刘晓磊, 贾永刚, 郑杰文. 波浪导致黄河口海床沉积物超孔压响应现场试验研究[J]. 岩土力学, 2015, 36(11): 3055−3062.

Liu Xiaolei, Jia Yonggang, Zheng Jiewen. In situ experiment of wave-induced excess pore pressure in the seabed sediment in Yellow River Estuary[J]. Rock and Soil Mechanics, 2015, 36(11): 3055−3062.

Prior D B, Yang Z S, Bornhold B D, et al. Active slope failure, sediment collapse, and silt flows on the modern subaqueous Huanghe (Yellow River) delta[J]. Geo-Marine Letters, 1986, 6(2): 85−95.

冯秀丽, 戚洪帅, 王腾, 等. 黄河三角洲埕岛海域地貌演化及其地质灾害分析[J]. 岩土力学, 2004, 25(S1): 17−20.

Feng Xiuli, Qi Hongshuai, Wang Teng, et al. Geomorphological evolution and geological disasters analysis in Chengdao sea area of the Yellow River Delta[J]. Rock and Soil Mechanics, 2004, 25(S1): 17−20.

许国辉, 孙永福, 于月倩, 等. 黄河水下三角洲浅表土体的风暴液化问题[J]. 海洋地质与第四纪地质, 2011, 31(2): 37−42.

Xu Guohui, Sun Yongfu, Yü Yueqian, et al. Storm-induced liquefaction of the surficial sediments in the Yellow River subaqueous delta[J]. Marine Geology & Quaternary Geology, 2011, 31(2): 37−42.

Wang Hu, Liu Hongjun. Evaluation of storm wave-induced silty seabed instability and geo-hazards: a case study in the Yellow River delta[J]. Applied Ocean Research, 2016, 58: 135−145.

马彬彬. 暴风浪作用下海底粉土液化研究[D]. 青岛: 国家海洋局第一海洋研究所, 2015.

Ma Binbin. Liquefaction of seabed silt under storm waves[D]. Qingdao: The First Institute of Oceanography, State Oceanic Administration, 2015.

孙永福, 董立峰, 宋玉鹏. 黄河水下三角洲粉质土扰动土层特征及成因探析[J]. 岩土力学, 2008, 29(6): 1494−1499.

Sun Yongfu, Dong Lifeng, Song Yupeng. Analysis of characteristics and formation of disturbed soil on subaqueous delta of Yellow River[J]. Rock and Soil Mechanics, 2008, 29(6): 1494−1499.

宋玉鹏, 孙永福, 杜星, 等. 黄河水下三角洲液化与未液化粉土的工程地质性质对比研究[J]. 海洋科学进展, 2019, 37(1): 55−64.

Song Yupeng, Sun Yongfu, Du Xing, et al. Comparative study on the difference of engineering geological characteristics between liquefied silt and non-liquefied silt in the Yellow River Delta[J]. Advances in Marine Science, 2019, 37(1): 55−64.

许国辉, 卫聪聪, 孙永福, 等. 黄河水下三角洲浅表局部扰动地层工程特性与成因[J]. 海洋地质与第四纪地质, 2008, 28(6): 19−25.

Xu Guohui, Wei Congcong, Sun Yongfu, et al. The engineering characteristics of shallow disturbed strata and analysis of their formation on the subaqueous Yellow River Delta[J]. Marine Geology & Quaternary Geology, 2008, 28(6): 19−25.

王虎. 波浪作用下黄河三角洲海床失稳机制与评价方法[D]. 青岛: 中国海洋大学, 2015.

Wang Hu. Mechanism and quantitative evaluation of wave-induced seabed instability in the Yellow River delta[D]. Qingdao: Ocean University of China, 2015.

王欣, 许国辉, 孙永福, 等. 黄河水下三角洲液化海底的重新层化及其试验求证[J]. 海洋地质与第四纪地质, 2013, 33(6): 29−40.

Wang Xin, Xu Guohui, Sun Yongfu, et al. Storm-waves-induced seabed sediment liquefaction and Re-stratification on the Yellow River subaqueous delta[J]. Marine Geology & Quaternary Geology, 2013, 33(6): 29−40.

张丽萍, 贾永刚, 侯伟, 等. 液化过程对海床土性质改造的波浪水槽试验[J]. 海洋地质与第四纪地质, 2013, 33(3): 171−180.

Zhang Liping, Jia Yonggang, Hou Wei, et al. Wave flume experiment on seabed reconstruction by liquefaction[J]. Marine Geology & Quaternary Geology, 2013, 33(3): 171−180.

刘晓磊, 贾永刚, 郑杰文. 现代黄河三角洲沉积物波浪动力响应研究评述[J]. 工程地质学报, 2015, 23(S1): 313−319.

Liu Xiaolei, Jia Yonggang, Zheng Jiewen. A review of sediment dynamic responses to waves in the modern Yellow River delta[J]. Journal of Engineering Geology, 2015, 23(S1): 313−319.

郭莹, 王琦. 落锥法确定粉土液限和塑限的试验研究[J]. 岩土力学, 2009, 30(9): 2569−2574.

Guo Ying, Wang Qi. Experimental research on fall cone test to determine liquid limit and plastic limit of silts[J]. Rock and Soil Mechanics, 2009, 30(9): 2569−2574.

王刚, 许国辉, 刘志钦, 等. 波致粉质土液化过程中物理力学性态变化试验研究[J]. 海洋地质与第四纪地质, 2017, 37(1): 176−183.

Wang Gang, Xu Guohui, Liu Zhiqin, et al. Experimental study on the physical and mechanical characteristics changes of wave-induced liquefied silt[J]. Marine Geology & Quaternary Geology, 2017, 37(1): 176−183.

常方强, 贾永刚. 波浪作用下埕岛海域粉质土海床的累积液化[J]. 华侨大学学报(自然科学版), 2013, 34(4): 434−438.

Chang Fangqiang, Jia Yonggang. Residual liquefaction of silt seabed induced by wave at the Chengdao sea area[J]. Journal of Huaqiao University (Natural Science), 2013, 34(4): 434−438.

唐亮, 凌贤长, 徐鹏举, 等. 土体液化动力分析数值模型[J]. 哈尔滨工业大学学报, 2010, 42(4): 521−524.

Tang Liang, Ling Xianzhang, Xu Pengju, et al. Numerical model for dynamic analysis of soil liquefaction[J]. Journal of Harbin Institute of Technology, 2010, 42(4): 521−524.

李飒, 孙兴松, 要明伦. 混黏土的粉土、粉砂室内试验液化判别标准的研究[J]. 岩土力学, 2006, 27(3): 360−364.

Li Sa, Sun Xingsong, Yao Minglun. Study of liquefaction evaluation used in indoor test of silt, silty sand mixed clay[J]. Rock and Soil Mechanics, 2006, 27(3): 360−364.

Seed H B. Soil liquefaction and cyclic mobility evaluation for level ground during earthquakes[J]. Journal of the Geotechnical Engineering Division, 1979, 105(2): 201−255.

谢定义. 土动力学[M]. 北京: 高等教育出版社, 2011.

Xie Dingyi. Soil Dynamics[M]. Beijing: Higher Education Press, 2011.

于濂洪, 王波. 饱和粉土振动孔隙水压力的试验研究[J]. 大连大学学报, 1999, 20(4): 59−62.

Yu Lianhong, Wang Bo. Study on pore water pressure of saturated and disturbed sandy loam during cyclic loading[J]. Journal of Dalian University, 1999, 20(4): 59−62.

曾长女, 刘汉龙, 丰土根, 等. 饱和粉土孔隙水压力性状试验研究[J]. 岩土力学, 2005, 26(12): 1963−1966.

Zeng Changnü, Liu Hanlong, Feng Tugen, et al. Test study on pore water pressure mode of saturated silt[J]. Rock and Soil Mechanics, 2005, 26(12): 1963−1966.

曾长女, 刘汉龙, 陈育民. 细粒含量对粉土动孔压发展模式影响的试验研究[J]. 岩土力学, 2008, 29(8): 2193−2198.

Zeng Changnü, Liu Hanlong, Chen Yumin. Test study on influence of fine particle content on dynamic pore water pressure development mode of silt[J]. Rock and Soil Mechanics, 2008, 29(8): 2193−2198.

郑刚, 霍海峰, 雷华阳, 等. 振动频率对饱和黏土动力特性的影响[J]. 天津大学学报, 2013, 46(1): 38−43.

Zheng Gang, Huo Haifeng, Lei Huayang, et al. Contrastive study on the dynamic characteristics of saturated clay in different vibration frequencies[J]. Journal of Tianjin University, 2013, 46(1): 38−43.

王树英, 阳军生, Luna R. 前期动载对低塑性粉土静态和动态强度的影响[J]. 岩石力学与工程学报, 2013, 32(2): 363−368.

Wang Shuying, Yang Junsheng, Luna R. Effect of previous dynamic loading on static and dynamic strengths of low-plasticity silt[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(2): 363−368.

2021年第43卷第6期

PDF下载

159

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2021100

接收时间：2020-06-18
首发时间：2026-02-26
出版时间：2021-06-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2020-06-18
修回日期：2020-10-19

基金

国家自然科学基金青年基金（41806079，42007274）；青岛海洋科学与技术试点国家实验室鳌山科技创新计划项目（2015ASKJ03）

作者信息

¹自然资源部第一海洋研究所海洋工程环境研究中心，山东青岛 266061

²国家深海基地管理中心，山东青岛 266237

³青岛海洋科学与技术试点国家实验室海洋地质过程与环境功能实验室，山东青岛 266237

通讯作者:

孙永福（1964-），男，山东省潍坊市人，研究员，博士，主要从事海洋工程地质方面研究。E-mail：sunyongfu@fio.org.cn

参考文献

分享链接

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

分享至

全文二维码

扫描看全文

引用本文

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

分类	海床以下深度/m	含水量/%	干密度/（g·cm⁻³）	比重	孔隙比	塑限/%	液限/%	塑性指数
未液化粉土	0.5	23.3	1.62	2.71	0.70	18.7	26.8	8.1
	1.0	23.4	1.62	2.70	0.69	19.1	27.0	7.9
	1.5	23.0	1.63	2.71	0.69	17.6	26.8	9.2
	2.5	26.7	1.58	2.70	0.71	21.3	28.5	7.2
	3.0	25.9	1.59	2.71	0.71	21.9	30.5	8.6
	3.5	25.3	1.60	2.70	0.69	21.7	29.9	8.2
	4.3	27.9	1.57	2.70	0.72	22.6	31.6	9.0
	4.7	25.6	1.59	2.70	0.69	21.6	29.7	8.1
已液化粉土	0.5	24.5	1.61	2.70	0.70	20.6	27.6	7.0
	1.0	24.2	1.61	2.70	0.70	20.6	26.3	5.7
	1.5	24.9	1.61	2.70	0.69	18.6	26.1	7.5
	2.5	23.4	1.62	2.70	0.66	20.6	28.1	7.5
	3.0	23.9	1.62	2.70	0.67	20.6	27.9	7.3
	3.5	23.6	1.62	2.71	0.67	19.7	27.7	8.0
	4.3	21.9	1.64	2.70	0.64	15.9	25.2	9.3
	4.7	22.1	1.64	2.71	0.65	16.7	24.5	7.8

分类

海床以下深度/m

含水量/%

干密度/（g·cm⁻³）

比重

孔隙比

塑限/%

液限/%

塑性指数

未液化粉土

0.5

23.3

1.62

2.71

0.70

18.7

26.8

8.1

1.0

23.4

1.62

2.70

0.69

19.1

27.0

7.9

1.5

23.0

1.63

2.71

0.69

17.6

26.8

9.2

2.5

26.7

1.58

2.70

0.71

21.3

28.5

7.2

3.0

25.9

1.59

2.71

0.71

21.9

30.5

8.6

3.5

25.3

1.60

2.70

0.69

21.7

29.9

8.2

4.3

27.9

1.57

2.70

0.72

22.6

31.6

9.0

4.7

25.6

1.59

2.70

0.69

21.6

29.7

8.1

已液化粉土

0.5

24.5

1.61

2.70

0.70

20.6

27.6

7.0

1.0

24.2

1.61

2.70

0.70

20.6

26.3

5.7

1.5

24.9

1.61

2.70

0.69

18.6

26.1

7.5

2.5

23.4

1.62

2.70

0.66

20.6

28.1

7.5

3.0

23.9

1.62

2.70

0.67

20.6

27.9

7.3

3.5

23.6

1.62

2.71

0.67

19.7

27.7

8.0

4.3

21.9

1.64

2.70

0.64

15.9

25.2

9.3

4.7

22.1

1.64

2.71

0.65

16.7

24.5

7.8

分类	钻孔编号	海床以下深度/m	固结围压/kPa	循环应力比	初始模拟波浪状态
未液化粉土	ZK1	0.75～1.60	20	0.35, 0.40, 0.50	6.0	60.0	10.0
2.55～3.60	30	0.23, 0.33, 0.37, 0.42, 0.47	4.7	60.0	10.0
3.75～4.70	40	0.23, 0.25, 0.30, 0.44, 0.48	5.2	60.0	10.0
已液化粉土	ZK3	0.95～1.85	20	0.37, 0.45, 0.50	6.4	60.0	10.0
2.48～3.50	30	0.37, 0.42, 0.50	7.5	60.0	10.0
3.70～4.70	40	0.43, 0.45, 0.50	9.0	60.0	10.0

分类

钻孔编号

海床以下深度/m

固结围压/kPa

循环应力比

初始模拟波浪状态

波高/m

波长/m

水深/m

未液化粉土

ZK1

0.75～1.60

0.35, 0.40, 0.50

6.0

60.0

10.0

2.55～3.60

0.23, 0.33, 0.37, 0.42, 0.47

4.7

60.0

10.0

3.75～4.70

0.23, 0.25, 0.30, 0.44, 0.48

5.2

60.0

10.0

已液化粉土

ZK3

0.95～1.85

0.37, 0.45, 0.50

6.4

60.0

10.0

2.48～3.50

0.37, 0.42, 0.50

7.5

60.0

10.0

3.70～4.70

0.43, 0.45, 0.50

9.0

60.0

10.0

分类	围压/kPa	模型	CSR	参数a	参数b	相关系数R²
未液化粉土	20	双曲线	0.35	0.512	0.918	0.99
0.40	0.818	0.761	0.99
0.50	0.375	1.494	0.98
30	双曲线	0.23	0.041	1.025	0.99
0.33	0.042	0.971	0.99
0.37	0.297	0.890	0.99
0.42	0.951	0.845	0.99
0.47	0.662	0.567	0.99
40	双曲线	0.23	0.133	0.973	0.99
0.25	0.154	0.955	0.99
0.30	0.394	1.149	0.99
0.44	0.275	0.890	0.99
0.48	0.365	1.090	0.99
已液化粉土	20	双曲线	0.37	0.146	0.906	0.99
0.45	0.128	1.072	0.99
0.50	0.022	0.969	0.99
30	双曲线	0.37	0.410	1.072	0.98
0.42	0.277	1.148	0.99
0.50	0.519	0.712	0.99
40	双曲线	0.43	0.577	0.427	0.99
0.45	0.292	0.913	0.99
0.50	0.410	0.952	0.98

分类

围压/kPa

模型

CSR

参数a

参数b