CN116579068B - 一种活动断层蠕滑变形量的确定方法和确定装置 - Google Patents
一种活动断层蠕滑变形量的确定方法和确定装置 Download PDFInfo
- Publication number
- CN116579068B CN116579068B CN202310848055.7A CN202310848055A CN116579068B CN 116579068 B CN116579068 B CN 116579068B CN 202310848055 A CN202310848055 A CN 202310848055A CN 116579068 B CN116579068 B CN 116579068B
- Authority
- CN
- China
- Prior art keywords
- target
- fault
- monomer
- stiffness
- time period
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000000178 monomer Substances 0.000 claims abstract description 173
- 239000011435 rock Substances 0.000 claims abstract description 108
- 230000006870 function Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 6
- 241000923606 Schistes Species 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229920001971 elastomer Polymers 0.000 description 1
- 239000000806 elastomer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 239000013598 vector Substances 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
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computer Hardware Design (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Geophysics (AREA)
- Food Science & Technology (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本申请提供了一种活动断层蠕滑变形量的确定方法和确定装置，包括：获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量。通过所述确定方法和确定装置，能够使得确定出的活动断层的蠕滑变形量更加准确。
Description
技术领域
本申请涉及地质技术领域，尤其是涉及一种活动断层蠕滑变形量的确定方法和确定装置。
背景技术
活动断层是指正在移动或者历史时期曾有移动或者在未来有移动倾向的断层。通常，活动断层在某一时间段内会存在蠕滑变形的现象，导致在这一时间段内存在一个蠕滑变形量。对于某些工程项目来说，如果能够确定出预定时间段内的蠕滑变形量，就能够方便后续工程项目的顺利进行。例如，对于隧道工程来说，许多隧道在建设时不可避免会穿过活动断层，但是，由于活动断层在一个时间段内存在蠕滑变形的现象，因此，在这一个时间段内会存在隧道空间被挤压的风险。如果能够确定出活动断层在某一时间段内的蠕滑变形量，就能够在建设隧道之前预设相应的设防量，从而避免出现隧道空间被挤压的情况。
现有技术中在确定活动断层的蠕滑变形量时，一般是将活动断层的刚度看作是定值，基于这一固定的刚度确定出活动断层的蠕滑变形量。但是，现有技术中的这种方法，忽略了在一个时间段内活动断层的刚度会随时间变化这一实际情况，从而导致确定出的活动断层的蠕滑变形量不准确。
发明内容
有鉴于此，本申请的目的在于提供一种活动断层蠕滑变形量的确定方法和确定装置，能够使得确定出的活动断层的蠕滑变形量更加准确。
第一方面，本申请实施例提供了一种活动断层蠕滑变形量的确定方法，所述确定方法包括：
获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；
基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；
基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量。
可选地，所述任一目标岩石单体在目标时间段的单体参数数据包括：
所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的。
可选地，所述基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据，包括：
分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数。
可选地，所述分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数，包括：
分别获取与所述单体参数数据中的所述第一单体刚度对应的第一预设量级、与所述第二单体刚度对应的第二预设量级和与所述目标岩石单体的单体弹性系数对应的第三预设量级；
将所述第一单体刚度与所述第一预设量级相乘，确定出所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度；
将所述第二单体刚度与所述第二预设量级相乘，确定出所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度；
将所述目标岩石单体的单体弹性系数与所述第三预设量级相乘，确定出所述目标活动断层的断层弹性系数。
可选地，所述预先构建的与目标活动断层对应的目标刚度公式为：
其中，表示目标活动断层在所述目标时间段的初始时刻的第一断层刚度；/>表示目标活动断层在所述目标时间段的结束时刻的第二断层刚度；/>表示所述目标时间段的时长；/>表示所述目标活动断层的断层弹性系数；/>表示所述目标活动断层的刚度。
第二方面，本申请实施例提供了一种活动断层蠕滑变形量的确定装置，所述确定装置包括：
单体参数获取模块，用于获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；
断层参数获取模块，用于基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
刚度确定模块，用于基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；
蠕滑变形量确定模块，用于基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量。
可选地，所述任一目标岩石单体在目标时间段的单体参数数据包括：
所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的。
可选地，所述断层参数获取模块，具体用于：
分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数。
第三方面，本申请实施例提供了一种电子设备，包括：处理器、存储器和总线，所述存储器存储有所述处理器可执行的机器可读指令，当电子设备运行时，所述处理器与所述存储器之间通过总线通信，所述机器可读指令被所述处理器执行时执行如上述的活动断层蠕滑变形量的确定方法的步骤。
第四方面，本申请实施例提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行如上述的活动断层蠕滑变形量的确定方法的步骤。
本申请实施例提供的一种活动断层蠕滑变形量的确定方法和确定装置，利用预先构建目标活动断层的目标刚度公式，确定出目标活动断层在所述目标时间段的刚度，使得确定出的目标活动断层在所述目标时间段的刚度为随时间变化的值，从而使得基于随时间变化的刚度确定出的活动断层的蠕滑变形量更加准确。
为使本申请的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1示出了本申请示例性实施例提供的一种活动断层蠕滑变形量的确定方法的流程图；
图2a示出了本申请示例性实施例提供的一种拟合曲线的示意图；
图2b示出了本申请另一示例性实施例提供的一种拟合曲线的示意图；
图3示出了本申请示例性实施例提供的一种活动断层蠕滑变形量的确定装置的结构示意图；
图4示出了本申请示例性实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的每个其他实施例，都属于本申请保护的范围。
现有技术中在确定活动断层的蠕滑变形量时，一般是将活动断层的刚度看作是定值，基于这一固定的刚度确定出活动断层的蠕滑变形量。但是，现有技术中的这种方法，忽略了在一个时间段内活动断层的刚度会随时间变化这一实际情况，从而导致确定出的活动断层的蠕滑变形量不准确。
基于此，本申请实施例提供了一种活动断层蠕滑变形量的确定方法和确定装置，能够使得确定出的活动断层的蠕滑变形量更加准确。
请参阅图1，图1示出了本申请示例性实施例提供的一种活动断层蠕滑变形量的确定方法的流程图。
如图1中所示，本申请实施例提供的一种活动断层蠕滑变形量的确定方法，包括以下步骤：
S1、获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据。
这里，需要说明的是，目标活动断层中包括多个岩石单体。其中，任一目标岩石单体可以是指多个岩石单体中的任意一个岩石单体。目标时间段是根据实际情况确定的，例如，目标时间段可以为1天、1年、3年等，在进行试验时，目标时间段可以是1天，在工程实际应用中，目标时间段可以是1年和3年等。
作为示例，所述任一目标岩石单体在目标时间段的单体参数数据可以包括：
所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的。
作为一示例，所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度可以是通过在所述目标时间段的结束时刻对所述目标岩石单体进行测量得到的。
作为另一示例，所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度可以是基于所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度被预测得到的。
这里，可以利用现有技术中的任意一种方式基于所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度预测所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度。
例如，可以将所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度和预设刚度变化量的差值，确定为所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度。其中，预设刚度变化量与目标岩石单体的岩石种类和目标时间段的时长有关，例如，当目标岩石单体的岩石种类为泥岩，目标时间段的时长为1年时，预设刚度变化量为1×pa。
这里，每个岩石种类都对应一个单体弹性系数，例如，当岩石种类为泥岩时，单体弹性系数可以为240；当岩石种类为绿片岩时，单体弹性系数可以为250。
S2、基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
作为示例，在步骤S2中，可以分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数。
具体地，可以首先，分别获取与所述单体参数数据中的所述第一单体刚度对应的第一预设量级、与所述第二单体刚度对应的第二预设量级和与所述目标岩石单体的单体弹性系数对应的第三预设量级；
这里，第一预设量级、第二预设量级和第三预设量级可以是根据实际情况确定的，例如，第一预设量级可以是，第二预设量级可以是/>，第三预设量级可以是/>。
然后，将所述第一单体刚度与所述第一预设量级相乘，确定出所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度；
将所述第二单体刚度与所述第二预设量级相乘，确定出所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度；
将所述目标岩石单体的单体弹性系数与所述第三预设量级相乘，确定出所述目标活动断层的断层弹性系数。
S3、基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；
其中，所述目标刚度公式为目标活动断层的刚度随时间变化的函数。
这里，所述预先构建的与目标活动断层对应的目标刚度公式为：
（1）
其中，表示目标活动断层在所述目标时间段的初始时刻的第一断层刚度；/>表示目标活动断层在所述目标时间段的结束时刻的第二断层刚度；/>表示所述目标时间段的时长；/>表示所述目标活动断层的断层弹性系数；/>表示所述目标活动断层的刚度。
这里，确定出的所述目标活动断层的刚度为随时间变化的值，并非是固定的值。
下面，将对上述公式（1）所对应的目标刚度公式的验证过程进行详细说明。
首先需要说明的是，公式（1）表示的是目标活动断层对应的目标刚度公式，可以理解，由于目标活动断层在所述目标时间段的断层参数数据是基于岩石单体在目标时间段的单体参数数据扩大量级得到的，实际上断层参数数据和单体参数数据的内容不变，因此，可以假设岩石单体对应的单体刚度公式为：
公式（2）
其中，表示岩石单体在目标时间段的初始时刻的第一单体刚度；/>表示岩石单体在目标时间段的结束时刻的第二单体刚度；/>表示目标时间段的时长；/>表示岩石单体的单体弹性系数；/>表示岩石单体的刚度。
具体地，由Maxwell（麦克斯韦）蠕变模型可知，任一活动断层中的任一岩石单体可以看作是由一个弹性体和一个黏性体串联组成的岩石体，因此，任一岩石单体在任一目标时间段内的蠕滑变形量可由下面公式（3）所对应的线性蠕变方程来表示：
公式（3）
其中，为目标时间段的时长；/>为岩石单体在目标时间段对应的刚度定值；/>为黏性系数；/>为蠕变应力；/>为岩石单体在目标时间段的蠕滑变形量。
将上述公式（2）所对应单体刚度公式代入到公式（3）所对应的蠕变模型方程中，具体地，利用公式（2）中的所对应的方程替换公式（3）中的E，可以得到下面公式（4）所对应的岩石单体的非线性蠕变方程：
公式（4）
在岩石蠕变试验中，当蠕变应力为10MPa时，岩石单体在1000h时间内没有出现明显的加速蠕变过程，在同等工况条件下加载时间至1900h，岩石单体出现明显的加速蠕变过程。即，当蠕变应力太小时，岩石单体出现加速蠕变过程的现象不明显，因此，在验证过程中应该使用较大的蠕变应力，至少要保证蠕变应力大于10Mpa。
作为示例，下面将以两个实例来说明对岩石单体的非线性蠕变方程的验证结果。其中，为了更好地验证非线性蠕变方程的普适性，第一个实例中采用的岩石单体为泥岩，第二个实例中采用的岩石单体为绿片岩。其中，第一个实例中蠕变应力为13Mpa，第二个实例中蠕变应力为100Mpa。
请参阅表1和图2a，表1示出了本申请示例性实施例提供的一种当岩石单体为泥岩时该泥岩的参数数据，图2a示出了本申请示例性实施例提供的一种拟合曲线的示意图，其中，图2a为将表1中的参数数据代入到上述公式（4）中所对应的拟合曲线。其中，未改进模型为公式（3）所对应的线性蠕变方程，改进模型为公式（4）所对应的非线性蠕变方程。
表1：
通过图2a可以看出，当岩石单体为泥岩时，通过岩石单体的非线性蠕变方程得到的应变（即蠕滑变形量）与岩石单体的线性蠕变方程得到的应变相比，岩石单体的非线性蠕变方程得到的应变与对该泥岩进行试验得到的实验数据更为接近。
请参阅表2和图2b，表2示出了本申请示例性实施例提供的一种当岩石单体为绿片岩时该绿片岩的参数数据，图2b示出了本申请另一示例性实施例提供的一种拟合曲线的示意图，其中，图2b示出了将表2中的参数数据代入到上述公式（4）中所对应的拟合曲线。
表2
通过图2b可以看出，当岩石单体为绿片岩时，通过岩石单体的非线性蠕变方程得到的应变（即蠕滑变形量）与岩石单体的线性蠕变方程得到的应变相比，岩石单体的非线性蠕变方程得到的应变与该泥岩进行试验得到的实验数据更为接近。
综上所述，通过公式（4）所对应的岩石单体的非线性蠕变方程能够较为准确地确定岩石单体的应变，从而证明公式（2）所对应的岩石单体对应的单体刚度公式较为准确合理，从而证明预先构建的公式（1）所对应的目标活动断层对应的目标刚度公式较为准确合理。
S4、基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量。
作为示例，在步骤S4中，可以基于所述目标活动断层在所述目标时间段的刚度，利用三维弹性半空间方程，确定目标活动断层在所述目标时间段的蠕滑变形量。
这里，三维弹性半空间方程为现有的方程，三维弹性半空间方程不是解析解，而是首先在Smith 3D-体力模型示意波数域解弹性方程，然后进行傅里叶逆变换得到空间域的解，该方案在保持计算精度与数值解一致的前提下大大地提高了计算效率。3D-体力模型能够得到在均匀、各向同性的弹性介质中，由直立断层上的体力矢量作用下，断层下界/>和上界/>之间引起的位移。这里，通过非线性蠕变方程（与公式（4）中的岩石单体的非线性蠕变方程对应的活动断层的非线性蠕变方程）计算蠕滑变形量时所需要施加的应力方向是一维的，而通过3D-体力模型计算蠕滑变形量时所需要施加的应力方向是三维的，因此，通过3D-体力模型计算蠕滑变形量相比较于通过非线性蠕变方程计算蠕滑变形量能够更加符合实际情况。
例如，三维弹性半空间方程可以表示如下：
其中，
；
本申请实施例提供的一种活动断层蠕滑变形量的确定方法，利用预先构建目标活动断层的目标刚度公式，确定出目标活动断层在所述目标时间段的刚度，使得确定出的目标活动断层在所述目标时间段的刚度为随时间变化的值，从而使得基于随时间变化的刚度确定出的活动断层的蠕滑变形量更加准确。
请参阅图3，图3示出了本申请示例性实施例提供的一种活动断层蠕滑变形量的确定装置的结构示意图。
如图3中所示，所述确定装置300包括：
单体参数获取模块310，用于获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；
断层参数获取模块320，用于基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
刚度确定模块330，用于基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；
蠕滑变形量确定模块340，用于基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量。
可选地，所述任一目标岩石单体在目标时间段的单体参数数据包括：
所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的。
可选地，所述断层参数获取模块320，具体用于：
分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数。
可选地，所述断层参数获取模块320，具体用于：
分别获取与所述单体参数数据中的所述第一单体刚度对应的第一预设量级、与所述第二单体刚度对应的第二预设量级和与所述目标岩石单体的单体弹性系数对应的第三预设量级；
将所述第一单体刚度与所述第一预设量级相乘，确定出所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度；
将所述第二单体刚度与所述第二预设量级相乘，确定出所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度；
将所述目标岩石单体的单体弹性系数与所述第三预设量级相乘，确定出所述目标活动断层的断层弹性系数。
可选地，所述预先构建的与目标活动断层对应的目标刚度公式为：
其中，表示目标活动断层在所述目标时间段的初始时刻的第一断层刚度；/>表示目标活动断层在所述目标时间段的结束时刻的第二断层刚度；/>表示所述目标时间段的时长；/>表示所述目标活动断层的断层弹性系数；/>表示所述目标活动断层的刚度。
本申请实施例提供的一种活动断层蠕滑变形量的确定装置，利用预先构建目标活动断层的目标刚度公式，确定出目标活动断层在所述目标时间段的刚度，使得确定出的目标活动断层在所述目标时间段的刚度为随时间变化的值，从而使得基于随时间变化的刚度确定出的活动断层的蠕滑变形量更加准确。
请参阅图4，图4为本申请实施例所提供的一种电子设备的结构示意图。如图4中所示，所述电子设备400包括处理器410、存储器420和总线430。
所述存储器420存储有所述处理器410可执行的机器可读指令，当电子设备400运行时，所述处理器410与所述存储器420之间通过总线430通信，所述机器可读指令被所述处理器410执行时，可以执行如上述方法实施例中的活动断层蠕滑变形量的确定方法的步骤，具体实现方式可参见方法实施例，在此不再赘述。
本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时可以执行如上述方法实施例中的活动断层蠕滑变形量的确定方法的步骤，具体实现方式可参见方法实施例，在此不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（Read-OnlyMemory，ROM）、随机存取存储器（Random Access Memory，RAM）、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上所述实施例，仅为本申请的具体实施方式，用以说明本申请的技术方案，而非对其限制，本申请的保护范围并不局限于此，尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。
Claims (5)
1.一种活动断层蠕滑变形量的确定方法，其特征在于，所述确定方法包括：
获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；所述任一目标岩石单体在目标时间段的单体参数数据包括：所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的；
基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；
基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量；
所述基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据，包括：
分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数；
所述预先构建的与目标活动断层对应的目标刚度公式为：
其中，表示目标活动断层在所述目标时间段的初始时刻的第一断层刚度；/>表示目标活动断层在所述目标时间段的结束时刻的第二断层刚度；/>表示所述目标时间段的时长；表示所述目标活动断层的断层弹性系数；/>表示所述目标活动断层的刚度。
2.根据权利要求1所述的确定方法，其特征在于，所述分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数，包括：
分别获取与所述单体参数数据中的所述第一单体刚度对应的第一预设量级、与所述第二单体刚度对应的第二预设量级和与所述目标岩石单体的单体弹性系数对应的第三预设量级；
将所述第一单体刚度与所述第一预设量级相乘，确定出所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度；
将所述第二单体刚度与所述第二预设量级相乘，确定出所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度；
将所述目标岩石单体的单体弹性系数与所述第三预设量级相乘，确定出所述目标活动断层的断层弹性系数。
3.一种活动断层蠕滑变形量的确定装置，其特征在于，所述确定装置包括：
单体参数获取模块，用于获取目标活动断层中的任一目标岩石单体在目标时间段的单体参数数据；所述任一目标岩石单体在目标时间段的单体参数数据包括：所述目标岩石单体在所述目标时间段的初始时刻的第一单体刚度、所述目标岩石单体在所述目标时间段的结束时刻的第二单体刚度和所述目标岩石单体的单体弹性系数；其中，所述单体弹性系数是基于所述目标岩石单体的岩石种类确定的；
断层参数获取模块，用于基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据；
刚度确定模块，用于基于所述断层参数数据，利用预先构建的与目标活动断层对应的目标刚度公式，确定所述目标活动断层在所述目标时间段的刚度；所述目标刚度公式为目标活动断层的刚度随时间变化的函数；
蠕滑变形量确定模块，用于基于所述目标活动断层在所述目标时间段的刚度，确定所述目标活动断层在所述目标时间段的蠕滑变形量；
所述基于所述单体参数数据，确定所述目标活动断层在所述目标时间段的断层参数数据，包括：
分别基于所述单体参数数据中的所述第一单体刚度、所述第二单体刚度和所述目标岩石单体的单体弹性系数，确定所述目标活动断层在所述目标时间段的初始时刻的第一断层刚度、所述目标活动断层在所述目标时间段的结束时刻的第二断层刚度和所述目标活动断层的断层弹性系数；
所述预先构建的与目标活动断层对应的目标刚度公式为：
其中，表示目标活动断层在所述目标时间段的初始时刻的第一断层刚度；/>表示目标活动断层在所述目标时间段的结束时刻的第二断层刚度；/>表示所述目标时间段的时长；表示所述目标活动断层的断层弹性系数；/>表示所述目标活动断层的刚度。
4.一种电子设备，其特征在于，包括：处理器、存储器和总线，所述存储器存储有所述处理器可执行的机器可读指令，当电子设备运行时，所述处理器与所述存储器之间通过所述总线进行通信，所述机器可读指令被所述处理器运行时执行如权利要求1至2任一所述的活动断层蠕滑变形量的确定方法的步骤。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器运行时执行如权利要求1至2任一所述的活动断层蠕滑变形量的确定方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310848055.7A CN116579068B (zh) | 2023-07-12 | 2023-07-12 | 一种活动断层蠕滑变形量的确定方法和确定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310848055.7A CN116579068B (zh) | 2023-07-12 | 2023-07-12 | 一种活动断层蠕滑变形量的确定方法和确定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116579068A CN116579068A (zh) | 2023-08-11 |
| CN116579068B true CN116579068B (zh) | 2023-10-20 |
Family
ID=87541666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310848055.7A Active CN116579068B (zh) | 2023-07-12 | 2023-07-12 | 一种活动断层蠕滑变形量的确定方法和确定装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116579068B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005061091A (ja) * | 2003-08-14 | 2005-03-10 | Seisuke Okubo | クリープモデルに従う増分型の構成関係式を用いた有限要素解析法 |
| CN110608054A (zh) * | 2019-09-30 | 2019-12-24 | 中铁第四勘察设计院集团有限公司 | 一种全断面法开挖隧道支护结构变形预测方法及装置 |
| CN112213179A (zh) * | 2019-07-09 | 2021-01-12 | 中国石油天然气集团有限公司 | 一种测量井下页岩层理面刚度的方法及装置 |
| CN114398805A (zh) * | 2022-03-25 | 2022-04-26 | 四川省水利水电勘测设计研究院有限公司 | 一种水岩耦合作用下裂隙岩石蠕变模型构建方法及系统 |
-
2023
- 2023-07-12 CN CN202310848055.7A patent/CN116579068B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005061091A (ja) * | 2003-08-14 | 2005-03-10 | Seisuke Okubo | クリープモデルに従う増分型の構成関係式を用いた有限要素解析法 |
| CN112213179A (zh) * | 2019-07-09 | 2021-01-12 | 中国石油天然气集团有限公司 | 一种测量井下页岩层理面刚度的方法及装置 |
| CN110608054A (zh) * | 2019-09-30 | 2019-12-24 | 中铁第四勘察设计院集团有限公司 | 一种全断面法开挖隧道支护结构变形预测方法及装置 |
| CN114398805A (zh) * | 2022-03-25 | 2022-04-26 | 四川省水利水电勘测设计研究院有限公司 | 一种水岩耦合作用下裂隙岩石蠕变模型构建方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 盐岩非线性蠕变损伤本构模型及其工程应用;陈卫忠;王者超;伍国军;杨建平;张保平;;岩石力学与工程学报(03);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116579068A (zh) | 2023-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109460561B (zh) | 采油设备故障诊断方法及装置 | |
| RU2016142123A (ru) | Способ оценки нормальности или ненормальности измеренного значения физического параметра двигателя летательного аппарата | |
| CN103983412A (zh) | 用于振动有限元模型修正航空电子设备工作模态测量方法 | |
| CN116579068B (zh) | 一种活动断层蠕滑变形量的确定方法和确定装置 | |
| CN116659803A (zh) | 基于天平零点在线监测获取连续式风洞气动力载荷的方法 | |
| CN103793032A (zh) | 用于确定上电复位的方法和装置 | |
| CN110309622B (zh) | 一种输电塔结构倒塌分析方法 | |
| JP6641678B2 (ja) | 異常診断装置 | |
| CN109740260B (zh) | 汽轮机转子动平衡处理方法及装置 | |
| CN108009582B (zh) | 一种电子设备标准工作指标设置的方法 | |
| CN111797483B (zh) | 马达均衡电信号的修正方法及设备、计算机可读存储介质 | |
| CN111582463B (zh) | 伺服电机故障识别及模型训练方法、装置、介质及终端 | |
| CN116306109B (zh) | 基于时域模型修正的海上风机土结相互作用状态识别方法 | |
| CN110807283B (zh) | 基于有限元计算模型的多载荷应变测量点获取方法及系统 | |
| CN104359765A (zh) | 测量钞票的弹性模量的方法以及最大弯曲挠度的装置 | |
| CN109101750B (zh) | 一种电子器件湿热环境下的加速寿命预测方法 | |
| CN113704871B (zh) | 车轮弯曲疲劳的确定方法、装置、终端设备及介质 | |
| JP5948998B2 (ja) | 異常診断装置 | |
| CN113419706A (zh) | 一种快速的任意分布随机数产生方法、系统及其检验方法、系统 | |
| CN115062424A (zh) | 航空发动机转子系统的数字孪生体建模方法及数字孪生体 | |
| CN110245425B (zh) | 空调压缩机激励辨识方法以及计算机装置 | |
| CN107958096B (zh) | 电机能耗处理方法及装置 | |
| CN110807284A (zh) | 基于有限元和应变测量的载荷谱识别方法、系统及装置 | |
| CN117174211B (zh) | 一种复合材料的力学性能分析方法、装置、设备及介质 | |
| JP6696742B2 (ja) | 免震部材応答推定装置及び免震部材応答推定方法 |
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 |