CN110188966B - 一种下层土体液化导致表层滑坡的塌块尺寸预测方法 - Google Patents

一种下层土体液化导致表层滑坡的塌块尺寸预测方法 Download PDF

Info

Publication number
CN110188966B
CN110188966B CN201910512623.XA CN201910512623A CN110188966B CN 110188966 B CN110188966 B CN 110188966B CN 201910512623 A CN201910512623 A CN 201910512623A CN 110188966 B CN110188966 B CN 110188966B
Authority
CN
China
Prior art keywords
layer
hard shell
shell layer
liquefaction
lying
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910512623.XA
Other languages
English (en)
Other versions
CN110188966A (zh
Inventor
刘文全
徐兴永
苏乔
付腾飞
陈广泉
李萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
First Institute of Oceanography MNR
Original Assignee
First Institute of Oceanography MNR
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by First Institute of Oceanography MNR filed Critical First Institute of Oceanography MNR
Priority to CN201910512623.XA priority Critical patent/CN110188966B/zh
Publication of CN110188966A publication Critical patent/CN110188966A/zh
Application granted granted Critical
Publication of CN110188966B publication Critical patent/CN110188966B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23Dune restoration or creation; Cliff stabilisation

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Marketing (AREA)
  • Educational Administration (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种下层土体液化导致表层滑坡的塌块尺寸预测方法,包括:获取海床表面的硬壳层的工程性质指标以及厚度;根据硬壳层的工程性质指标以及厚度,计算出硬壳层的抗拉强度以及位于硬壳层下方的下卧液化层顶面处的自重应力;获取硬壳层下方的下卧软弱层的工程性质指标;根据下卧软弱层的工程性质指标,确定出下卧液化层的工程性质指标;根据下卧液化层顶面处的自重应力以及下卧液化层的工程性质指标,确定出下卧液化层与硬壳层之间的界面剪切力;根据硬壳层的抗拉强度、厚度以及界面剪切力,预测出硬壳层在出现滑坡时所将形成的塌块的尺寸。本发明可以针对海床因下层土体液化而导致其表层出现滑坡时所产生的塌块,进行塌块尺寸的有效预测。

Description

一种下层土体液化导致表层滑坡的塌块尺寸预测方法
技术领域
本发明属于海洋工程地质灾害预测技术领域,具体地说,是涉及一种用于预测海底滑坡所产生的塌块的尺寸的计算方法。
背景技术
在我国渤海湾发现的海底滑坡中,曾经观测到在滑坡体表面存在较多厚度不大的塌块,据现场调查和研究表明:这些塌块是由于波浪冲击导致位于海床表面的硬壳层下方的软弱土体出现液化,液化土体沿硬壳层的坡面向下滑动,引发上覆未液化的硬壳层撕裂后产生的。明确塌块的尺寸大小,有助于加深了解海底滑坡的机理。另外,塌块下滑时,会对沿线的海底管道产生冲击力,并且冲击力的大小与塌块的尺寸密切相关。一般来讲,塌块的尺寸越大,其产生的冲击力就越大,继而对海底管道可能造成的危害就越严重。
在目前的海底观测技术领域,虽然已经充分认识到了塌块形成的机理,但是,迄今为止尚未出现有关塌块尺寸的预测方法。
发明内容
本发明的目的在于提供一种塌块尺寸预测方法,可以针对海床因下层土体液化而导致其表层出现滑坡时所产生的塌块,实现塌块尺寸的有效预测。
为解决上述技术问题,本发明采用以下技术方案予以实现:
一种下层土体液化导致表层滑坡的塌块尺寸预测方法,包括:获取海床表面的硬壳层的工程性质指标以及厚度z;根据所述硬壳层的工程性质指标以及厚度z,计算出硬壳层的抗拉强度σt以及位于硬壳层下方的下卧液化层顶面处的自重应力σ;获取位于硬壳层下方的下卧软弱层的工程性质指标;根据所述下卧软弱层的工程性质指标,确定出下卧液化层的工程性质指标;根据所述下卧液化层顶面处的自重应力σ以及下卧液化层的工程性质指标,确定出下卧液化层与硬壳层之间的界面剪切力τ;根据所述硬壳层的抗拉强度σt、厚度z以及所述界面剪切力τ,预测出硬壳层在出现滑坡时所将形成的塌块的尺寸。
其中,所述下卧液化层由所述下卧软弱层液化后形成,位于所述硬壳层与下卧软弱层之间。
优选的,所述硬壳层的工程性质指标包括重度γ1、粘聚力c1和内摩擦角
Figure BDA0002093990150000021
所述硬壳层的抗拉强度
Figure BDA0002093990150000022
所述下卧液化层顶面处的自重应力σ=γ1z。
优选的,所述硬壳层的工程性质指标的获取方式为:在硬壳层处进行钻探,钻透硬壳层取原状样本;针对所述硬壳层的原状样本,采用环刀法确定出硬壳层的密度;将硬壳层的密度乘以重力加速度得到硬壳层的重度γ1;对所述硬壳层的原状样本进行直剪试验,测出硬壳层的粘聚力c1和内摩擦角
Figure BDA0002093990150000023
优选的,所述下卧软弱层的工程性质指标包括粘聚力c2和内摩擦角
Figure BDA0002093990150000024
其获取方式为:在下卧软弱层处进行钻探,取原状样本;对所述下卧软弱层的原状样本进行直剪试验,测出下卧软弱层的粘聚力c2和内摩擦角
Figure BDA0002093990150000025
优选的,所述下卧液化层的工程性质指标包括粘聚力c3和内摩擦角
Figure BDA0002093990150000026
其确定方式优选采用以下两种方式:
其一是,利用公式c3=mcc2
Figure BDA0002093990150000027
计算获得;其中,mc为液化后土体粘聚力折减系数,
Figure BDA0002093990150000028
为液化后土体内摩擦角折减系数;所述mc
Figure BDA0002093990150000029
可以根据经验取值;
其二是,将所述下卧软弱层的原状样本放到动三轴试验机上,经过模拟波浪的循环荷载作用液化后,放入直剪试验仪上测试出下卧液化层的粘聚力c3和内摩擦角
Figure BDA00020939901500000210
优选的,所述下卧液化层与硬壳层之间的界面剪切力τ可以采用以下公式计算获得:
Figure BDA00020939901500000211
优选的,所述塌块的尺寸L的预测方式为
Figure BDA00020939901500000212
其中,所述塌块的尺寸L为塌块在硬壳层滑坡轴线方向上的面积。
与现有技术相比,本发明的优点和积极效果是:本发明针对海床因下层土体液化而导致其表层出现滑坡时所产生的塌块,提出塌块尺寸预测方法。该方法设计简单,容易实现,预测结果误差小,可以协助海洋观测人员准确地预测出海底滑坡对沿线的海底管道可能引发的危害的严重程度,继而帮助海洋观测人员及时、有效地制定应对措施,以降低海底滑坡对海洋工程造成的损失。
结合附图阅读本发明实施方式的详细描述后,本发明的其他特点和优点将变得更加清楚。
附图说明
图1是海床的硬壳层、下卧液化层和下卧软弱层的一种实施例的分层示意图;
图2是本发明所提出的下层土体液化导致表层滑坡的塌块尺寸预测方法的一种实施例的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步详细地说明。
海床在未经历波浪冲击前,其表层通常称之为硬壳层1,硬壳层1下方的土体通常称之为下卧软弱层3,如图1所示。当海床在经历了长时间的波浪冲击后,位于硬壳层1下方的下卧软弱层可能会发生液化,形成下卧液化层2。下卧液化层2沿硬壳层1的坡面向下流动,对上覆硬壳层1会产生向下的剪切力,该剪切力存在于两者之间的界面上。当剪切力超过了硬壳层1的抗拉强度时,硬壳层1就会出现撕裂,进而崩塌成块状,形成塌块。
本实施例基于上述塌块形成机理,提出了一种可以有效预测塌块尺寸的预测方法,如图2所示,具体包括以下过程:
S201、获取海床表面的硬壳层1的工程性质指标;
在本实施例中,所述硬壳层1的工程性质指标具体包括硬壳层的重度γ1、粘聚力c1和内摩擦角
Figure BDA0002093990150000031
为了获取这些工程性质指标,本实施例提出以下优选获取方法:
首先,在硬壳层1处进行钻探,并钻透硬壳层1,以进行原状取样,继而获取硬壳层1的原状样本;
其次,针对获取到的硬壳层1的原状样本,采用环刀法确定出硬壳层1的密度;
然后,将硬壳层1的密度乘以重力加速度,以得到硬壳层1的重度γ1
最后,对硬壳层1的原状样本进行直剪试验,以测量出硬壳层1的粘聚力c1和内摩擦角
Figure BDA0002093990150000044
环刀法和直剪试验均是目前的常规试验方法,因此,本实施例对其具体步骤不做详细说明。
S202、获取硬壳层1的厚度z;
在本实施例中,可以对硬壳层1的原状样本进行厚度测量,以获取到硬壳层1的厚度z。
S203、计算下卧液化层2顶面处的自重应力σ;
在本实施例中,可以根据硬壳层1的重度γ1和厚度z,计算出下卧液化层2顶面处的自重应力σ,其计算公式为:
σ=γ1z。
S204、计算硬壳层1的抗拉强度σt
在本实施例中,可以根据硬壳层1的粘聚力c1和内摩擦角
Figure BDA0002093990150000041
计算出硬壳层1的抗拉强度σt,其计算公式为:
Figure BDA0002093990150000042
S205、获取下卧软弱层3的工程性质指标;
下卧软弱层3为硬壳层1下方未被液化的土体,本实施例主要获取下卧软弱层3的两项工程性质指标,分别为:下卧软弱层3的粘聚力c2和内摩擦角
Figure BDA0002093990150000043
其获取方式优选包括以下过程:
首先,在下卧软弱层3处进行钻探,获取原状样本;
然后,对获取到的下卧软弱层3的原状样本进行直剪试验,测出下卧软弱层3的粘聚力c2和内摩擦角
Figure BDA0002093990150000051
S206、确定下卧液化层2的工程性质指标;
本实施例的下卧液化层2是所述下卧软弱层3液化后形成的流动土体层,位于海床的硬壳层1与所述下卧软弱层3之间,如图1所示。本实施例需要确定出所述下卧液化层2的两项工程性质指标,分别为:下卧液化层2的粘聚力c3和内摩擦角
Figure BDA0002093990150000052
对于这两项工程性质指标的确定,本实施例提出以下两种优选方案:
方案一、采用以下公式计算获得:
c3=mcc2
Figure BDA0002093990150000053
其中,mc为液化后土体粘聚力折减系数,
Figure BDA0002093990150000054
为液化后土体内摩擦角折减系数;所述折减系数mc
Figure BDA0002093990150000055
可以根据经验取值。
方案二、采用动三轴试验和直剪试验方法确定,具体过程为:将下卧软弱层3的原状样本放到动三轴试验机上,经过模拟波浪的循环荷载作用液化后,放入直剪试验仪上测试出下卧液化层2的粘聚力c3和内摩擦角
Figure BDA0002093990150000056
动三轴试验和直剪试验均是目前的常规试验方法,因此,本实施例对其具体试验步骤不做详细说明。
S207、确定下卧液化层2与硬壳层1之间的界面剪切力τ;
在本实施例中,可以根据下卧液化层2顶面处的自重应力σ以及下卧液化层2的粘聚力c3和内摩擦角
Figure BDA0002093990150000057
计算出下卧液化层2与硬壳层1之间的界面剪切力τ。其计算公式为:
Figure BDA0002093990150000058
S208、预测硬壳层1在出现滑坡时所将形成的塌块的尺寸L;
在本实施例中,可以根据硬壳层1的抗拉强度σt、硬壳层1的厚度z以及下卧液化层2与硬壳层1之间的界面剪切力τ,计算出塌块的尺寸L。其计算公式为:
Figure BDA0002093990150000061
在本实施例中,所述塌块的尺寸L是硬壳层1在出现滑坡时,塌块在滑坡轴线方向A上的面积,结合图1所示。
根据预测出的塌块尺寸L的大小,可以预估出塌块所能产生的冲击力的大小,进而对塌块可能造成的危害进行预测,并可以进一步生成海床滑坡预警。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。

Claims (10)

1.一种下层土体液化导致表层滑坡的塌块尺寸预测方法,其特征在于,包括:
获取海床表面的硬壳层的工程性质指标以及厚度z;
根据所述硬壳层的工程性质指标以及厚度z,计算出硬壳层的抗拉强度σt以及位于硬壳层下方的下卧液化层顶面处的自重应力σ;
获取位于硬壳层下方的下卧软弱层的工程性质指标;
根据所述下卧软弱层的工程性质指标,确定出下卧液化层的工程性质指标;
根据所述下卧液化层顶面处的自重应力σ以及下卧液化层的工程性质指标,确定出下卧液化层与硬壳层之间的界面剪切力τ;
根据所述硬壳层的抗拉强度σt、厚度z以及所述界面剪切力τ,预测出硬壳层在出现滑坡时所将形成的塌块的尺寸。
2.根据权利要求1所述的塌块尺寸预测方法,其特征在于,
所述硬壳层的工程性质指标包括重度γ1、粘聚力c1和内摩擦角
Figure FDA0002093990140000011
所述硬壳层的抗拉强度
Figure FDA0002093990140000012
所述下卧液化层顶面处的自重应力σ=γ1z。
3.根据权利要求2所述的塌块尺寸预测方法,其特征在于,所述硬壳层的工程性质指标的获取方式为:
在硬壳层处进行钻探,钻透硬壳层取原状样本;
针对所述硬壳层的原状样本,采用环刀法确定出硬壳层的密度;
将硬壳层的密度乘以重力加速度得到硬壳层的重度γ1
对所述硬壳层的原状样本进行直剪试验,测出硬壳层的粘聚力c1和内摩擦角
Figure FDA0002093990140000013
4.根据权利要求1所述的塌块尺寸预测方法,其特征在于,所述下卧软弱层的工程性质指标包括粘聚力c2和内摩擦角
Figure FDA0002093990140000014
其获取方式为:
在下卧软弱层处进行钻探,取原状样本;
对所述下卧软弱层的原状样本进行直剪试验,测出下卧软弱层的粘聚力c2和内摩擦角
Figure FDA0002093990140000021
5.根据权利要求4所述的塌块尺寸预测方法,其特征在于,所述下卧液化层的工程性质指标包括粘聚力c3和内摩擦角
Figure FDA0002093990140000022
其确定方式为:
c3=mcc2
Figure FDA0002093990140000023
其中,mc为液化后土体粘聚力折减系数,
Figure FDA0002093990140000024
为液化后土体内摩擦角折减系数;所述mc
Figure FDA0002093990140000025
根据经验取值。
6.根据权利要求4所述的塌块尺寸预测方法,其特征在于,所述下卧液化层的工程性质指标包括粘聚力c3和内摩擦角
Figure FDA0002093990140000026
其确定方式为:
将所述下卧软弱层的原状样本放到动三轴试验机上,经过模拟波浪的循环荷载作用液化后,放入直剪试验仪上测试出下卧液化层的粘聚力c3和内摩擦角
Figure FDA0002093990140000029
7.根据权利要求5或6所述的塌块尺寸预测方法,其特征在于,所述下卧液化层与硬壳层之间的界面剪切力τ的确定方式为:
Figure FDA0002093990140000027
8.根据权利要求7所述的塌块尺寸预测方法,其特征在于,所述塌块的尺寸L的预测方式为:
Figure FDA0002093990140000028
9.根据权利要求8所述的塌块尺寸预测方法,其特征在于,所述塌块的尺寸L为塌块在硬壳层滑坡轴线方向上的面积。
10.根据权利要求1至6中任一项所述的塌块尺寸预测方法,其特征在于,所述下卧液化层由所述下卧软弱层液化后形成,位于所述硬壳层与下卧软弱层之间。
CN201910512623.XA 2019-06-13 2019-06-13 一种下层土体液化导致表层滑坡的塌块尺寸预测方法 Active CN110188966B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910512623.XA CN110188966B (zh) 2019-06-13 2019-06-13 一种下层土体液化导致表层滑坡的塌块尺寸预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910512623.XA CN110188966B (zh) 2019-06-13 2019-06-13 一种下层土体液化导致表层滑坡的塌块尺寸预测方法

Publications (2)

Publication Number Publication Date
CN110188966A CN110188966A (zh) 2019-08-30
CN110188966B true CN110188966B (zh) 2022-12-02

Family

ID=67721732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910512623.XA Active CN110188966B (zh) 2019-06-13 2019-06-13 一种下层土体液化导致表层滑坡的塌块尺寸预测方法

Country Status (1)

Country Link
CN (1) CN110188966B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112734081B (zh) * 2020-12-15 2022-06-21 福建工程学院 一种海底滑坡塌块对海底管道冲击力的预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105224704A (zh) * 2014-06-25 2016-01-06 中国科学院海洋研究所 一种海底不稳定性的评价方法
CN108228961A (zh) * 2017-11-27 2018-06-29 华侨大学 一种软弱岩脉海蚀洞洞顶滑塌厚度预测方法
CN108627137A (zh) * 2018-05-09 2018-10-09 中国石油天然气集团公司 一种滑坡变形预测计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101564214B1 (ko) * 2014-08-18 2015-10-29 한국지질자원연구원 간극수압-지표변위-전단면 완전연계 산사태 조기탐지 시험장치

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105224704A (zh) * 2014-06-25 2016-01-06 中国科学院海洋研究所 一种海底不稳定性的评价方法
CN108228961A (zh) * 2017-11-27 2018-06-29 华侨大学 一种软弱岩脉海蚀洞洞顶滑塌厚度预测方法
CN108627137A (zh) * 2018-05-09 2018-10-09 中国石油天然气集团公司 一种滑坡变形预测计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于CDEM对大涵地震滑坡稳定性分析;秦顺海等;《工程技术研究》;20181023(第12期);全文 *
基于拉剪强度同等折减法的滑坡稳定性分析;曾一芳等;《探矿工程(岩土钻掘工程)》;20180210(第02期);全文 *

Also Published As

Publication number Publication date
CN110188966A (zh) 2019-08-30

Similar Documents

Publication Publication Date Title
Liu et al. Laboratory study for pore pressures in sandy deposit under wave loading
CN104164860A (zh) 用于海底浅层土体的自落式孔压动力触探装置
Huang et al. Uplifting behavior of shallow buried pipe in liquefiable soil by dynamic centrifuge test
CN106442937A (zh) 一种新型海洋浅层土体特性探测系统及其评估方法
CN109556653A (zh) 一种带有液压吸力筒基础的原位管土作用测试系统及其测试方法
CN110188966B (zh) 一种下层土体液化导致表层滑坡的塌块尺寸预测方法
CN112734081A (zh) 一种海底滑坡塌块对海底管道冲击力的预测方法
Simpson et al. Internal seiche modes and bottom boundary‐layer dissipation in a temperate lake from acoustic measurements
Kovalev et al. Physical and mechanical characteristics of sea ice in the Kara and Laptev Seas
CN107747306A (zh) 一种用于海洋超软土原位测试的十字形全流触探探头
Bridge Effects of seabed interaction on steel catenary risers
Wu et al. Cyclic resistance evaluation of marine clay based on CPTu data: a case study of Shaba Wind Farm
Burmester et al. Comparing different approaches for calculating wave impacts on a monopile turbine foundation
Abelev et al. Numerical investigation of dynamic free-fall penetrometers in soft cohesive marine sediments using a finite difference approach
Sinha et al. On borehole indentor (BHI) measurements and analysis
CN206177937U (zh) 一种海洋浅层土体特性探测系统
Ge et al. Recent improvements in subsea wellhead fatigue monitoring algorithm and accuracy using verification and calibration techniques
Abyaneh et al. Undrained capacity of suction piles subjected to moment loading
Hossain et al. Spudcan foundations on multi-layered soils with interbedded sand and stiff clay layers
YANG et al. First In-situ Full Deep Sea Sediment Strength Testing Device: MEGE.
Rossiter et al. Assessment of ice-soil interactions: Continuum modeling in clays
Smirnov et al. Large-scale ice crushing experiments with icebreaker
Mumtaz Investigation of pore pressures during high-velocity impact by a free fall penetrometer
CN112432882B (zh) 一种近海水下高压稀释淤泥孔稳定角预测方法
CN115618593A (zh) 一种波浪作用下硬壳层塌块下陷深度预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant