CN107480349B - 基于三维地质模型和反分析的高铁沿线地面沉降预测方法 - Google Patents
基于三维地质模型和反分析的高铁沿线地面沉降预测方法 Download PDFInfo
- Publication number
- CN107480349B CN107480349B CN201710619280.8A CN201710619280A CN107480349B CN 107480349 B CN107480349 B CN 107480349B CN 201710619280 A CN201710619280 A CN 201710619280A CN 107480349 B CN107480349 B CN 107480349B
- Authority
- CN
- China
- Prior art keywords
- ground settlement
- water level
- monitoring
- ground
- inverse analysis
- 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 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 64
- 238000012544 monitoring process Methods 0.000 claims abstract description 59
- 239000002689 soil Substances 0.000 claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 10
- 238000011160 research Methods 0.000 claims abstract description 7
- 238000007906 compression Methods 0.000 claims description 40
- 230000006835 compression Effects 0.000 claims description 39
- 239000003673 groundwater Substances 0.000 claims description 15
- 238000007596 consolidation process Methods 0.000 claims description 14
- 230000035699 permeability Effects 0.000 claims description 8
- 238000011161 development Methods 0.000 claims description 6
- 238000004062 sedimentation Methods 0.000 claims description 5
- 230000002123 temporal effect Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 description 12
- 239000004927 clay Substances 0.000 description 6
- YSIITVVESCNIPR-UHFFFAOYSA-N Troxipide Chemical compound COC1=C(OC)C(OC)=CC(C(=O)NC2CNCCC2)=C1 YSIITVVESCNIPR-UHFFFAOYSA-N 0.000 description 4
- 229960001341 troxipide Drugs 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 239000012528 membrane Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- XHEFDIBZLJXQHF-UHFFFAOYSA-N fisetin Chemical compound C=1C(O)=CC=C(C(C=2O)=O)C=1OC=2C1=CC=C(O)C(O)=C1 XHEFDIBZLJXQHF-UHFFFAOYSA-N 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 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/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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于三维地质模型和反分析的高铁沿线地面沉降预测方法,该方法包括以下步骤:S1、获取研究区域内地质勘察数据、地下水位和地面沉降监测数据;S2、区域三维地质模型构建;S3、区域地下水位变化计算;S4、区域地面沉降模型构建;S5、地层力学参数反分析计算;S6、区域地面沉降量预测计算。该预测方法利用高速铁路工程地质勘察数据和地下水位、地面沉降监测数据,结合反分析构建能考虑土层厚度和土体力学参数空间不均匀分布的三维地面沉降预测模型,具有较高的预测精度,可广泛用于高速铁路沿线水位变化引发的地面沉降量的预测。
Description
技术领域
本发明涉及工程地质勘察领域,特别是涉及一种基于三维地质模型和反分析的高铁沿线地面沉降预测方法。
背景技术
目前京沪高速铁路、京津城际等10余条高速铁路经过严重的地面沉降区,地面沉降将影响高铁结构稳定和线路平顺性,严重威胁高速铁路的快速、平稳、安全运营,因此需要对高速铁路沿线地面沉降的发展趋势进行准确预测,以采取相应的防治措施。
当前主要采用统计法估算地面沉降的发展趋势或采用商业数值软件计算局部范围的地面沉降变化,其中统计法无法考虑三维地质力学条件,难以准确预测沉降长期发展;商业数值软件通常适用于模拟局部均质地层下的地面沉降,不能自动对土层参数进行反分析,预测结果无法达到工程应用的精度要求。由于高速铁路是线状工程,沿线穿越多种复杂的工程地质条件,地层分布和土体力学参数具有显著的不均匀性,因此已有地面沉降预测手段难以满足高速铁路沿线的地面沉降预测需求。
发明内容
针对高速铁路沿线开采地下水引发的地面沉降问题,为解决现有地面沉降预测方法的局限性,本发明提供一种能够快速建立考虑地层厚度和力学参数空间不均匀分布的地面沉降预测模型、准确预测计算高速铁路沿线地面沉降量随时间的发展状况的基于三维地质模型和反分析的高铁沿线地面沉降预测方法。
为此,本发明的技术方案如下:
一种基于三维地质模型和反分析的高铁沿线地面沉降预测方法,根据高铁工程勘察数据、地下水位和地面沉降监测数据并结合反分析建立高铁沿线三维地面沉降预测模型,根据预期地下水位变化状况预测计算高速铁路沿线地面沉降量的发展状况,其具体步骤如下:
S1、获取研究区域内地质勘察数据、地下水位和地面沉降监测数据:获取区域内高速铁路工程勘探点地层数据,获取区域内地下水监测点处的水位变化数据,获取区域内地面沉降量的变化数据;
S2、区域三维地质模型构建:根据步骤S1中获取的数据识别区域地层,对区域进行离散网格划分,离散化计算区域内各地层厚度值的空间分布,形成区域三维地质模型;
S3、区域地下水位场计算:根据步骤S1中获取的水位变化数据,按照步骤S2中的离散网格,插值计算区域内各含水层在各监测时刻的水位场;
S4、区域地面沉降预测模型构建:根据步骤S2中构建的三维地质模型,按离散网格识别含水层压缩系统和弱透水层固结压缩系统,根据步骤S3中的地下水位场确定压缩系统的边界条件,形成区域地面沉降计算模型;
S5、地层力学参数反分析计算:根据步骤S4中构建的区域地面沉降预测模型,确定沉降监测点处的地面沉降量s(x,t)的计算方法,根据步骤S1中的地面沉降监测数据,建立反分析目标函数:
式中:x为土层力学参数向量,ti为第i个监测时刻,s(x,ti)为ti时的沉降计算值,si为ti时的沉降监测值,在限定条件下对土层力学参数进行反分析计算,利用反分析得到的力学参数更新区域地面沉降预测模型;
S6、区域地面沉降量预测计算:设定地下水位随时间的预期变化,按步骤S3计算区域内地下水位场预期变化数据,更新步骤S4中地面沉降预测模型的边界条件,利用更新后的地面沉降预测模型计算区域内地面沉降量随时间的变化。
在上述方法中:
步骤S1中所述的高速铁路工程勘探点地层数据包括:勘探点的位置坐标、高程,各层土的层厚、土体类别、土体的压缩系数、渗透系数;所述的地下水位监测数据包括:地下水位监测点的位置坐标、水位监测地层、水位各监测时刻及相应的水位深度;所述的地面沉降监测数据包括:地面沉降监测点的位置坐标、各监测时刻及相应的地面沉降量。
步骤S2所述的识别区域地层包括:整理统计钻孔地层数据、根据土性判断土层是否为含水层或弱透水层、确定区域标准地层分布;所述的离散网格划分是根据所需模型计算精度,并结合已知的勘探点、监测点位置,将平面区域划分为四边形网格;所述的离散化计算区域内各地层厚度值的空间分布是根据已知勘探点处各地层的厚度按克里格法插值计算各网格点处各地层的厚度。
步骤S3所述插值计算地下水位场为采用径向基函数法、反距离加权法对已知监测点处地下水位按离散网格点进行插值。
步骤S4所述的含水层压缩系统包括所有发生水位变化的含水层,弱透水层固结压缩系统包括所有被含水层分隔的弱透水层组,根据系统边界处地下水位变化,按弹性压缩计算含水层系统的压缩量,按固结理论计算所述弱透水层固结压缩系统的固结压缩量,将各系统压缩量累加后获得地面沉降量。
步骤S5所述的地层力学参数包括:含水层的压缩系数、弱透水层的压缩系数和渗透系数,步骤S5所述的在限定条件下对地层力学参数进行反分析计算包括:设定地层力学参数范围为0.1~10倍的初始值、按Levenberg-Marquardt法对目标函数进行迭代计算直至满足收敛条件;步骤S5所述的更新区域地面沉降模型为采用克里格法将反分析得到的力学参数在离散网格点处进行插值,更新初始地层力学参数。
本发明具有的优点和积极效果是:
通过综合利用高铁工程勘察数据、地下水位和地面沉降监测数据,能够结合反分析快速建立考虑地层厚度和力学参数空间不均匀分布的高速铁路沿线三维地面沉降预测模型,针对预期地下水位变化,能够准确预测高铁沿线地面沉降变化,为合理设计高铁沉降防治措施提供可靠的数据基础。
附图说明
图1为本发明的基于三维地质模型和反分析的高铁沿线地面沉降预测方法的流程图;
图2为网格点处土层系统识别示意图。
图中:
1、含水层压缩系统; 2、第一弱透水层组; 3、第二弱透水层组;
4、弱透水层固结压缩系统。
具体实施方式
下面结合附图对本发明的基于三维地质模型和反分析的高铁沿线地面沉降预测方法做进一步说明。
参见图1,本发明的基于三维地质模型和反分析的高铁沿线地面沉降预测方法主要包括如下步骤:
S1、获取研究区域内地质勘察数据、地下水位和地面沉降监测数据;S2、区域三维地质模型构建;S3、区域地下水位变化计算;S4、区域地面沉降模型构建;S5、地层力学参数反分析计算;S6、区域地面沉降量预测计算。具体如下:
S1、获取研究区域内地质勘察数据、地下水位和地面沉降监测数据:
获取高速铁路工程勘探地层数据,包括勘探点的位置坐标、高程,各层土的层厚、土体类别、土体的压缩系数、渗透系数;获取地下水位监测数据,包括地下水位监测点的位置坐标、水位监测层位、水位各监测时刻及相应的水位深度;获取地面沉降监测数据,包括地面沉降监测点的位置坐标、各监测时刻及相应的地面沉降量。
S2、区域三维地质模型构建:
首先根据S1中的勘察数据识别区域地层:
①整理钻孔地层数据,对室内试验得到的土体力学参数进行统计,确定地层的标准力学参数值;
②根据土性判断土层是否为含水层或弱透水层,当土体渗透系数小于10-6m/s时土层为弱透水层,当土体渗透系数大于10-6m/s时土层为含水层;
③根据地质成因、地质时代分析钻孔地层从上至下的分布,确定区域标准地层分布;
其次对研究区域进行离散网格划分,根据模型计算精度需求确定网格标准长度,对研究区域进行初步划分,利用已知的勘探点、监测点的位置坐标进一步对网格进行加密,获得最终的离散网格。
然后离散化计算区域内各地层厚度值的空间分布:
①计算网格节点与步骤S1中各勘探点的距离值ri,计算各勘探点之间的距离rij;
②根据区域标准地层分布,选取当前待离散地层为标准地层的最顶部土层,根据勘探点数据确定各勘探点处土层厚度为bi;
④将距离值rij代入拟合后的变差公式计算变差值vij,将距离值ri代入拟合后的变差公式中计算变差值vi;
⑦重复④~⑥计算所有网格节点处的土层厚度;
⑧选取待离散地层为下一层标准地层,重新确定各勘探点处土层厚度bi,重复③~⑦直至标准地层最底部土层。
S3、区域地下水位场计算:
①计算步骤S2中网格节点与步骤S1中各水位监测点的距离值ri,计算各水位监测点之间的距离rij;
⑤选取下一时刻的地下水位监测值pi,重复③~④直至最终的监测时刻。
⑥重复①~⑤计算所有网格节点处的地下水位值。
S4、区域地面沉降预测模型构建:
①根据S2中构建的区域三维地质模型确定网格点处地层分布如图2所示,从上至下为细砂层、黏土层、粉质黏土层、中砂层、粉质黏土层、中砂层;
②从上至下对图2中土层进行识别后,从上至下确定细砂层、第一个中砂层、第二个中砂层为含水层,共同构成含水层压缩系统1,黏土层和第一个粉质黏土层构成第一弱透水层组2,第二个粉质黏土层单独构成第二弱透水层组3,由第一、第二弱透水层组2和3共同构成弱透水层固结压缩系统4;
③根据S3中的地下水位场确定三个含水层(即细砂层、第一个中砂层、第二个中砂层)中的水位变化分别为h1(t)、h2(t)、h3(t),其中t为监测时刻;
⑤按固结理论计算弱透水层系统中各个弱透水层组的固结压缩量,建立控制方程为:
式中:
z为竖向坐标,以土层顶部为原点,
对控制方程和边界条件进行Laplace变换后可得:
利用传递矩阵法可得:
式中:
L-1表示Laplace逆变换,可采用Crump法。
⑥将所有土层的压缩量进行累加后得出网格点处的沉降量计算式为:
式中:
se(t)为含水层压缩系统的沉降值,
sc(t)为弱透水层固结压缩系统的沉降值,
ne为含水层的个数,
nc为弱透水层组的个数,
nl,j为第j个弱透水层组中弱透水层的个数。
S5、地层力学参数反分析计算:
①根据步骤S4中构建的区域地面沉降预测模型,确定沉降监测点处的地面沉降量s(x,t)的计算方法,根据步骤S1中的地面沉降监测数据,建立反分析目标函数:
式中:x为土层力学参数向量,ti为第i个监测时刻,s(x,ti)为ti时的沉降计算值,si为ti时的沉降监测值;
②设定地层力学参数范围为0.1~10倍的初始值、按Levenberg-Marquardt法对目标函数进行迭代计算直至满足收敛条件;
③采用克里格法将反分析得到的力学参数在离散网格点处进行插值,更新地面沉降预测模型中的地层力学参数。
S6、区域地面沉降量预测计算:
设定地下水位随时间的预期变化,按步骤S3计算区域内地下水位场预期变化数据,更新步骤S4中地面沉降预测模型的边界条件,利用更新后的地面沉降预测模型计算区域内地面沉降量随时间的变化。
Claims (6)
1.一种基于三维地质模型和反分析的高铁沿线地面沉降预测方法,根据高铁工程勘察数据、地下水位和地面沉降监测数据并结合反分析建立高铁沿线三维地面沉降预测模型,根据预期地下水位变化状况预测计算高速铁路沿线地面沉降量的发展状况,该方法具有下述步骤:
S1、获取研究区域内地质勘察数据、地下水位和地面沉降监测数据:获取区域内高速铁路工程勘探点地层数据,获取区域内地下水监测点处的水位变化数据,获取区域内地面沉降量的变化数据;
S2、区域三维地质模型构建:根据步骤S1中的地层资料识别区域地层,离散化计算区域内各地层厚度值的空间分布,形成区域三维地质模型;
S3、区域地下水位场计算:根据步骤S1中的地下水位监测数据,按照步骤S2中的离散网格,插值计算区域内各含水层在各监测时刻的水位场;
S4、区域地面沉降预测模型构建:根据步骤S2中构建的三维地质模型,按离散网格识别含水层压缩系统和弱透水层固结压缩系统,根据步骤S3中的地下水位场确定压缩系统的边界条件,形成区域地面沉降预测模型;
S5、地层力学参数反分析计算:根据步骤S4中构建的区域地面沉降预测模型,确定沉降监测点处的地面沉降量s(x,t)的计算方法,根据步骤S1中的地面沉降监测数据,建立反分析目标函数在限定条件下对土层力学参数进行反分析计算,利用反分析得到的力学参数更新区域地面沉降预测模型,式中:x为土层力学参数向量,ti为第i个监测时刻,s(x,ti)为ti时的沉降计算值,si为ti时的沉降监测值;
S6、区域地面沉降量预测计算:设定地下水位随时间的预期变化,按步骤S3计算区域内地下水位场预期变化数据,更新步骤S4中地面沉降预测模型的边界条件,利用更新后的地面沉降预测模型计算区域内地面沉降量随时间的变化。
2.根据权利要求1所述的基于三维地质模型和反分析的高铁沿线地面沉降预测方法,其特征在于:步骤S1所述的高速铁路工程勘探点地层数据,包括勘探点的位置坐标、高程,各层土的层厚、土体类别、土体的压缩系数、渗透系数;所述的地下水位监测数据,包括地下水位监测点的位置坐标、水位监测地层、水位各监测时刻及相应的水位深度;所述的地面沉降监测数据,包括地面沉降监测点的位置坐标、各监测时刻及相应的地面沉降量。
3.根据权利要求1所述的基于三维地质模型和反分析的高铁沿线地面沉降预测方法,其特征在于:步骤S2所述的识别区域地层包括整理统计钻孔地层数据、根据土性判断土层是否为含水层或弱透水层、确定区域标准地层分布;所述的离散化计算区域内各地层厚度值的空间分布包括划分区域网格、根据已知勘探点处的地层厚度按克里格法插值计算各网格点处的土层厚度。
4.根据权利要求1所述的基于三维地质模型和反分析的高铁沿线地面沉降预测方法,其特征在于:步骤S3中,插值计算区域内各含水层在各监测时刻的水位场时采用径向基函数法、反距离加权法对已知监测点处地下水位按离散网格点进行插值。
5.根据权利要求1所述的基于三维地质模型和反分析的高铁沿线地面沉降预测方法,其特征在于:步骤S4所述的含水层压缩系统包括所有发生水位变化的含水层,弱透水层固结压缩系统包括所有被含水层分隔的弱透水层组,根据系统边界处地下水位变化,按弹性压缩计算含水层系统的压缩量,按固结理论计算所述弱透水层组的固结压缩量,将各系统压缩量累加后获得地面沉降量。
6.根据权利要求1所述的基于三维地质模型和反分析的高铁沿线地面沉降预测方法,其特征在于:步骤S5所述的地层力学参数包括含水层的压缩系数、弱透水层的压缩系数和渗透系数,步骤S5所述的在限定条件下对地层力学参数进行反分析计算包括:设定地层力学参数范围为0.1~10倍的初始值、按Levenberg-Marquardt法进行对目标函数进行迭代计算直至满足收敛条件,步骤S5所述的更新区域地面沉降预测模型为采用克里格法将反分析得到的力学参数在离散网格点处进行插值,更新初始地层力学参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710619280.8A CN107480349B (zh) | 2017-07-26 | 2017-07-26 | 基于三维地质模型和反分析的高铁沿线地面沉降预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710619280.8A CN107480349B (zh) | 2017-07-26 | 2017-07-26 | 基于三维地质模型和反分析的高铁沿线地面沉降预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107480349A CN107480349A (zh) | 2017-12-15 |
CN107480349B true CN107480349B (zh) | 2020-07-24 |
Family
ID=60597999
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710619280.8A Active CN107480349B (zh) | 2017-07-26 | 2017-07-26 | 基于三维地质模型和反分析的高铁沿线地面沉降预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107480349B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108876918B (zh) * | 2018-07-02 | 2022-07-22 | 长安大学 | 一种高填方黄土体沉降变形预测方法 |
CN109098753B (zh) * | 2018-08-09 | 2019-07-26 | 中国矿业大学 | 一种浅埋煤层开采潜水漏失致灾程度的划分方法 |
CN108827233B (zh) * | 2018-09-17 | 2019-12-10 | 中国地质大学(北京) | 一种两层采空区地面沉降的预测方法 |
CN109446675B (zh) * | 2018-11-01 | 2023-04-07 | 河北昕佳工程勘查设计有限公司 | 一种地热资源量计算修正方法及沉降预警系统 |
CN109165480B (zh) * | 2018-11-01 | 2023-08-11 | 河北昕佳工程勘查设计有限公司 | 一种基于粘性土应力路径预测矿坑涌水量的方法及系统 |
CN110398227A (zh) * | 2019-07-29 | 2019-11-01 | 中国电建集团铁路建设有限公司 | 一种可溶岩场地近接地铁车站异常差异沉降监测方法 |
CN110851991B (zh) * | 2019-11-18 | 2023-07-14 | 核工业二〇八大队 | 一种地下水流数值模拟方法 |
CN111046586B (zh) * | 2019-12-29 | 2023-04-07 | 河南省地质矿产勘查开发局第二地质环境调查院 | 开采深部松散层地热引发地面沉降量的预测方法 |
CN113128106B (zh) * | 2021-04-06 | 2023-03-21 | 汕头大学 | 一种岩溶地层盾构施工引起地表沉降的确定方法 |
CN113378406A (zh) * | 2021-06-30 | 2021-09-10 | 江苏旭辰交通科技发展有限公司 | 多层土工格室加固的非粘性土圆形基础沉降预测方法 |
CN115540816A (zh) * | 2022-10-12 | 2022-12-30 | 广东电网有限责任公司 | 变电站设备的沉降监测方法、装置、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014020167A (ja) * | 2012-07-23 | 2014-02-03 | Kokuka Corp | 建造物周囲の地盤構築工法並びに該工法に用いられる浸透側溝 |
CN105320830A (zh) * | 2014-07-28 | 2016-02-10 | 中国石油天然气集团公司 | 构建双溶腔盐穴储库地面沉降预测模型的方法和装置 |
KR101610343B1 (ko) * | 2015-08-17 | 2016-04-08 | (주)한국이엔지 | 터널굴착용 지하수 또는 토사유출방지장치 |
-
2017
- 2017-07-26 CN CN201710619280.8A patent/CN107480349B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014020167A (ja) * | 2012-07-23 | 2014-02-03 | Kokuka Corp | 建造物周囲の地盤構築工法並びに該工法に用いられる浸透側溝 |
CN105320830A (zh) * | 2014-07-28 | 2016-02-10 | 中国石油天然气集团公司 | 构建双溶腔盐穴储库地面沉降预测模型的方法和装置 |
KR101610343B1 (ko) * | 2015-08-17 | 2016-04-08 | (주)한국이엔지 | 터널굴착용 지하수 또는 토사유출방지장치 |
Non-Patent Citations (2)
Title |
---|
Strength and deformation of silt foundation effected by groundwater level fluctuation;Long Zhang等;《Applied mechanics and materials》;20121230;第204-208卷;第571-576页 * |
京沪高速铁路沿线地面沉降与地下水位变化关系探讨;李国和等;《水文地质工程地质》;20080630(第6期);第90-94页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107480349A (zh) | 2017-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107480349B (zh) | 基于三维地质模型和反分析的高铁沿线地面沉降预测方法 | |
KR101618713B1 (ko) | 지반정보에 기초한 3차원 공간 모델링 방법 | |
CN104730596B (zh) | 一种基于多尺度因素约束的离散裂缝建模方法 | |
KR101642951B1 (ko) | Gis 기반 실시간 지진피해 예측 방법 | |
CN106401574B (zh) | 一种钻前高温地热井地层压力的预测方法 | |
CN105242307B (zh) | 复杂碳酸盐储层地震孔隙度获取方法及装置 | |
CN103077558A (zh) | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 | |
GB2528384A (en) | Completion design based on logging while drilling (LWD) data | |
CN110056346B (zh) | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 | |
WO2014151440A1 (en) | Targeted survey design under uncertainty | |
CN104616353A (zh) | 储层随机地质模型的建模及优选方法 | |
CN105137482A (zh) | 一种沉积体古坡度的计算方法 | |
Burbery et al. | Study of connectivity of open framework gravel facies in the Canterbury Plains aquifer using smoke as a tracer | |
Mahmod et al. | Analysis of groundwater flow in arid areas with limited hydrogeological data using the Grey Model: a case study of the Nubian Sandstone, Kharga Oasis, Egypt | |
CN113552621A (zh) | 页岩气地应力确定方法和装置 | |
CN108563832A (zh) | 一种针对地浸采铀流场的岩相物性模拟方法 | |
US10401808B2 (en) | Methods and computing systems for processing and transforming collected data to improve drilling productivity | |
Albarrán-Ordás et al. | The D i models method: geological 3-D modeling of detrital systems consisting of varying grain fractions to predict the relative lithological variability for a multipurpose usability | |
CN110309597A (zh) | 基于构型界面的阻流带模型确定方法、装置及存储介质 | |
Frick et al. | The geothermal field below the city of Berlin, Germany: Results from structurally and parametrically improved 3D Models | |
CN115422789B (zh) | 一种考虑全过程优化断块油藏水驱采收率预测方法及系统 | |
Hamdulla et al. | Building geological model for tertiary reservoir of exploration Ismail oil field, North Iraq | |
Witty et al. | Application of Geostatistical Sequential Simulation Methods for Probabilistic 3D Subsoil Modeling and Uncertainty Quantification Concept and Examples | |
Rogiers et al. | Groundwater model parameter identification using a combination of cone-penetration tests and borehole data | |
Carpio | HIGH-RESOLUTION FACIES INTERPRETATION OF QUATERNARY SEDIMENTARY SERIES FOR HYDROGEOLOGICAL MODELLING OF THE NYÍRSÉG-HAJDÚSÁG GROUNDWATER BODY, NE HUNGARY. |
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 |