CN116296820B - 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 - Google Patents
双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 Download PDFInfo
- Publication number
- CN116296820B CN116296820B CN202310339341.0A CN202310339341A CN116296820B CN 116296820 B CN116296820 B CN 116296820B CN 202310339341 A CN202310339341 A CN 202310339341A CN 116296820 B CN116296820 B CN 116296820B
- Authority
- CN
- China
- Prior art keywords
- stress
- crack
- sigma
- stress sigma
- strain
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000005336 cracking Methods 0.000 title claims abstract description 16
- 230000002457 bidirectional effect Effects 0.000 title claims abstract description 13
- 230000000977 initiatory effect Effects 0.000 claims abstract description 118
- 238000004364 calculation method Methods 0.000 claims description 49
- 238000004590 computer program Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 14
- 230000006870 function Effects 0.000 claims description 11
- 238000012669 compression test Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 9
- 230000009471 action Effects 0.000 description 3
- 230000005489 elastic deformation Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002927 high level radioactive waste Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0064—Initiation of crack
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品，涉及岩石力学技术领域。本发明实施例中，通过判断裂纹体积应变法得出的裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)是否满足容许误差，并对假设的岩石弹性段的下限应力和上限应力进行迭代更新，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，从而从下限应力和上限应力两个方向逼近试样的弹性段，避免了人为主观因素的影响，提高了弹性参数、裂纹闭合应力和起裂应力的精确度。此外，本发明实施例通过双向迭代可以同时获取试样的弹性参数，裂纹闭合应力和起裂应力，提升工作效率。
Description
技术领域
本发明实施例涉及岩石力学技术领域，尤其涉及一种双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品。
背景技术
随着我国地下工程的大规模发展，如石油天然气开采、地下隧洞工程和高放废物深埋处置等，岩石力学性能的研究越来越受到重视。弹性模量和泊松比作为研究岩体变形性质最基本的参数，准确确定其取值具有重要意义。岩石的渐进破坏过程可以划分为五个阶段：裂纹闭合-弹性变形-裂纹起裂和稳定扩展-非稳定扩展-峰后变形阶段。其中，当岩石所受轴向应力达到裂纹闭合应力(σcc)时，岩石进入弹性变形阶段；继续加载到裂纹起裂应力(σci)，微裂纹将开始以稳定的速度扩展，被测试样的弹性变形结束。因此，识别裂纹闭合应力和起裂应力是准确获取应力-应变曲线的弹性段从而确定岩石弹性参数的关键。此外，当工程围岩在应力达到其裂纹起裂强度时，岩体强度及性能通过产生裂纹、断裂等形式退化，而岩体强度脆性劣化往往会给地下工程带来突水、岩爆等问题，因此准确识别岩石起裂强度对围岩稳定性评价和支护设计等方面具有重要意义。
可见，目前亟需一种可以确定岩石弹性参数与起裂应力的技术方案。
发明内容
本发明实施例提供一种双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品，以至少部分解决相关技术中存在的问题。
本发明实施例第一方面提供了一种双向逼近确定岩石弹性参数与起裂应力的迭代方法，所述方法包括：
确定下限应力σa(k)和上限应力σb(k)；
根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk；
根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)；
在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新；
根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力。
可选地，在首次计算时，将初始下限应力σa(1)确定为所述应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)确定为应力-应变曲线中峰值应力的一半。
可选地，根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，包括：
根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和下式确定裂纹体积应变随轴向应变变化的曲线：
其中，εV,cr为裂纹体积应变，εz为轴向应变，ε1为横向应变，σz为施加的轴向应力。
可选地，所述方法还包括：
通过下式计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
通过下式计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差。
可选地，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，包括：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；
在χcc(k)不大于0的情况下，采用所述下限应力σa(k)作为下一次迭代计算的下限应力。
可选地，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新，包括：
在χci(k)大于0的情况下，采用所述上限应力σb(k)作为下一次迭代计算的上限应力；
在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
可选地，所述方法还包括：
采集所述岩石的试样的常温单轴压缩试验数据，得到所述应力-应变曲线。
基于同一发明构思，本发明实施例第二方面提供了一种双向逼近确定岩石弹性参数与起裂应力的迭代装置，所述装置包括：
第一确定模块，用于确定下限应力σa(k)和上限应力σb(k)；
第二确定模块，用于根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk；
第三确定模块，用于根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)；
更新模块，用于在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新；
第四确定模块，用于根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力。
可选地，在首次计算时，将初始下限应力σa(1)确定为所述应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)确定为应力-应变曲线中峰值应力的一半。
可选地，所述第二确定模块，具有用于：
根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和下式确定裂纹体积应变随轴向应变变化的曲线：
其中，εV,cr为裂纹体积应变，εz为轴向应变，ε1为横向应变，σz为施加的轴向应力。
可选地，所述装置还包括：
第一计算模块，用于通过下式计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
第二计算模块，用于通过下式计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
判断模块，用于基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差。
可选地，所述更新模块，具体用于：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；
在χcc(k)不大于0的情况下，采用所述下限应力σa(k)作为下一次迭代计算的下限应力。
可选地，所述更新模块，具体用于：
在χci(k)大于0的情况下，采用所述上限应力σb(k)作为下一次迭代计算的上限应力；
在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
可选地，所述装置还包括：
采集模块，用于采集所述岩石的试样的常温单轴压缩试验数据，得到所述应力-应变曲线。
基于同一发明构思，本发明实施例第三方面提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如本发明第一方面所述的方法中的步骤。
基于同一发明构思，本发明实施例第四方面提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如本发明第一方面所述的方法中的步骤。
本发明实施例中，从下限应力和上限应力两个方向逼近试样的弹性段，避免了人为主观因素的影响，提高了弹性参数、裂纹闭合应力和起裂应力的精确度。此外，本发明实施例通过双向迭代可以同时获取试样的弹性参数，裂纹闭合应力和起裂应力，提升工作效率。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例的描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的流程图；
图2是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的步骤流程图；
图3是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的示例实施例中页岩试样的应力-应变曲线；
图4是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的示例实施例中页岩试样的应力-应变曲线及第一次计算过程中对应绘制的裂纹体积应变-轴向应变图曲线示意图；
图5是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的示例实施例中页岩试样的应力-应变曲线及第二次计算过程中对应绘制的裂纹体积应变-轴向应变图曲线示意图；
图6是本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代装置的结构框图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
岩石内部微裂纹扩展引起的体积应变εV由弹性体积应变εV,e和裂纹体积应变εV,cr组成，体积应变εV可以通过以下的公式表示：
εV＝εV,e+εV,cr (1)
对于单轴压缩下的各向同性岩样，假设轴向应变为εz，横向应变为ε1，则总的体积应变εV和弹性体积应变εV,e可表示为：
εV＝εz+2ε1 (2)
其中，E为弹性模量，ν为泊松比，σz为施加的轴向应力。由式(1)可知，总的体积应变εV减去弹性体积应变εV,e得到裂纹体积应变εV,cr，代入式(2)、式(3)可得裂纹体积应变εV,cr的计算公式：
根据式(4)可以得到理想岩石裂纹体积应变-轴向应变曲线。
参照图1，示出了本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的流程图，本发明实施例提供的双向逼近确定岩石弹性参数与起裂应力的迭代方法可以由双向逼近确定岩石弹性参数与起裂应力的迭代装置执行，所述双向逼近确定岩石弹性参数与起裂应力的迭代装置可以体现为一个程序或者一个软件。具体，所述方法可以包括以下步骤：
S101，确定确定下限应力σa(k)和上限应力σb(k)。
在具体实施时，可以先对目标岩石取试样进行常温单轴压缩试验，得到该目标岩石对应的应力-应变曲线。然后可以取该应力-应变曲线上峰前的任意一段曲线，将该曲线假设为弹性段，并以此确定假设的岩石弹性段的下限应力σa(k)和上限应力σb(k)，即假设的弹性段的下限(裂纹闭合应力)和上限(裂纹起裂应力)。本发明实施例中，应力-应变曲线上的弹性段与岩石的弹性变形阶段相对应。
在具体实施时，为了节省计算过程，本发明实施例中，基于相关经验研究，在首次计算时，可以假设：初始下限应力σa(1)为所述应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)为应力-应变曲线中峰值应力的一半。
本发明实施例中，在步骤S101之前，还可以采集所述岩石的试样的常温单轴压缩试验数据，得到所述应力-应变曲线。
S102，根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk。
本发明实施例中，可以根据下限应力σa(k)和上限应力σb(k)、以及试验得到的应力-应变曲线确定出弹性参数，具体的，采用σa(k)和σb(k)之间的割线值作为弹性模量Ek，泊松比＝-横向应变差/轴向应变差。
S103，根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)。
具体的，本发明实施例中，可以根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和公式(4)确定裂纹体积应变随轴向应变变化的数据，并绘制得到裂纹体积应变随轴向应变变化的曲线。基于该曲线，可以得到裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)。
从而，本发明实施例中，可以先基于假设的下限应力σa(k)和上限应力σb(k)，确定弹性参数，再基于弹性参数和裂纹体积应变法得出裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)。从而判断用于计算弹性参数的下限应力σa(k)和上限应力σb(k)，分别与基于弹性参数和裂纹体积应变法得出裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)之间的误差，以确定假设的弹性段是否正确。
S104，在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新。
本发明实施例中，在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，表示S101中假设的下限应力σa(k)和上限应力σb(k)与计算得到的裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)偏差较大，相应的步骤S101中假设的弹性段并不是应力-应变曲线实际的弹性阶段，需要对假设的弹性段进行迭代更新，并重新计算，以逐步逼近应力-应变曲线的弹性阶段。
本发明实施例中，在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)均满足容许误差的情况下，表明假设的弹性段与实际弹性段偏差较小，裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)被成功获取，可以结束计算。
具体的，本发明实施例中，所述步骤S104之前可以对所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)是否满足容许误差进行判断，本发明实施例中，对所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)是否满足容许误差进行判断，是指判断裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差是否满足容许误差，以及判断裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差是否满足容许误差。
具体的判断步骤可以包括：
S1，通过下式(5)计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
S2，通过下式(6)计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
S3，基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差。
具体的，可以采用下式(7)判断所述裂纹闭合应力σcc(k)是否满足容许误差，采用下式(8)判断裂纹起裂应力σci(k)是否满足容许误差：
-χcc≤χcc(k)≤χcc (7)
-χci≤χci(k)≤χci (8)
本发明实施例中，χcc为容许闭合应力误差，χci为容许起裂应力误差，可以根据实际精度要求设置，通常可以设置为1％或2％。
本发明实施例中，基于上述判断步骤，具体的，在步骤S104中包括以下两部分：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；在χcc(k)不大于0的情况下，采用下限应力σa(k)作为下一次迭代计算的下限应力。
本发明实施例中，若由式(5)计算出的误差χcc(k)大于0，说明裂纹闭合应力σcc(k)在σa(k)的上部，则认为σcc(k)更靠近应力-应变曲线的弹性段，将裂纹闭合应力σcc(k)作为假设的下限应力进行下一次迭代计算，以逐步逼近实际的弹性段。若由式(5)计算出的误差χcc(k)不大于0，说明裂纹闭合应力σcc(k)在σa(k)的下部，则认为σa(k)更靠近应力-应变曲线的弹性段，将下限应力σa(k)作为假设的下限应力继续进行下一次迭代计算。
在χci(k)大于0的情况下，采用上限应力σb(k)作为下一次迭代计算的上限应力；在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
本发明实施例中，若由式(6)计算出的误差χci(k)大于0，说明裂纹起裂应力σci(k)在σb(k)的上部，则认为σb(k)更靠近应力-应变曲线的弹性段，将上限应力σb(k)作为假设的上限应力继续进行下一次迭代计算。若由式(6)计算出的误差χci(k)不大于0，说明裂纹起裂应力σci(k)在σb(k)的下部，则认为裂纹起裂应力σci(k)更靠近应力-应变曲线的弹性段，将计算得到的裂纹起裂应力σci(k)作为假设的上限应力进行下一次迭代计算，以逐步逼近实际的弹性段。
由此，本发明实施例中，通过对假设的下限应力和上限应力的不断迭代更新，可以从两端双向调整，逐渐逼近实际的弹性段，并同时确定出弹性参数和裂纹起裂应力。
S105，根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力。
参照图2，示出了本发明实施例的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法的步骤流程图，所述方法具体可以包括：
步骤1，定义σa(k)和σb(k)作为计算弹性参数的下限应力和上限应力。
其中，σa(k)小于σb(k)。首次计算时，k＝1，则初始下限应力σa(1)可以为应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)可以为应力-应变曲线中峰值应力的一半。
步骤2，使用σa(k)和σb(k)计算Ek、νk、σcc(k)、σci(k)。
步骤3，判断-χcc≤χcc(k)≤χcc和-χci≤χci(k)≤χci是否同时成立；若成立，则结束计算；若不成立，转步骤4。
步骤4，判断χcc(k)是否大于0，在χcc(k)大于0的情况下，σa(k+1)＝σcc(k)，在χcc(k)不大于0的情况下，σa(k+1)＝σa(k)；
判断χci(k)是否大于0，在χci(k)大于0的情况下，σb(k+1)＝σb(k)，在χci(k)不大于0的情况下，σb(k+1)＝σci(k)；
返回步骤2，利用σa(k+1)和σb(k+1)进行迭代计算，直至-χcc≤χcc(k)≤χcc和-χci≤χci(k)≤χci同时成立。
为使本领域技术人员更好地理解本发明，以下通过具体的实施例来说明本发明提供的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品。
本发明实施例中的岩石为层理面角度β＝0°(即垂直于轴向加载方向)的页岩，按照国际岩石力学学会(ISRM)标准将页岩加工成高100mm，长和宽均为50mm的长方体试样，对试样进行常温单轴压缩试验，加载速率为0.002mm/s，直至岩样破坏，试验结束。试验所得的页岩试样的应力-应变曲线如图3所示，图3中σa(1)和σb(1)分别是下限应力σa(k)和上限应力σb(k)的起始值。
步骤1：确定用于计算弹性参数的初始下限应力σa(1)和初始上限应力σb(1)。σa(1)取应力-应变曲线0.01％横向应变所对应的应力，σa(1)＝13.52MPa；σb(1)取峰值应力的一半，σb(1)＝60.01MPa。
步骤2：计算σa(1)至σb(1)应力阶段间的割线值作为弹性模量E1，E1＝20.88GPa，泊松比ν1＝0.287，代入公式(4)计算试样的裂纹体积应变，绘制轴向应变-裂纹体积应变曲线，得出裂纹闭合应力σcc(1)＝25.53MPa，起裂应力σci(1)＝46.12MPa，如图4所示，其中，图4中，σcc(1)和σci(1)分别是通过裂纹体积应变法获取的裂纹闭合应力和裂纹起裂应力。
步骤3：根据公式(5)和(6)计算出χcc(1)＝88.83％，χci(1)＝-45.12％，设置裂纹闭合应力和裂纹起裂应力的容许误χcc和χcc都为2％，由于χcc(1)大于2％，χcc小于-2％，不满足公式(7)和公式(8)，不在容许误差范围之内，需要继续进行迭代计算，进入步骤4。
步骤4：由于χcc(1)＝88.83％>0，因此χcc(1)比σa(1)更接近试样的弹性段，使σa(2)＝σcc(1)＝25.53MPa；由于χci(1)＝-45.12％<0，因此σci(1)比σb(1)更接近试样的弹性段，使σb(2)＝σci(1)＝46.12MPa，重复步骤2和步骤3进行弹性参数、闭合应力和起裂应力的计算、以及容许误差的判断，为了便于理解，下面以步骤2-2和步骤3-2进行表示第二次计算过程。
步骤2-2：按照步骤2的运算规则，计算σa(2)与σb(2))间的弹性模量E2＝20.99GPa，泊松比ν2＝0.286，代入公式(4)计算试样的裂纹体积应变，采用裂纹体积应变法绘制试样轴向应变-裂纹体积应变曲线，得出闭合应力σcc(2)＝25.36MPa，起裂应力σci(2)＝45.67MPa，如图5所示，其中，σcc(2)和σci(2)分别是基于迭代后的弹性参数以及裂纹体积应变法获取的裂纹闭合应力和裂纹起裂应力。
步骤3-2：按照步骤3的运算规则，根据公式(5)和(6)计算出χcc(2)＝-0.66％，χci(2)＝-0.98％，由于设置的裂纹闭合应力和起裂应力的容许误差χcc和χcc都为2％，判断χcc(2)和χci(2)分别满足公式(7)与(8)，成功获取试样的弹性参数E2，ν2，裂纹闭合应力σcc(2)和起裂应力σci(2)。因此，页岩试样的弹性模量E＝20.99GPa，泊松比ν＝0.286，闭合应力σcc＝25.36MPa，起裂应力σci＝45.67MPa，迭代结束。
基于同一发明构思，本发明实施例提供一种双向逼近确定岩石弹性参数与起裂应力的迭代装置，参考图6，图6是本发明实施例提供的双向逼近确定岩石弹性参数与起裂应力的迭代装置的示意图，所述装置包括：
第一确定模块601，用于确定下限应力σa(k)和上限应力σb(k)；
第二确定模块602，用于根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk；
第三确定模块603，用于根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)；
更新模块604，用于在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新；
第四确定模块605，用于根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力。
可选地，在首次计算时，将初始下限应力σa(1)确定为所述应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)确定为应力-应变曲线中峰值应力的一半。
可选地，所述第二确定模块，具有用于：
根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和下式确定裂纹体积应变随轴向应变变化的曲线：
其中，εV,cr为裂纹体积应变，εz为轴向应变，ε1为横向应变，σz为施加的轴向应力。
可选地，所述装置还包括：
第一计算模块，用于通过下式计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
第二计算模块，用于通过下式计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
判断模块，用于基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差。
可选地，所述更新模块，具体用于：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；
在χcc(k)不大于0的情况下，采用所述下限应力σa(k)作为下一次迭代计算的下限应力。
可选地，所述更新模块，具体用于：
在χci(k)大于0的情况下，采用所述上限应力σb(k)作为下一次迭代计算的上限应力；
在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
可选地，所述装置还包括：
采集模块，用于采集所述岩石的试样的常温单轴压缩试验数据，得到所述应力-应变曲线。
对于装置实施例而言，由于其与方法实施例基本相似，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
基于同一发明构思，本发明实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法中的步骤。
基于同一发明构思，本发明实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行计算机程序时实现上述任一实施例所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法中的步骤。
本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
本领域内的技术人员应明白，本发明实施例的实施例可提供为方法、装置、或计算机程序产品。因此，本发明实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明实施例是参照根据本发明实施例的方法、终端设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程双向逼近确定岩石弹性参数与起裂应力的迭代终端设备的处理器以产生一个机器，使得通过计算机或其他可编程双向逼近确定岩石弹性参数与起裂应力的迭代终端设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程双向逼近确定岩石弹性参数与起裂应力的迭代终端设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程双向逼近确定岩石弹性参数与起裂应力的迭代终端设备上，使得在计算机或其他可编程终端设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程终端设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明实施例的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。
最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
以上对本发明所提供的一种双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品，进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种双向逼近确定岩石弹性参数与起裂应力的迭代方法，其特征在于，所述方法包括：
确定下限应力σa(k)和上限应力σb(k)；
根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk；
根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)；
在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新；
根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力；
根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，包括：
根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和下式确定裂纹体积应变随轴向应变变化的曲线：
其中，εV,cr为裂纹体积应变，εz为轴向应变，ε1为横向应变，σz为施加的轴向应力；
所述方法还包括：
通过下式计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
通过下式计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差；
根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，包括：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；
在χcc(k)不大于0的情况下，采用所述下限应力σa(k)作为下一次迭代计算的下限应力；
根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新，包括：
在χci(k)大于0的情况下，采用所述上限应力σb(k)作为下一次迭代计算的上限应力；
在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
2.根据权利要求1所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法，其特征在于，在首次计算时，将初始下限应力σa(1)确定为所述应力-应变曲线0.01％横向应变所对应的应力，初始上限应力σb(1)确定为应力-应变曲线中峰值应力的一半。
3.根据权利要求1或2所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法，其特征在于，所述方法还包括：
采集所述岩石的试样的常温单轴压缩试验数据，得到所述应力-应变曲线。
4.一种双向逼近确定岩石弹性参数与起裂应力的迭代装置，其特征在于，所述装置包括：
第一确定模块，用于确定下限应力σa(k)和上限应力σb(k)；
第二确定模块，用于根据所述下限应力σa(k)和所述上限应力σb(k)确定岩石的弹性参数，所述弹性参数包括：弹性模量Ek和泊松比νk；
第三确定模块，用于根据应力-应变曲线、弹性模量Ek和所述泊松比νk确定裂纹体积应变随轴向应变变化的曲线，根据所述裂纹体积应变随轴向应变变化的曲线，确定裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)；
更新模块，用于在所述裂纹闭合应力σcc(k)和裂纹起裂应力σci(k)任一不满足容许误差的情况下，根据计算得到的裂纹闭合应力σcc(k)，对所述下限应力σa(k)进行迭代更新，根据计算得到的裂纹起裂应力σci(k)，对所述上限应力σb(k)进行迭代更新；
第四确定模块，用于根据迭代更新后的下限应力和上限应力重新确定弹性参数，并确定新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)，直至迭代后计算得到的新的裂纹闭合应力σcc(k+1)和新的裂纹起裂应力σci(k+1)均满足所述容许误差，将最后一次迭代计算出的弹性参数和裂纹起裂应力确定为确定岩石的弹性参数与起裂应力；
所述第二确定模块，具有用于：
根据所述应力-应变曲线、所述弹性模量Ek、所述泊松比νk和下式确定裂纹体积应变随轴向应变变化的曲线：
其中，εV,cr为裂纹体积应变，εz为轴向应变，ε1为横向应变，σz为施加的轴向应力；
所述装置还包括：
第一计算模块，用于通过下式计算裂纹体积应变法得出的裂纹闭合应力σcc(k)与用于计算弹性参数的下限应力σa(k)之间的误差χcc(k)：
第二计算模块，用于通过下式计算裂纹体积应变法得出的裂纹起裂应力σci(k)与用于计算弹性参数的上限应力σb(k)之间的误差χci(k)：
判断模块，用于基于χcc(k)判断所述裂纹闭合应力σcc(k)是否满足容许误差，基于χci(k)判断裂纹起裂应力σci(k)是否满足容许误差；
所述更新模块，具体用于：
在χcc(k)大于0的情况下，将计算得到的裂纹闭合应力σcc(k)作为下一次迭代计算的下限应力；
在χcc(k)不大于0的情况下，采用所述下限应力σa(k)作为下一次迭代计算的下限应力；
所述更新模块，具体用于：
在χci(k)大于0的情况下，采用所述上限应力σb(k)作为下一次迭代计算的上限应力；
在χci(k)不大于0的情况下，将计算得到的裂纹起裂应力σci(k)作为下一次迭代计算的上限应力。
5.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-3任一项所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法的步骤。
6.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1-3任一项所述的双向逼近确定岩石弹性参数与起裂应力的迭代方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310339341.0A CN116296820B (zh) | 2023-03-31 | 2023-03-31 | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310339341.0A CN116296820B (zh) | 2023-03-31 | 2023-03-31 | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116296820A CN116296820A (zh) | 2023-06-23 |
| CN116296820B true CN116296820B (zh) | 2024-04-05 |
Family
ID=86828716
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310339341.0A Active CN116296820B (zh) | 2023-03-31 | 2023-03-31 | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116296820B (zh) |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2682522A1 (fr) * | 1991-10-11 | 1993-04-16 | Sgs Thomson Microelectronics | Procede pour verifier le contenu apres effacement d'une memoire permanente effacable, notamment de type eprom, dispositif pour sa mise en óoeuvre et memoire integrant ce dispositif. |
| CN104359776A (zh) * | 2014-11-14 | 2015-02-18 | 四川大学 | 一种爆炸载荷下i型裂纹动态断裂韧度测试方法 |
| CN104897322A (zh) * | 2015-06-23 | 2015-09-09 | 中国矿业大学 | 一种实时采集煤岩体端面各点垂直应力的承压板 |
| CN106124311A (zh) * | 2016-07-29 | 2016-11-16 | 中国科学院武汉岩土力学研究所 | 基于应变测试判别裂纹扩展演化过程的方法 |
| CN107907409A (zh) * | 2017-11-10 | 2018-04-13 | 中国地质大学（武汉） | 一种确定岩石起裂应力的方法、设备及存储设备 |
| CN108343433A (zh) * | 2018-02-28 | 2018-07-31 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
| CN108614035A (zh) * | 2018-03-23 | 2018-10-02 | 长江水利委员会长江科学院 | 大埋深洞室硬岩岩爆孕灾风险识别方法 |
| CN109117573A (zh) * | 2018-08-28 | 2019-01-01 | 中国海洋大学 | 一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法 |
| CN110017135A (zh) * | 2019-02-15 | 2019-07-16 | 西南石油大学 | 一种裂缝性地层井壁裂缝扩展压力预测方法 |
| CN110031304A (zh) * | 2019-04-25 | 2019-07-19 | 四川大学 | 一种岩石变形破坏模式预测方法 |
| CN110658067A (zh) * | 2019-10-14 | 2020-01-07 | 中国科学院武汉岩土力学研究所 | 一种基于全场应变测量判别岩石裂纹应力门槛值的方法 |
| CN112364535A (zh) * | 2020-11-06 | 2021-02-12 | 华东理工大学 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
| CN112986020A (zh) * | 2021-02-03 | 2021-06-18 | 长安大学 | 一种基于应力与声波变化联合表征岩石渐进破坏的方法 |
| CN113742963A (zh) * | 2021-08-18 | 2021-12-03 | 同济大学 | 数物空间实时融合的足尺整体结构试验方法 |
| CN114139310A (zh) * | 2021-11-25 | 2022-03-04 | 成都飞机工业（集团）有限责任公司 | 一种机加反拉方案变形预测及控制方法 |
| CN114279842A (zh) * | 2021-12-27 | 2022-04-05 | 中国科学院武汉岩土力学研究所 | 一种确定岩石裂纹起裂应力和损伤应力的方法及系统 |
| CN115127914A (zh) * | 2022-06-28 | 2022-09-30 | 生态环境部固体废物与化学品管理技术中心 | 一种构建损伤本构模型分析易泥化矿岩损伤演化的方法 |
| CN115292999A (zh) * | 2022-08-04 | 2022-11-04 | 上海交通大学 | SiCf/SiC超声振动辅助切削有限元仿真方法 |
| CN115639057A (zh) * | 2022-10-14 | 2023-01-24 | 中国矿业大学 | 岩石爆破等效加载应力曲线的反演方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9029523B2 (en) * | 2000-04-26 | 2015-05-12 | Ceres, Inc. | Promoter, promoter control elements, and combinations, and uses thereof |
| WO2006006412A1 (ja) * | 2004-07-09 | 2006-01-19 | Kyushu Tlo Company, Limited | 疲労き裂成長曲線の推定法、推定プログラムおよび推定装置 |
| US20160370499A1 (en) * | 2013-07-08 | 2016-12-22 | Schlumberger Technology Corporation | Determining Geomechanics Completion Quality |
| CA2945742C (en) * | 2014-06-04 | 2018-10-30 | Halliburton Energy Services, Inc. | Identifying wellbore location based on seismic data |
| CN109885874B (zh) * | 2019-01-11 | 2022-12-23 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
| CN111859739B (zh) * | 2020-07-02 | 2022-09-20 | 大连理工大学 | 一种Sneddon模型拟合细胞弹性模量的修正方法 |
| US20220341835A1 (en) * | 2021-04-22 | 2022-10-27 | Saudi Arabian Oil Company | Methods and Systems for Determining Reservoir and Fracture Properties |
-
2023
- 2023-03-31 CN CN202310339341.0A patent/CN116296820B/zh active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2682522A1 (fr) * | 1991-10-11 | 1993-04-16 | Sgs Thomson Microelectronics | Procede pour verifier le contenu apres effacement d'une memoire permanente effacable, notamment de type eprom, dispositif pour sa mise en óoeuvre et memoire integrant ce dispositif. |
| CN104359776A (zh) * | 2014-11-14 | 2015-02-18 | 四川大学 | 一种爆炸载荷下i型裂纹动态断裂韧度测试方法 |
| CN104897322A (zh) * | 2015-06-23 | 2015-09-09 | 中国矿业大学 | 一种实时采集煤岩体端面各点垂直应力的承压板 |
| CN106124311A (zh) * | 2016-07-29 | 2016-11-16 | 中国科学院武汉岩土力学研究所 | 基于应变测试判别裂纹扩展演化过程的方法 |
| CN107907409A (zh) * | 2017-11-10 | 2018-04-13 | 中国地质大学（武汉） | 一种确定岩石起裂应力的方法、设备及存储设备 |
| CN108343433A (zh) * | 2018-02-28 | 2018-07-31 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
| CN108614035A (zh) * | 2018-03-23 | 2018-10-02 | 长江水利委员会长江科学院 | 大埋深洞室硬岩岩爆孕灾风险识别方法 |
| CN109117573A (zh) * | 2018-08-28 | 2019-01-01 | 中国海洋大学 | 一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法 |
| CN110017135A (zh) * | 2019-02-15 | 2019-07-16 | 西南石油大学 | 一种裂缝性地层井壁裂缝扩展压力预测方法 |
| CN110031304A (zh) * | 2019-04-25 | 2019-07-19 | 四川大学 | 一种岩石变形破坏模式预测方法 |
| CN110658067A (zh) * | 2019-10-14 | 2020-01-07 | 中国科学院武汉岩土力学研究所 | 一种基于全场应变测量判别岩石裂纹应力门槛值的方法 |
| CN112364535A (zh) * | 2020-11-06 | 2021-02-12 | 华东理工大学 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
| CN112986020A (zh) * | 2021-02-03 | 2021-06-18 | 长安大学 | 一种基于应力与声波变化联合表征岩石渐进破坏的方法 |
| CN113742963A (zh) * | 2021-08-18 | 2021-12-03 | 同济大学 | 数物空间实时融合的足尺整体结构试验方法 |
| CN114139310A (zh) * | 2021-11-25 | 2022-03-04 | 成都飞机工业（集团）有限责任公司 | 一种机加反拉方案变形预测及控制方法 |
| CN114279842A (zh) * | 2021-12-27 | 2022-04-05 | 中国科学院武汉岩土力学研究所 | 一种确定岩石裂纹起裂应力和损伤应力的方法及系统 |
| CN115127914A (zh) * | 2022-06-28 | 2022-09-30 | 生态环境部固体废物与化学品管理技术中心 | 一种构建损伤本构模型分析易泥化矿岩损伤演化的方法 |
| CN115292999A (zh) * | 2022-08-04 | 2022-11-04 | 上海交通大学 | SiCf/SiC超声振动辅助切削有限元仿真方法 |
| CN115639057A (zh) * | 2022-10-14 | 2023-01-24 | 中国矿业大学 | 岩石爆破等效加载应力曲线的反演方法 |
Non-Patent Citations (5)
| Title |
|---|
| "Elasticity and Characteristic Stress Thresholds of Shale under Deep In Situ Geological Conditions";Nie, XF et al;《MATERIALS》;第16卷(第19期);第1-22页 * |
| "EXPERIMENTAL AND NUMERICAL SIMULATION OF THE MICROCRACK COALESCENCE MECHANISM IN ROCK-LIKE MATERIALS";Haeri, H et al;《STRENGTH OF MATERIALS》;第47卷(第5期);第740-754页 * |
| "不同赋存深度岩石力学行为差异及本构模型研究";张茹等;《中国科学基金》;第36卷(第6期);第1008-1015页 * |
| "拉剪作用下类岩石试样的力学特性与损伤破坏机制";周超彪等;《岩土力学》;第42卷(第12期);第3335-3344页 * |
| "裂隙砂岩破坏力学实验与数值模拟研究";付金磊;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》.(第6期);第1-89页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116296820A (zh) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gao et al. | A non-coaxial critical-state model for sand accounting for fabric anisotropy and fabric evolution | |
| Liu et al. | Constitutive modeling of coarse-grained materials incorporating the effect of particle breakage on critical state behavior in a framework of generalized plasticity | |
| Geranmayeh Vaneghi et al. | Strength and damage response of sandstone and granodiorite under different loading conditions of multistage uniaxial cyclic compression | |
| Basu et al. | Shaft resistance and setup factors for piles jacked in clay | |
| Sakurai et al. | Back analysis for tunnel engineering as a modern observational method | |
| Jivkov et al. | Pore space and brittle damage evolution in concrete | |
| Walton | Initial guidelines for the selection of input parameters for cohesion-weakening-friction-strengthening (CWFS) analysis of excavations in brittle rock | |
| CN107340183B (zh) | 结构性软土次固结系数描述方法 | |
| Zhang et al. | Elasto-plastic model of structured marine clay under general loading conditions | |
| Motamedi et al. | Numerical simulation of mixed mode (I and II) fracture behavior of pre-cracked rock using the strong discontinuity approach | |
| Krasiński et al. | Static load test on instrumented pile–field data and numerical simulations | |
| Fang et al. | Geomechanical constitutive modelling of gas hydrate-bearing sediments by a state-dependent multishear bounding surface model | |
| CN105628486A (zh) | 测量页岩岩石力学性质的方法及装置 | |
| Li et al. | Unified solution to drained expansion of a spherical cavity in clay and sand | |
| Liu et al. | Variable parameter creep model based on the separation of viscoelastic and viscoplastic deformations | |
| Karlsson et al. | Modeling aging of displacement piles in natural soft clay | |
| Peng et al. | Analysis of the effect of porosity in concrete under compression based on DIP technology | |
| CN116296820B (zh) | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 | |
| Kamdem et al. | New description of the mechanical creep response of rocks by fractional derivative theory | |
| Liu et al. | Poroelastic response of spherical indentation into a half space with an impermeable surface via step displacement | |
| Titscher et al. | Implicit–explicit integration of gradient-enhanced damage models | |
| Wang et al. | Transversely isotropic creep model for rocks | |
| CN111783332B (zh) | 压硬性非线性变化和剪缩突变特性材料振动累积变形的有限元仿真方法 | |
| Li et al. | Limit deformation analysis of unsaturated expansive soils during wetting and drying cycles | |
| Peng et al. | Research on rock creep characteristics based on the fractional calculus meshless method |
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 |