CN109190291B - 获取动力触探锤击数修正系数的方法 - Google Patents

获取动力触探锤击数修正系数的方法 Download PDF

Info

Publication number
CN109190291B
CN109190291B CN201811124871.9A CN201811124871A CN109190291B CN 109190291 B CN109190291 B CN 109190291B CN 201811124871 A CN201811124871 A CN 201811124871A CN 109190291 B CN109190291 B CN 109190291B
Authority
CN
China
Prior art keywords
correction coefficient
hammering
hammering number
number correction
rod length
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
CN201811124871.9A
Other languages
English (en)
Other versions
CN109190291A (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.)
PowerChina Chengdu Engineering Co Ltd
Original Assignee
PowerChina Chengdu Engineering Co Ltd
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 PowerChina Chengdu Engineering Co Ltd filed Critical PowerChina Chengdu Engineering Co Ltd
Priority to CN201811124871.9A priority Critical patent/CN109190291B/zh
Publication of CN109190291A publication Critical patent/CN109190291A/zh
Application granted granted Critical
Publication of CN109190291B publication Critical patent/CN109190291B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Abstract

本发明涉及动力触探领域，公开了一种获取动力触探锤击数修正系数的方法，解决杆长与实测锤击数同时或者部分缺失的情况下，运用内插法获取锤击数修正系数造成误差大的问题。包括步骤：a.从锤击数修正系数推荐表上获取锤击数所对应的锤击数修正系数值；b.对锤击数修正系数值进行拟合，得到关于锤击数与锤击数修正系数的拟合曲线及关系式；c.对步骤b中所有关系式进行分析，得到关于锤击数与锤击数修正系数的通用关系式；d.对步骤b中所有关系式的系数进行拟合，得到关于杆长与待定系数的拟合曲线及关系式；e.将步骤d得到的关系式带入步骤c中的通用关系式，得到锤击数修正系数计算公式。本发明适用于超重型圆锥动力触探锤击数修正系数获取。

Description

