CN114692391A - 用于模拟岩溶地区地层的相似试验装置及其制作的方法 - Google Patents
用于模拟岩溶地区地层的相似试验装置及其制作的方法 Download PDFInfo
- Publication number
- CN114692391A CN114692391A CN202210227771.9A CN202210227771A CN114692391A CN 114692391 A CN114692391 A CN 114692391A CN 202210227771 A CN202210227771 A CN 202210227771A CN 114692391 A CN114692391 A CN 114692391A
- Authority
- CN
- China
- Prior art keywords
- karst
- stratum
- module
- evolution
- simulating
- 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.)
- Granted
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 31
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 32
- 238000001514 detection method Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 238000010276 construction Methods 0.000 claims abstract description 21
- 238000002474 experimental method Methods 0.000 claims abstract description 15
- 230000008676 import Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 45
- 230000004992 fission Effects 0.000 claims description 27
- 238000004458 analytical method Methods 0.000 claims description 23
- 238000011161 development Methods 0.000 claims description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- 230000002706 hydrostatic effect Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000010291 electrical method Methods 0.000 claims description 6
- 238000010587 phase diagram Methods 0.000 claims description 6
- 239000013535 sea water Substances 0.000 claims description 6
- 239000007789 gas Substances 0.000 claims description 4
- 239000003345 natural gas Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000007635 classification algorithm Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 239000002689 soil Substances 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 239000004576 sand Substances 0.000 abstract description 2
- 238000009825 accumulation Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- 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
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Food Science & Technology (AREA)
- Algebra (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Operations Research (AREA)
- Remote Sensing (AREA)
- Immunology (AREA)
- Geology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于模拟岩溶地区地层技术领域，公开了一种用于模拟岩溶地区地层的相似试验装置及其制作的方法，所述用于模拟岩溶地区地层的相似试验装置包括：地层倾角检测模块、岩溶地层参数采集模块、参数导入模块、模型构建模块、相似实验模拟模块、塌陷评价模块。本发明通过地层倾角检测模块利用预设窗口长度的深度窗口来其他区相应的演变波形数据，这样相较于现有的地层倾角检测方法，能够明显减少同相轴稀疏性的影响，能够得到更为稳定、精度更高的倾角剖面；同时，通过塌陷评价模块对岩溶上覆砂层塌陷的严重性等级进行判定，在实际工程的施工有很强的指导作用，也为类似岩溶环境提供了计算参考。
Description
技术领域
本发明属于模拟岩溶地区地层技术领域，尤其涉及一种用于模拟岩溶地区地层的相似试验装置及其制作的方法方法。
背景技术
地层结构(stratigraphic texture)是地层序列内岩层的叠覆与堆积型式。通常所说的旋回性、韵律性、加积型、退积型、进积型沉积等，都属于地层结构的概念范畴。地层结构主要是用于研究和描述相当于或略小于一个体系域的地层间隔内岩层的纵、横向总体(或优势)堆积方式，并可根据其成因或形态，将总体退积、加积、进积，或超覆、叠覆、退覆堆积的地层，分别称为退积、加积、进积,或超覆、叠覆、退覆结构的地层。对于地层分析、预测来说，地层结构概念是非常重要的。然而，现有用于模拟岩溶地区地层的相似试验装置及其制作的方法在对低信噪比数据时的地层倾角进行检测时，会出现检测结果不准确的问题；同时，不能对岩溶地面塌陷进行评价。
综上所述，现有技术存在的问题是：现有用于模拟岩溶地区地层的相似试验装置及其制作的方法在对低信噪比数据时的地层倾角进行检测时，会出现检测结果不准确的问题；同时，不能对岩溶地面塌陷进行评价。
发明内容
针对现有技术存在的问题，本发明提供了一种用于模拟岩溶地区地层的相似试验装置及其制作的方法。
本发明是这样实现的，一种用于模拟岩溶地区地层的相似试验装置包括：
地层倾角检测模块、岩溶地层参数采集模块、参数导入模块、模型构建模块、相似实验模拟模块、塌陷评价模块；
地层倾角检测模块，与岩溶地层参数采集模块连接，用于检测岩溶地区地层倾角；
岩溶地层参数采集模块，与地层倾角检测模块、参数导入模块连接，用于采集岩溶地区地层物理参数、力学参数；
参数导入模块，与岩溶地层参数采集模块、模型构建模块连接，用于将采集的岩溶地层参数导入到模型构建模块；
模型构建模块，与参数导入模块、相似实验模拟模块连接，用于通过模型构建程序根据导入的参数构建岩溶地层模型；
相似实验模拟模块，与模型构建模块、塌陷评价模块连接，用于通过模拟程序对模拟岩溶地区地层的相似试验进行模拟；
所述模拟方法包括：
选择岩溶发育且演变资料上BSR特征明显的演变剖面；
确定演变剖面上岩溶地区地层底部界面的双程走时Tsb、观测BSR的双程走时Tbsr和岩溶地区地层顶部边界的双程走时Tsalt；
根据分析区天然气岩溶的气体组分信息，选择合适的岩溶相平衡曲线；
演变剖面上BSR发育深度处的强度值为静水强度；
BSR的深度值通过下面公式求得：
Hbsr＝Vsw×Tbsr/2
式中，Vsw为海水速度，为1500m/s；Tbsr为BSR位置的双程走时，单位是s，通过演变剖面读取；
确定演变剖面上岩溶地区地层底部各点的裂变值T；
通过二维稳态热传导方程，计算演变剖面的裂变场分布：
将岩溶稳定底界移向低的裂变条件；
通过设置不同的裂变梯度，模拟对应的BSR深度，将模拟的BSR位置与演变剖面上的观测BSR进行对比。
进一步，纯岩溶的相平衡公式如下：
In(P)＝a+bT+cT2+dT3+fT4+gT5
式中，P和T分别为岩溶的稳定强度条件和稳定裂变条件，a、b、c、d、f、g为经验常数，分别为a＝-1.94138504464560×105,b＝3.31018213397926×103,c＝-2.25540264493806×101,d＝7.67559117787059×10-2,f＝-1.30465829788791×10-4,g＝8.8606531668757×10-8。
进一步，静水强度通过以下公式求取：
Pbsr＝ρsw g Hbsr
式中，ρsw为海水密度，为1028kg/m3，g是重力加速度，为9.81m/s2，Hbsr是BSR的深度值。
进一步，计算演变剖面的裂变场分布：
式中，T为裂变℃，x为横向距离km，z为垂向距离km，kx为横向导热率W m-1K-1，kz为垂向导热率W m-1K-1；一般认为沉积物各向同性且均质性较强，因此kx＝kz，但沉积物和盐的热导率差别较大，分别设置为2.5W m-1K-1和5.9W m-1K-1。
进一步，岩溶稳定底界移向低的裂变条件：
式中，m为裂变对岩溶稳定条件的影响参数，计算公式如下：
式中，Sw为裂变值。
本发明的另一目的在于提供一种用于模拟岩溶地区地层的相似试验装置的制作的方法包括以下步骤：
步骤一，通过地层倾角检测模块检测岩溶地区地层倾角；通过岩溶地层参数采集模块采集岩溶地区地层物理参数、力学参数；
步骤二，通过参数导入模块将采集的岩溶地层参数导入到模型构建模块；通过模型构建模块利用模型构建程序根据导入的参数构建岩溶地层模型；
步骤三，通过相似实验模拟模块利用模拟程序对模拟岩溶地区地层的相似试验进行模拟；
步骤四，通过塌陷评价模块利用评价程序对岩溶塌陷进行评价。
进一步，所述地层倾角检测模块检测方法如下：
1)通过检测设备根据当前道分析点和各个预设扫描倾角，确定对应于所述各个预设扫描倾角的临近道分析点，根据所述当前道分析点和临近道分析点，分别获取当前道的演变波形数据和对应于所述各个预设扫描倾角的临近道的演变波形数据；
2)沿层位获得演变波形数据，并采用基于线性变换的SSDR算法为演变波形数据降维；利用演变波形数据中的标签数据训练一个距离衡量矩阵；采用半监督Kmeans分类算法对演变波形数据进行分类，生成演变相图；
3)根据所述当前道的演变波形数据、演变相图和对应于所述各个预设扫描倾角的临近道的演变波形数据，确定所述各个预设扫描倾角的相似能量谱；
4)从各个相似能量谱中提取最大相似能量谱，从各个相似能量谱中提取与所述最大相似能量谱临近的若干相似能量谱，根据所述最大相似能量谱和所述若干临近的能量谱确定地层倾角插值曲线，并根据所述地层倾角插值曲线确定地层倾角；
其中，从各个相似能量谱中提取与所述最大相似能量谱相邻的两个相似能量谱；计算所述地层倾角插值曲线的一阶导数；将所述一阶导数等于零时所对应的倾角作为地层倾角；
根据如下公式计算所述地层倾角：
其中，x表示地层倾角，y1表示最大相似能量谱，x1表示最大相似能量谱对应的预设扫描倾角，y0和y2表示与所述最大相似能量谱相邻的两个相似能量谱。
进一步，在当前道中，获取以所述当前道分析点为中心的预设深度窗口的演变波形数据，作为所述当前道的演变波形数据；在各个临近道中，获取以相应临近道分析点为中心的预设深度窗口的演变波形数据，作为所述临近道的演变波形数据。
进一步，所述塌陷评价模块评价方法如下：
(1)通过探地雷达法、瞬变电磁法、跨孔高密度电法和地表高密度电法岩溶进行探测；获取岩溶基本数据以及岩溶换算数据；
(2)建立沙漏型溶洞-地层三维数值模型，引入所获取的数据模拟计算岩溶上覆砂土漏失塌陷过程；
(3)提取岩溶塌陷后的底层沉降剖面图，根据熔岩塌陷对周围底层及建筑构筑物的影响确定岩溶沉降等级。
本发明的优点及积极效果为：本发明通过地层倾角检测模块利用预设窗口长度的深度窗口来其他区相应的演变波形数据，这样相较于现有的地层倾角检测方法，能够明显减少同相轴稀疏性的影响，能够得到更为稳定、精度更高的倾角剖面；同时，通过塌陷评价模块紧密结合沙漏型岩溶塌陷机理，能够更好地模拟出沙漏型岩溶塌陷后地层及周边重要建(构)筑物的变形特性；通过该理论能够对在岩溶发育地区进行城市地下建设过程中遇到沙漏型岩溶塌陷的问题进行计算，对岩溶上覆砂层塌陷的严重性等级进行判定，在实际工程的施工有很强的指导作用，也为类似岩溶环境提供了计算参考。
本发明所述模拟方法模拟效果真实逼真。选择岩溶发育且演变资料上BSR特征明显的演变剖面；
确定演变剖面上岩溶地区地层底部界面的双程走时Tsb、观测BSR的双程走时Tbsr和岩溶地区地层顶部边界的双程走时Tsalt；
根据分析区天然气岩溶的气体组分信息，选择合适的岩溶相平衡曲线；
演变剖面上BSR发育深度处的强度值为静水强度；实现了获取真实数据的依据。
附图说明
图1是本发明实施例提供的用于模拟岩溶地区地层的相似试验装置的制作的方法流程图。
图2是本发明实施例提供的用于模拟岩溶地区地层的相似试验装置结构框图。
图3是本发明实施例提供的地层倾角检测模块检测方法流程图。
图4是本发明实施例提供的塌陷评价模块评价方法流程图。
图2中：1、地层倾角检测模块；2、岩溶地层参数采集模块；3、参数导入模块；4、模型构建模块；5、相似实验模拟模块；6、塌陷评价模块。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效，兹例举以下实施例，并配合附图详细说明如下。
下面结合附图对本发明的结构作详细的描述。
如图1所示，本发明提供的用于模拟岩溶地区地层的相似试验装置的制作的方法包括以下步骤：
S101，通过地层倾角检测模块检测岩溶地区地层倾角；通过岩溶地层参数采集模块采集岩溶地区地层物理参数、力学参数；
S102，通过参数导入模块将采集的岩溶地层参数导入到模型构建模块；通过模型构建模块利用模型构建程序根据导入的参数构建岩溶地层模型；
S103，通过相似实验模拟模块利用模拟程序对模拟岩溶地区地层的相似试验进行模拟；
S104，通过塌陷评价模块利用评价程序对岩溶塌陷进行评价。
如图2所示，本发明实施例提供的用于模拟岩溶地区地层的相似试验装置包括：地层倾角检测模块1、岩溶地层参数采集模块2、参数导入模块3、模型构建模块4、相似实验模拟模块5、塌陷评价模块6。
地层倾角检测模块1，与岩溶地层参数采集模块2连接，用于检测岩溶地区地层倾角；
岩溶地层参数采集模块2，与地层倾角检测模块1、参数导入模块3连接，用于采集岩溶地区地层物理参数、力学参数；
参数导入模块3，与岩溶地层参数采集模块2、模型构建模块4连接，用于将采集的岩溶地层参数导入到模型构建模块4；
模型构建模块4，与参数导入模块3、相似实验模拟模块5连接，用于通过模型构建程序根据导入的参数构建岩溶地层模型；
相似实验模拟模块5，与模型构建模块4、塌陷评价模块6连接，用于通过模拟程序对模拟岩溶地区地层的相似试验进行模拟；
所述模拟方法包括：
选择岩溶发育且演变资料上BSR特征明显的演变剖面；
确定演变剖面上岩溶地区地层底部界面的双程走时Tsb、观测BSR的双程走时Tbsr和岩溶地区地层顶部边界的双程走时Tsalt；
根据分析区天然气岩溶的气体组分信息，选择合适的岩溶相平衡曲线；
演变剖面上BSR发育深度处的强度值为静水强度；
BSR的深度值通过下面公式求得：
Hbsr＝Vsw×Tbsr/2
式中，Vsw为海水速度，为1500m/s；Tbsr为BSR位置的双程走时，单位是s，通过演变剖面读取；
确定演变剖面上岩溶地区地层底部各点的裂变值T；
通过二维稳态热传导方程，计算演变剖面的裂变场分布：
将岩溶稳定底界移向低的裂变条件；
通过设置不同的裂变梯度，模拟对应的BSR深度，将模拟的BSR位置与演变剖面上的观测BSR进行对比。
在本发明中，纯岩溶的相平衡公式如下：
In(P)＝a+bT+cT2+dT3+fT4+gT5
式中，P和T分别为岩溶的稳定强度条件和稳定裂变条件，a、b、c、d、f、g为经验常数，分别为a＝-1.94138504464560×105,b＝3.31018213397926×103,c＝-2.25540264493806×101,d＝7.67559117787059×10-2,f＝-1.30465829788791×10-4,g＝8.8606531668757×10-8。
在本发明中，静水强度通过以下公式求取：
Pbsr＝ρsw g Hbsr
式中，ρsw为海水密度，为1028kg/m3，g是重力加速度，为9.81m/s2，Hbsr是BSR的深度值。
在本发明中，计算演变剖面的裂变场分布：
式中，T为裂变℃，x为横向距离km，z为垂向距离km，kx为横向导热率W m-1K-1，kz为垂向导热率W m-1K-1；一般认为沉积物各向同性且均质性较强，因此kx＝kz，但沉积物和盐的热导率差别较大，分别设置为2.5W m-1K-1和5.9W m-1K-1。
在本发明中，岩溶稳定底界移向低的裂变条件：
式中，m为裂变对岩溶稳定条件的影响参数，计算公式如下：
式中，Sw为裂变值。
塌陷评价模块6，与相似实验模拟模块5连接，用于通过评价程序对岩溶塌陷进行评价。
如图3所示，本发明提供的地层倾角检测模块1检测方法如下：
S201，通过检测设备根据当前道分析点和各个预设扫描倾角，确定对应于所述各个预设扫描倾角的临近道分析点，根据所述当前道分析点和临近道分析点，分别获取当前道的演变波形数据和对应于所述各个预设扫描倾角的临近道的演变波形数据；
S202，沿层位获得演变波形数据，并采用基于线性变换的SSDR算法为演变波形数据降维；利用演变波形数据中的标签数据训练一个距离衡量矩阵；采用半监督Kmeans分类算法对演变波形数据进行分类，生成演变相图；
S203，根据所述当前道的演变波形数据、演变相图和对应于所述各个预设扫描倾角的临近道的演变波形数据，确定所述各个预设扫描倾角的相似能量谱；
S204，从各个相似能量谱中提取最大相似能量谱，从各个相似能量谱中提取与所述最大相似能量谱临近的若干相似能量谱，根据所述最大相似能量谱和所述若干临近的能量谱确定地层倾角插值曲线，并根据所述地层倾角插值曲线确定地层倾角；
其中，从各个相似能量谱中提取与所述最大相似能量谱相邻的两个相似能量谱；计算所述地层倾角插值曲线的一阶导数；将所述一阶导数等于零时所对应的倾角作为地层倾角；
根据如下公式计算所述地层倾角：
其中，x表示地层倾角，y1表示最大相似能量谱，x1表示最大相似能量谱对应的预设扫描倾角，y0和y2表示与所述最大相似能量谱相邻的两个相似能量谱。
本发明提供的在当前道中，获取以所述当前道分析点为中心的预设深度窗口的演变波形数据，作为所述当前道的演变波形数据；在各个临近道中，获取以相应临近道分析点为中心的预设深度窗口的演变波形数据，作为所述临近道的演变波形数据。
如图4所示，本发明提供的塌陷评价模块6评价方法如下：
S301，通过探地雷达法、瞬变电磁法、跨孔高密度电法和地表高密度电法岩溶进行探测；获取岩溶基本数据以及岩溶换算数据；
S302，建立沙漏型溶洞-地层三维数值模型，引入所获取的数据模拟计算岩溶上覆砂土漏失塌陷过程；
S303，提取岩溶塌陷后的底层沉降剖面图，根据熔岩塌陷对周围底层及建筑构筑物的影响确定岩溶沉降等级。
以上所述仅是对本发明的较佳实施例而已，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所做的任何简单修改，等同变化与修饰，均属于本发明技术方案的范围内。
Claims (10)
1.一种用于模拟岩溶地区地层的相似试验装置，其特征在于，所述用于模拟岩溶地区地层的相似试验装置包括：
地层倾角检测模块，与岩溶地层参数采集模块连接，用于检测岩溶地区地层倾角；
岩溶地层参数采集模块，与地层倾角检测模块、参数导入模块连接，用于采集岩溶地区地层物理参数、力学参数；
参数导入模块，与岩溶地层参数采集模块、模型构建模块连接，用于将采集的岩溶地层参数导入到模型构建模块；
模型构建模块，与参数导入模块、相似实验模拟模块连接，用于通过模型构建程序根据导入的参数构建岩溶地层模型；
塌陷评价模块，与相似实验模拟模块连接，用于通过评价程序对岩溶塌陷进行评价；
相似实验模拟模块，与模型构建模块、塌陷评价模块连接，用于通过模拟程序对模拟岩溶地区地层的相似试验进行模拟；所述模拟方法包括：
选择岩溶发育且演变资料上BSR特征明显的演变剖面；
确定演变剖面上岩溶地区地层底部界面的双程走时Tsb、观测BSR的双程走时Tbsr和岩溶地区地层顶部边界的双程走时Tsalt；
根据分析区天然气岩溶的气体组分信息，选择合适的岩溶相平衡曲线；
演变剖面上BSR发育深度处的强度值为静水强度；
BSR的深度值通过下面公式求得：
Hbsr＝Vsw×Tbsr/2
式中，Vsw为海水速度，为1500m/s；Tbsr为BSR位置的双程走时，单位是s，通过演变剖面读取；
确定演变剖面上岩溶地区地层底部各点的裂变值T；
通过二维稳态热传导方程，计算演变剖面的裂变场分布：
将岩溶稳定底界移向低的裂变条件；
通过设置不同的裂变梯度，模拟对应的BSR深度，将模拟的BSR位置与演变剖面上的观测BSR进行对比。
2.如权利要求1所述的用于模拟岩溶地区地层的相似试验装置，其特征在于，纯岩溶的相平衡公式如下：
In(P)＝a+bT+cT2+dT3+fT4+gT5
式中，P和T分别为岩溶的稳定强度条件和稳定裂变条件，a、b、c、d、f、g为经验常数，分别为a＝-1.94138504464560×105,b＝3.31018213397926×103,c＝-2.25540264493806×101,d＝7.67559117787059×10-2,f＝-1.30465829788791×10-4,g＝8.8606531668757×10-8。
3.如权利要求1所述的用于模拟岩溶地区地层的相似试验装置，其特征在于，静水强度通过以下公式求取：
Pbsr＝ρsw g Hbsr
式中，ρsw为海水密度，为1028kg/m3，g是重力加速度，为9.81m/s2，Hbsr是BSR的深度值。
4.如权利要求1所述的用于模拟岩溶地区地层的相似试验装置，其特征在于，计算演变剖面的裂变场分布：
式中，T为裂变℃，x为横向距离km，z为垂向距离km，kx为横向导热率W m-1K-1，kz为垂向导热率W m-1K-1；一般认为沉积物各向同性且均质性较强，因此kx＝kz，但沉积物和盐的热导率差别较大，分别设置为2.5W m-1K-1和5.9W m-1K-1。
5.如权利要求1所述的用于模拟岩溶地区地层的相似试验装置，其特征在于，岩溶稳定底界移向低的裂变条件：
式中，m为裂变对岩溶稳定条件的影响参数，计算公式如下：
式中，Sw为裂变值。
6.一种如权利要求1～5任意一项所述的用于模拟岩溶地区地层的相似试验装置的制作的方法，其特征在于，所述用于模拟岩溶地区地层的相似试验装置的制作的方法包括以下步骤：
步骤一，通过地层倾角检测模块检测岩溶地区地层倾角；通过岩溶地层参数采集模块采集岩溶地区地层物理参数、力学参数；
步骤二，通过参数导入模块将采集的岩溶地层参数导入到模型构建模块；通过模型构建模块利用模型构建程序根据导入的参数构建岩溶地层模型；
步骤三，通过相似实验模拟模块利用模拟程序对模拟岩溶地区地层的相似试验进行模拟；
步骤四，通过塌陷评价模块利用评价程序对岩溶塌陷进行评价。
7.如权利要求6所述用于模拟岩溶地区地层的相似试验装置，其特征在于，所述地层倾角检测模块检测方法如下：
1)通过检测设备根据当前道分析点和各个预设扫描倾角，确定对应于所述各个预设扫描倾角的临近道分析点，根据所述当前道分析点和临近道分析点，分别获取当前道的演变波形数据和对应于所述各个预设扫描倾角的临近道的演变波形数据；
2)沿层位获得演变波形数据，并采用基于线性变换的SSDR算法为演变波形数据降维；利用演变波形数据中的标签数据训练一个距离衡量矩阵；采用半监督Kmeans分类算法对演变波形数据进行分类，生成演变相图；
3)根据所述当前道的演变波形数据、演变相图和对应于所述各个预设扫描倾角的临近道的演变波形数据，确定所述各个预设扫描倾角的相似能量谱；
4)从各个相似能量谱中提取最大相似能量谱，从各个相似能量谱中提取与所述最大相似能量谱临近的若干相似能量谱，根据所述最大相似能量谱和所述若干临近的能量谱确定地层倾角插值曲线，并根据所述地层倾角插值曲线确定地层倾角；
其中，从各个相似能量谱中提取与所述最大相似能量谱相邻的两个相似能量谱；计算所述地层倾角插值曲线的一阶导数；将所述一阶导数等于零时所对应的倾角作为地层倾角；
根据如下公式计算所述地层倾角：
其中，x表示地层倾角，y1表示最大相似能量谱，x1表示最大相似能量谱对应的预设扫描倾角，y0和y2表示与所述最大相似能量谱相邻的两个相似能量谱。
8.如权利要求7所述用于模拟岩溶地区地层的相似试验装置，其特征在于，在当前道中，获取以所述当前道分析点为中心的预设深度窗口的演变波形数据，作为所述当前道的演变波形数据；在各个临近道中，获取以相应临近道分析点为中心的预设深度窗口的演变波形数据，作为所述临近道的演变波形数据。
9.如权利要求6所述用于模拟岩溶地区地层的相似试验装置，其特征在于，所述塌陷评价模块评价方法包括：
(1)通过探地雷达法、瞬变电磁法、跨孔高密度电法和地表高密度电法岩溶进行探测；获取岩溶基本数据以及岩溶换算数据；
(2)建立沙漏型溶洞-地层三维数值模型，引入所获取的数据模拟计算岩溶上覆砂土漏失塌陷过程。
10.如权利要求9所述用于模拟岩溶地区地层的相似试验装置，其特征在于，步骤(2)后，还进行：提取岩溶塌陷后的底层沉降剖面图，根据熔岩塌陷对周围底层及建筑构筑物的影响确定岩溶沉降等级。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210227771.9A CN114692391B (zh) | 2022-03-08 | 2022-03-08 | 用于模拟岩溶地区地层的相似试验装置及其使用方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210227771.9A CN114692391B (zh) | 2022-03-08 | 2022-03-08 | 用于模拟岩溶地区地层的相似试验装置及其使用方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114692391A true CN114692391A (zh) | 2022-07-01 |
| CN114692391B CN114692391B (zh) | 2024-06-25 |
Family
ID=82137328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210227771.9A Active CN114692391B (zh) | 2022-03-08 | 2022-03-08 | 用于模拟岩溶地区地层的相似试验装置及其使用方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114692391B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115392137A (zh) * | 2022-10-27 | 2022-11-25 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | 一种基于岩溶塌陷水土耦合作用的三维模拟系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112329103A (zh) * | 2020-11-04 | 2021-02-05 | 西南交通大学 | 一种岩溶上覆砂层塌陷对地层扰动的评价方法 |
| CN113075381A (zh) * | 2021-03-15 | 2021-07-06 | 华南理工大学 | 一种模拟岩溶地区地铁振动引发岩溶覆土塌陷的试验系统 |
| CN113311136A (zh) * | 2021-05-20 | 2021-08-27 | 柳州铁道职业技术学院 | 用于模拟岩溶地区地层的相似试验装置及其制作的方法 |
| CN215768563U (zh) * | 2021-05-20 | 2022-02-08 | 柳州铁道职业技术学院 | 一种用于模拟岩溶地区地层的相似试验装置 |
-
2022
- 2022-03-08 CN CN202210227771.9A patent/CN114692391B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112329103A (zh) * | 2020-11-04 | 2021-02-05 | 西南交通大学 | 一种岩溶上覆砂层塌陷对地层扰动的评价方法 |
| CN113075381A (zh) * | 2021-03-15 | 2021-07-06 | 华南理工大学 | 一种模拟岩溶地区地铁振动引发岩溶覆土塌陷的试验系统 |
| CN113311136A (zh) * | 2021-05-20 | 2021-08-27 | 柳州铁道职业技术学院 | 用于模拟岩溶地区地层的相似试验装置及其制作的方法 |
| CN215768563U (zh) * | 2021-05-20 | 2022-02-08 | 柳州铁道职业技术学院 | 一种用于模拟岩溶地区地层的相似试验装置 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115392137A (zh) * | 2022-10-27 | 2022-11-25 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | 一种基于岩溶塌陷水土耦合作用的三维模拟系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114692391B (zh) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110795878B (zh) | 一种隧道涌水量预测方法 | |
| Eiken et al. | Gravimetric monitoring of gas production from the Troll field | |
| CN100429531C (zh) | 目标最小化的三维电磁快速反演方法 | |
| US20150066460A1 (en) | Stratigraphic function | |
| Liu et al. | Improved inverse distance weighting method application considering spatial autocorrelation in 3D geological modeling | |
| CN112465191B (zh) | 隧道突涌水灾害预测方法、装置、电子设备及存储介质 | |
| Hare et al. | The 4D microgravity method for waterflood surveillance: Part IV—Modeling and interpretation of early epoch 4D gravity surveys at Prudhoe Bay, Alaska | |
| Lepillier et al. | From outcrop scanlines to discrete fracture networks, an integrative workflow | |
| Shean et al. | GPS-derived estimates of surface mass balance and ocean-induced basal melt for Pine Island Glacier ice shelf, Antarctica | |
| Su et al. | Experimental investigation of hydrate accumulation distribution in gas seeping system using a large scale three-dimensional simulation device | |
| CN114692391B (zh) | 用于模拟岩溶地区地层的相似试验装置及其使用方法 | |
| CN109145467B (zh) | 一种适用于台风区域的沙波运移预测方法 | |
| Berrino et al. | Repeated absolute gravity measurements on a dense network at Campi Flegrei–a reliable tool for volcano monitoring | |
| CN107526108B (zh) | 碳酸盐岩缝洞体积校正方法 | |
| CN105628576B (zh) | 灰质白云岩孔隙度测量方法及装置 | |
| CN111399037B (zh) | 高速顶界面提取的方法和装置 | |
| CN116398118A (zh) | 一种基于水平井随钻跟踪快速成像处理方法 | |
| CN116628990A (zh) | 一种模拟污染物在土中运移的模型试验方法 | |
| Perozzi et al. | Quantitative uncertainty analysis of gravity disturbance. The case of the Geneva Basin (Switzerland) | |
| CN106846481B (zh) | 一种地质剖面图的生成方法 | |
| CN110989034B (zh) | 一种回归-分形插值法反演测井横波时差方法 | |
| CN110687591B (zh) | 基于先验数据的波形匹配确定煤层及围岩物性参数的方法 | |
| CN112649880A (zh) | 砾岩下伏构造变速成图方法及装置 | |
| Verscheure et al. | History matching of a stochastic multifractal subesismic fault model | |
| Garner et al. | The Quest for innovative technology solutions for in-situ development of challenging oil sands reservoirs in Alberta |
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 |