获取动力触探锤击数修正系数的方法
技术领域
本发明涉及动力触探领域，特别涉及获取动力触探锤击数修正系数的方法。
背景技术
动力触探试验与标准贯入试验是岩土工程勘察中常规的原位测试方法，其试验方法为利用一定质量的落锤，以一定高度的自由落距将标准规格的圆锥形探头(标贯器)打入土层中，根据探头贯入的难易程度判定土层的性质。标准贯入试验主要针对于对粘性土、砂土进行原位测试，重型、超重型动力触探一般针对于大颗粒砂类土和碎石类土进行原位测试。根据测试击数，依据相关规范进行修正计算后，可用于定量判断各类土层的物理力学性质，如:砂类土和碎石土的密度,粘性土的状态，以及砂类土、碎石土、粘性土的承载力，砂类土的地震液化现象等。
动力触探试验成果的整理上，前人有较多的研究，各种类型的圆锥动力触探试验均是以贯入一定深度的锤击数作为触探指标，通过与其他室内试验和原位测试指标建立相关关系来获得地基土的物理力学性质指标，从而评价地基土的性质。这种方法简单、直观、使用方便，因此被国内外广泛采用。
目前国内外对实测触探基数是否修正，或采取何种方法修正历来有不同的观点，缺乏统一认识，一般需要进行修正的内容主要包括探杆长度修正、地下水修正、以及上覆土自重修正等。根据工程手册和现行规范，动力触探应用一般局限于深厚覆盖层浅表部(约20m范围内)，根据《岩土工程勘察规范》(GB50021－2001)附录B圆锥动力触探锤击数修正表B.0.1、表B.0.2，重型、超重型动力触探锤击数修正系数如表1和表2。
表1重型圆锥动力触探锤击数修正系数
Figure BDA0001812114950000011
表2超重型圆锥动力触探锤击数修正系数
Figure BDA0001812114950000021
从表1和2不难看出，该方法取值具有一定的局限性：①杆长L不连续，需要采用内插法才能获得；②锤击数实测值N′65、N′120不连续，同样需要采用内插法获得；③杆长L与实测锤击数N′65、N′120同时缺失的情况，运用内插法获取修正系数造成误差较大。
发明内容
本发明要解决的技术问题是：提供一种获取动力触探锤击数修正系数的方法，解决杆长与实测锤击数同时或者部分缺失的情况下，运用内插法获取锤击数修正系数造成误差大的问题。
为解决上述问题，本发明采用的技术方案是：获取动力触探锤击数修正系数的方法，包括步骤：
步骤1：从动力触探的锤击数修正系数推荐表上，获取各种锤击数所对应的锤击数修正系数值；
步骤2：对步骤1获得的锤击数修正系数值进行拟合，得到每种固定杆长下关于锤击数与锤击数修正系数的拟合曲线及关系式；
步骤3：对步骤2中得到的所有关系式进行分析，得到关于锤击数与锤击数修正系数的通用关系式，该通用关系式的系数待定；
步骤4：对步骤2中得到的所有关系式的系数进行拟合，得到关于杆长与待定系数的拟合曲线及关系式；
步骤5：将步骤4得到的关系式带入步骤3中的通用关系式，从而得到基于杆长和锤击数的锤击数修正系数计算公式。
进一步的，步骤3中，所述通用关系式为：α＝AN′120 B，其中，α为锤击数修正系数，N′120为锤击数实测值，A和B分别为待定系数。
进一步的，步骤4得到的关于杆长与待定系数的关系式为：A＝1.006L-0.099，B＝-0.04lnL+0.0174，其中，L为杆长。
进一步的，步骤5得到的锤击数修正系数计算公式为：
Figure BDA0001812114950000031
其中，L为杆长，N′120为锤击数实测值，α为锤击数修正系数。
进一步的，步骤3和步骤4可以采用MATLAB软件进行非线性拟合。
进一步的，步骤3和步骤4可以基于信赖域法进行拟合。
本发明的有益效果是：本发明利用信赖域法、双平方法以及控制变量法结合的优化算法，基于MATLAB数学优化软件，对现行规范中的超重型动力触探锤击数修正系数进行了非线性拟合研究，从而得到了超重型圆锥动力触探锤击数修正系数与杆长、锤击数的函数关系式，经过实施例验证，拟合结果与原始数据具有较高的吻合度，详见表3，说明项目建立的动力触探锤击数修正系数具有一定的准确性和适用性。
附图说明
图1为本发明的流程图。
图2为L＝2m时锤击数修正系数与锤击数的关系图；
图3为L＝3m时锤击数修正系数与锤击数的关系图；
图4为L＝5m时锤击数修正系数与锤击数的关系图；
图5为L＝7m时锤击数修正系数与锤击数的关系图；
图6为L＝9m时锤击数修正系数与锤击数的关系图；
图7为L＝11m时锤击数修正系数与锤击数的关系图；
图8为L＝13m时锤击数修正系数与锤击数的关系图；
图9为L＝15m时锤击数修正系数与锤击数的关系图；
图10为L＝17m时锤击数修正系数与锤击数的关系图；
图11为L＝19m时锤击数修正系数与锤击数的关系图；
图12为系数A与杆长L的拟合曲线图；
图13为系数B与杆长L的拟合曲线图。
具体实施方式
针对现行规范给出的杆长修正系数是非连续的，实际运用过程中存在很大的局限性的问题，本发明提供一种获取动力触探锤击数修正系数的方法，利用信赖域法、双平方法以及控制变量法结合的优化算法，基于MATLAB数学优化软件，对现行规范中的超重型动力触探锤击数修正系数进行了非线性拟合研究，从而得到了超重型圆锥动力触探锤击数修正系数与杆长、锤击数的函数关系式。如图1所示，本发明的具体步骤如下：
步骤1：从动力触探的锤击数修正系数推荐表上，获取各种锤击数所对应的锤击数修正系数值；
步骤2：对步骤1获得的锤击数修正系数值进行回归函数拟合，得到每种固定杆长下关于锤击数与锤击数修正系数的拟合曲线及关系式；
步骤3：对步骤2中得到的所有关系式进行分析，得到关于锤击数与锤击数修正系数的通用关系式，该通用关系式的系数待定；
步骤4：对步骤2中得到的所有关系式的系数进行回归函数拟合，得到关于杆长与待定系数的拟合曲线及关系式；
步骤3和步骤4可基于信赖域法进行拟合。信赖域方法的研究起始于Powell。他提出了一个求解无约束优化问题的算法，该算法在每次迭代时强制性地要求新的迭代点与当前的迭代点之间的距离不超过某一个控制量。引入控制步长是因为传统的线性搜索方法常常由于步长过大而导致算法失败，特别是当问题是病态时尤为如此。控制步长实质上等价于在以当前迭代点为中心的一个邻域内对一个近似于原问题的简单模型求极值。这种技巧可理解为只在一个邻域内对近似模型信赖，所以此邻域被称为信赖域，利用这一技巧的方法也就被称为信赖域法。信赖域的大小通过迭代逐步调节。一般来说，如果在当前迭代模型较好地逼近原问题，则信赖域可扩大，否则信赖域应缩小。
信赖域方法的关键组成部分是如何得信赖域试探步以及怎样决定试探步是否可以接受，试探步一般是子问题的解。所以，如何求得信赖域试探步，实质上归结于子问题的构造，决定试探步是否可被接受通常是利用某一价值函数，看试探步是否使价值函数下降，对于无约束优化问题，价值函数显然就是目标函数，对于约束优化问题，价值函数常常是一罚函数。
步骤3和步骤4可采用MATLAB软件进行非线性拟合。MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂(矩阵实验室)。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式，代表了当今国际科学计算软件的先进水平。
步骤5：将步骤4得到的关系式带入步骤3中的通用关系式，从而得到基于杆长和锤击数的锤击数修正系数计算公式。利用得到的锤击数修正系数计算公式即可对锤击数实测值进行修正。
实施例
以下以获取超重型圆锥动力触探锤击数修正系数的方法为例，对本发明做进一步说明。
从表2不难看出，表2取值具有一定的局限性：①杆长L不连续，需要采用内插法才能获得；②锤击数实测值N′120不连续，同样需要采用内插法获得；③杆长L与锤击数实测值N′120同时缺失的情况，运用内插法获取修正系数造成误差较大。因此，考虑到以上的因素，实施例通过总结规范里面修正系数的规律，利用统计分析原理对其进行了简化归纳，使其运用简便和可靠。
考虑到该表杆长取值(杆长小于20m)是非连续、运用方面存在较大的局限性等因素，本项目对《岩土工程勘察规范》中超重型圆锥动力触探锤击数修正系数与杆长、锤击数的关系进行分析，并得到相应连续的二元函数。
根据前人研究可知，超重型圆锥动力触探锤击数修正系数主要与杆长、锤击数有关系，通过控制变量法，控制其余变量，找出单一变量条件下单个因子(锤击数)与超重型圆锥动力触探锤击数修正系数的函数关系，为下一步将各因子与超重型圆锥动力触探锤击数修正系数进行多元函数拟合提供依据。
在相同杆长L的条件下，将超重型圆锥动力触探锤击数修正系数与锤击数进行线性回归分析。从而，可以得到不同杆长条件下超重型圆锥动力触探锤击数修正系数与锤击数的关系式。
综合上述分析，实施例提供了一种获取超重型圆锥动力触探锤击数修正系数的方法，其具体步骤如下：
S1：从超重型圆锥动力触探的锤击数修正系数推荐表(表2)上，获取各种锤击数所对应的锤击数修正系数值,即获取L＝1、2、3、5、7、9、11、13、15、17、19时，锤击数N′120为分别1、3、5、7、9、10、15、20、25、30、35、40所对应的锤击数修正系数值。
S2：利用信赖域法，对步骤1获得的锤击数修正系数值进行拟合，得到每种固定杆长下关于锤击数与锤击数修正系数的拟合曲线及关系式。L＝2、3、5、7、9、11、13、15、17、19时拟合的图形分别对应图2-图11所示，图2-图11对应的关系式分别为：α＝0.935N′120 -0.017、α＝0.9034N′120 -0.030、α＝0.8589N′120 -0.048、α＝0.8325N′120 -0.061、α＝0.8110N′120 -0.073、α＝0.7954N′120 -0.084、α＝0.7847N′120 -0.094、α＝0.7653N′120 -0.097、α＝0.7590N′120 -0.107、α＝0.7461N′120 -0.112。
S3：对步骤S2中得到的所有关系式进行分析，可得到关于锤击数与锤击数修正系数的通用关系式：α＝A N′120 B，其中，α为锤击数修正系数，N′120为锤击数实测值，A和B分别为待定系数。
S4：对步骤S2中得到的所有关系式的系数进行拟合，得到关于杆长与待定系数A、B的拟合曲线及关系式。
将步骤S2中得到的所有关系式的系数取出，即有：
A＝(0.935，0.9034，0.8589，0.8325，0.8110，0.7954，0.7847，0.7653，0.7590，0.7461)
B＝(-0.017，-0.030，-0.048，-0.061，-0.073，-0.084，-0.094，-0.097，-0.107，-0.112)
考虑到杆长对超重型圆锥动力触探锤击数修正系数影响较大且系数A、B随杆长变化呈规律性变化，将系数A和B分别与杆长L(m)进行线性回归，得到以下式子以及图12和图13：
A＝1.006L-0.099
B＝-0.04lnL+0.0174
S5：将步骤S4得到的关系式带入步骤S3中的通用关系式，从而得到基于杆长和锤击数的锤击数修正系数计算公式。
将A、B与杆长L的拟合关系式代入α＝A N′120 B中可得1<L≤20m的杆长修正系数：
Figure BDA0001812114950000061
因此，可将《岩土工程勘察规范》中超重型圆锥动力触探锤击数修正系数计算公式总结为连续的二元函数：
Figure BDA0001812114950000062
由于N120＝αN′120，因此，修正后锤击数N120可以表示为：
Figure BDA0001812114950000063
其中，N′120为锤击数实测值。
采用信赖域法，利用Matlab软件进行拟合分析，从实施例拟合结果可以得出实施例的SEE(残差平方和)和RMSE(均方差)均较小，且相关性系数基本上均接近1，说明拟合精度较高。
我们选取《岩土工程勘察规范》中超重型圆锥动力触探锤击数修正系数中杆长L＝3m，5m，7m的数据对上述锤击数修正系数计算公式进行误差分析。
根据上述计算公式得到了不同样品的修正系数计算值α，同时，根据以下公式分别计算得到了α与α0之间的绝对误差和相对误差。
绝对误差Δα可以表示为：
Δα＝α-α0
相对误差δ可以表示为：
δ＝(α-α0)/α
公式中α0为由通过查表2得到的修正系数，根据绝对误差和相对误差公式对拟合结果进行误差分析，可以看出，拟合结果与原始数据具有较高的吻合度，详见表3。说明项目建立的动力触探锤击数修正系数具有一定的准确性和适用性。
表3采用实施例方法计算锤击数修正系数误差统计表
Figure BDA0001812114950000071

Claims (5)

1.获取动力触探锤击数修正系数的方法，其特征在于，包括步骤：
步骤1：从动力触探的锤击数修正系数推荐表上，获取各种锤击数所对应的锤击数修正系数值；
步骤2：对步骤1获得的锤击数修正系数值进行拟合，得到每种固定杆长下关于锤击数与锤击数修正系数的拟合曲线及关系式；
步骤3：对步骤2中得到的所有关系式进行分析，得到关于锤击数与锤击数修正系数的通用关系式，该通用关系式的系数待定；
步骤4：对步骤2中得到的所有关系式的系数进行拟合，得到关于杆长与待定系数的拟合曲线及关系式；
步骤5：将步骤4得到的关系式带入步骤3中的通用关系式，从而得到基于杆长和锤击数的锤击数修正系数计算公式：
Figure FDA0003671658560000011
其中，L为杆长，N′120为锤击数实测值，α为锤击数修正系数。
2.如权利要求1所述的获取动力触探锤击数修正系数的方法，其特征在于，步骤3中，所述通用关系式为：α＝AN′120 B，其中，α为锤击数修正系数，N′120为锤击数实测值，A和B分别为待定系数。
3.如权利要求2所述的获取动力触探锤击数修正系数的方法，其特征在于，步骤4得到的关于杆长与待定系数的关系式为：A＝1.006L-0.099，B＝-0.04lnL+0.0174，其中，L为杆长。
4.如权利要求1所述的获取动力触探锤击数修正系数的方法，其特征在于，步骤3和步骤4采用MATLAB软件进行拟合。
5.如权利要求1所述的获取动力触探锤击数修正系数的方法，其特征在于，步骤3和步骤4基于信赖域法进行拟合。
CN201811124871.9A 2018-09-26 2018-09-26 获取动力触探锤击数修正系数的方法 Active CN109190291B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811124871.9A CN109190291B (zh) 2018-09-26 2018-09-26 获取动力触探锤击数修正系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811124871.9A CN109190291B (zh) 2018-09-26 2018-09-26 获取动力触探锤击数修正系数的方法

Publications (2)

Publication Number Publication Date
CN109190291A CN109190291A (zh) 2019-01-11
CN109190291B true CN109190291B (zh) 2022-07-08

Family

ID=64906168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811124871.9A Active CN109190291B (zh) 2018-09-26 2018-09-26 获取动力触探锤击数修正系数的方法

Country Status (1)

Country Link
CN (1) CN109190291B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112942290B (zh) * 2021-01-29 2023-09-08 北京城建勘测设计研究院有限责任公司 一种实时计算原位测试结果的工程勘察外业采集系统
CN113591305B (zh) * 2021-07-30 2023-04-28 中国电建集团成都勘测设计研究院有限公司 超重型动力触探试验的锤击数修正方法
CN113569414B (zh) * 2021-07-30 2023-05-26 中国电建集团成都勘测设计研究院有限公司 超重型动力触探试验的锤击数修正系数模型的构建方法
CN115046726B (zh) * 2022-08-16 2022-11-01 中国航空工业集团公司沈阳空气动力研究所 一种模型和天平的自重修正方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101922159A (zh) * 2010-07-21 2010-12-22 江苏省交通科学研究院股份有限公司 碎石桩质量快速检测方法
CN102879551A (zh) * 2012-10-24 2013-01-16 中国地震局工程力学研究所 一种砂砾土力学性能评价方法
CN103911980A (zh) * 2013-01-08 2014-07-09 朱文凯 动力触探自动记录仪
CN105044776A (zh) * 2015-07-31 2015-11-11 中国能源建设集团江苏省电力设计院有限公司 基于abaqus的土体地基液化研究方法
CN106198266A (zh) * 2016-06-30 2016-12-07 长江三峡勘测研究院有限公司（武汉） 一种超长杆重型、超重型圆锥动力触探锤击数修正方法
CN106483018A (zh) * 2016-11-08 2017-03-08 中国水利水电科学研究院 考虑原位结构效应确定深厚覆盖层土体的动强度参数的方法
CN107169224A (zh) * 2017-06-08 2017-09-15 天津大学 基于cptu测试进行大直径超长管桩可打入性分析方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140302914A1 (en) * 2013-03-15 2014-10-09 Bram Weinstein Computer-based Methods and Systems for Fantasy Sports Gaming
DE102016008819A1 (de) * 2016-07-19 2018-01-25 Liebherr-Werk Nenzing Gmbh Schlagzahloptimierung

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101922159A (zh) * 2010-07-21 2010-12-22 江苏省交通科学研究院股份有限公司 碎石桩质量快速检测方法
CN102879551A (zh) * 2012-10-24 2013-01-16 中国地震局工程力学研究所 一种砂砾土力学性能评价方法
CN103911980A (zh) * 2013-01-08 2014-07-09 朱文凯 动力触探自动记录仪
CN105044776A (zh) * 2015-07-31 2015-11-11 中国能源建设集团江苏省电力设计院有限公司 基于abaqus的土体地基液化研究方法
CN106198266A (zh) * 2016-06-30 2016-12-07 长江三峡勘测研究院有限公司（武汉） 一种超长杆重型、超重型圆锥动力触探锤击数修正方法
CN106483018A (zh) * 2016-11-08 2017-03-08 中国水利水电科学研究院 考虑原位结构效应确定深厚覆盖层土体的动强度参数的方法
CN107169224A (zh) * 2017-06-08 2017-09-15 天津大学 基于cptu测试进行大直径超长管桩可打入性分析方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Die Energiewende - Blueprints for the new energy age》.2016,1-6. *
Philipp Erlinghagen 等.Comparison of Dynamic Simulations and Steady-State Calculations of Fault Currents in Distribution Grids with high DG Penetration.《International ETG Congress 2015 *
动力触探杆长适应性及其修正试验;李会中 等;《地球科学》;20160715;第41卷(第6期);1249-1258 *
计算动力触探修正系数的多项式拟合法;鲁兴社 等;《长春工程学院学报(自然科学版)》;20110328;第12卷(第1期);39-43 *
超重型圆锥动力触探杆长修正优化与成果转换研究;单诗涵 等;《地下空间与工程学报》;20201115;第16卷(第S2期);885-890,933 *

Also Published As

Publication number Publication date
CN109190291A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109190291B (zh) 获取动力触探锤击数修正系数的方法
Douglas et al. Recent and future developments in earthquake ground motion estimation
Wang et al. Bayesian perspective on geotechnical variability and site characterization
Zalachoris et al. Evaluation of one-dimensional site response techniques using borehole arrays
Aimene et al. Geomechanical modeling of hydraulic fractures interacting with natural fractures—Validation with microseismic and tracer data from the Marcellus and Eagle Ford
CN113901681B (zh) 一种全寿命周期页岩气储层双甜点三维可压性评估方法
Ching et al. Impact of statistical uncertainty on geotechnical reliability estimation
Bahrampouri et al. Mapping the uncertainty in modulus reduction and damping curves onto the uncertainty of site amplification functions
Chartier et al. SHERIFS: Open‐source code for computing earthquake rates in fault systems and constructing hazard models
CN113820750A (zh) 基于弹塑性力学定量预测泥岩构造裂缝的方法
Kim et al. Probabilistic approaches to estimating variation in the mechanical properties of rock masses
CN109034641A (zh) 管道缺陷预测方法及装置
US11126694B2 (en) Automatic calibration for modeling a field
CN109085643A (zh) 早至波的分步联合反演方法
Tran et al. Probabilistic models for uncertainty quantification of soil properties on site response analysis
Wang et al. Probabilistic seismic displacement hazard assessment of earth slopes incorporating spatially random soil parameters
Zetterlund et al. Framework for value of information analysis in rock mass characterization for grouting purposes
CN112231874A (zh) 地下管线模型的建立方法、装置、计算机设备和存储介质
Vardanega et al. Stiffness of clays and silts: Modeling considerations
Han et al. Borehole cross-sectional shape analysis under in situ stress
Amjadi et al. Stochastic nonlinear ground response analysis considering existing boreholes locations by the geostatistical method
Zhao et al. SPT-CPT correlation and its application for liquefaction evaluation in China
Ku et al. Liquefaction and cyclic softening potential of soils–a unified piezocone penetration testing-based approach
Jia et al. A model of undrained stress–strain curves considering stress path and strain softening
WO2023049138A1 (en) Automatic sensor data validation on a drilling rig site

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