CN111401716A - 一种基于叠加放水试验获取含水层连续参数的方法 - Google Patents
一种基于叠加放水试验获取含水层连续参数的方法 Download PDFInfo
- Publication number
- CN111401716A CN111401716A CN202010161020.2A CN202010161020A CN111401716A CN 111401716 A CN111401716 A CN 111401716A CN 202010161020 A CN202010161020 A CN 202010161020A CN 111401716 A CN111401716 A CN 111401716A
- Authority
- CN
- China
- Prior art keywords
- aquifer
- water
- water discharge
- discharge test
- parameters
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 167
- 238000012360 testing method Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 238000011160 research Methods 0.000 claims abstract description 20
- 238000004587 chromatography analysis Methods 0.000 claims abstract description 8
- 238000005315 distribution function Methods 0.000 claims abstract description 6
- 238000007599 discharging Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000005086 pumping Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 238000005553 drilling Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 230000002159 abnormal effect Effects 0.000 abstract description 9
- 238000005065 mining Methods 0.000 abstract description 8
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 12
- 238000009530 blood pressure measurement Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 239000003245 coal Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Agronomy & Crop Science (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及矿山水文地质监测技术领域，具体涉及一种基于叠加放水试验获取含水层连续参数的方法，步骤一：在研究区域实施叠加放水试验；步骤二：建立研究区域的数值模型并划分网格离散模型；步骤三：设定含水层参数的随机分布函数值；步骤四：选择适用的协方差函数描述含水层参数的变异程度；步骤五：利用放水层析法迭代计算含水层参数的连续分布。本发明的方法操作简单，适用性广，相对于传统分析方法，获取的水文地质参数更加详细具体，可准确刻画含水层的非均质性，精确定位富水异常区的位置及范围，有效降低矿井突水灾害风险，为矿区安全生产提供综合技术保障。
Description
技术领域
本发明涉及矿山水文地质监测技术领域，具体涉及一种基于叠加放水试验获取含水层连续参数的方法。
背景技术
公开该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不必然被视为承认或以任何形式暗示该信息构成已经成为本领域一般技术人员所公知的现有技术。
含水层水文地质参数的空间分布特征对于采矿工作影响巨大，是分析矿区含水层富水性强弱的基础。富水异常区的位置及范围，直接制约着井巷系统的布局，影响矿井安全。查明含水层参数的详细信息，可为井巷系统的设计提供依据，提前采取有效防治措施并封堵导水通道，保证采矿安全。
目前矿区水文地质调查的主要手段有抽水试验、放水试验、地球物理勘探等方法。大流量、大降深、大范围的放水试验是获取水文地质信息的重要手段，传统放水试验分析方法得到的含水层水文地质参数是含水层的分区等效参数，反映出的主要地质构造范围粗略，界限模糊。地球物理勘探获取的是地层的电导率、介电常数等物理信号，无法直接获得水文地质参数，数据解译时容易增加人为的不确定性。
发明内容
为了解决现有技术中存在的技术问题，本发明的目的是提供一种基于叠加放水试验获取含水层连续参数的方法，该方法操作简单，适用性广，获取的水文地质参数详细具体，可准确刻画含水层的非均质性，精确定位富水异常区的位置及范围，有效降低矿井突水灾害风险，为矿区安全生产提供综合技术保障。
为了解决以上技术问题，本发明一个或多个实施例提供的技术方案为：
一种基于叠加放水试验获取含水层连续参数的方法，包括如下步骤：
在研究区域实施叠加放水试验；
建立研究区域的数值模型并划分网格离散模型；
设定含水层参数的随机分布函数值；
选择适用的协方差函数描述含水层参数的变异程度；
利用放水层析法迭代计算含水层参数的连续分布。
本发明的以上一个或多个实施例的有益效果为：
通过一系列分阶段分区域叠加放水试验，设定随机分布的含水层参数，分析试验的水压数据，以获取含水层水文地质参数的连续分布，为判断富水异常区提供详细的参数信息。该方法操作简单，适用性广，获取的水文地质参数详细具体，可准确刻画含水层的非均质性，精确定位富水异常区的位置及范围，有效降低矿井突水灾害风险，为矿区安全生产提供综合技术保障。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1是本发明实施例的基于叠加放水试验获取含水层连续参数的方法流程图。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
一种基于叠加放水试验获取含水层连续参数的方法，包括如下步骤：
在研究区域实施叠加放水试验；
建立研究区域的数值模型并划分网格离散模型；
设定含水层参数的随机分布函数值；
选择适用的协方差函数描述含水层参数的变异程度；
利用放水层析法迭代计算含水层参数的连续分布。
此处的叠加放水试验，涉及的是群孔放水试验，叠加的含义体现在时间和空间两方面。在某一段时间，某一区域的放水孔实施放水，再接着的另一段时间，另一区域的放水孔开启放水，与此同时，前一段时间开启的放水孔仍在放水。这种放水方案设计在煤矿是普遍的，在这里重点是强调对试验数据的分析，试验是分阶段分区域实施的，代入不同阶段的水压数据，逐步反演出参数的分布情况，体现水力“层析”的内涵。
在一些实施例中，实施叠加放水试验时，潜在富水区的钻孔加倍布设。
进一步的，在潜在富水区的钻孔密度加密2-3倍。保证钻孔在使用过程中不塌孔、不堵塞，保障试验的正常进行和数据测量的质量。依据现场实际情况确定放水的次数，每次放水运用到的一个或多个放水孔及放水持续时间。
进一步的，在研究区域分阶段分区域依次实施放水，所有观测孔同步测量水压，每次试验的前期加密观测水压，中后期定时观测。
进一步的，叠加放水试验的次数大于或等于3次。以作验证之用。在每个阶段放水试验的前期，水压响应强烈，需加密观测；在试验进行的中后期，水压变动缓慢并趋于稳定，定时观测即可。
在一些实施例中，网格离散模型中的网格尺寸小于或等于相关长度(Correlationscale)的1/2。
相关长度反映的是某一参数在空间上的连续性，相关长度大表征参数各个位置的值变化相对较小，连续性好。
进一步的，设定研究区域为非均质饱和含水层，根据研究区域的规模及对反演精度的需求，确定离散网格的规格。
在一些实施例中，设定含水层参数的随机分布函数值时，根据已知的地质信息预估水文地质参数的均值，将参数的初始估计值设定为均质分布；
或，在均值基础上，为每个离散元的参数附加符合正态分布的扰动值，参数的初始估计值设定为非均质分布。
进一步的，地质信息包括地质填图、钻孔岩芯取样及抽水试验结果。
在一些实施例中，选择适用的协方差函数描述含水层参数的变异程度时，通过方差和相关长度描述含水层参数的变异程度。
进一步的，研究区域在无经验协方差函数时，根据含水层类型，方差和相关长度的赋值如下：对于孔隙含水层，方差设为0.5-1，相关长度设为20-50m；对于裂隙含水层，方差设为1-2，相关长度设为50-200m；对于岩溶含水层，方差设为1-2，相关长度设为100-500m。
在一些实施例中，利用放水层析法迭代计算含水层参数的连续分布的具体步骤如下：
将水压数据整理为相等时间间隔的数据集，并进行平滑降噪，提取数据的明显变化规律；
确定钻孔在反演区域中的相对位置，平滑处理过后的水压值按照时间顺序依次赋值于对应钻孔结点处，边界包括定水头边界和定流量边界，定水头边界保证含水层始终维持饱和；
设定试验开始时间为反演计算的初始时间，初始时间区域内部结点的水头值与定水头边界的水头值数值相等；
根据放水试验的顺序逐次反演并迭代计算，动态显示连续参数的刻画过程，迭代停止的条件包括两方面，一方面是限制迭代次数，另一方面是结果位于误差允许范围之内。
如图1所示，本发明的基于叠加放水试验获取含水层连续参数的方法，包括下述步骤：
步骤一：在研究区域实施叠加放水试验，潜在富水区钻孔加倍布设，试验次数不少于3次，在每次试验的前期加密观测水压，中后期定时观测。具体过程为：根据所研究含水层的空间位置及范围，均匀布设地面或地下钻孔，在潜在富水异常区钻孔密度加密2-3倍。保证钻孔在使用过程中不塌孔、不堵塞，保障试验的正常进行和数据测量的质量。依据现场实际情况确定放水的次数，每次放水运用到的一个或多个放水孔及放水持续时间。试验次数不少于3次，以作验证之用。在每个阶段放水试验的前期，水压响应强烈，需加密观测；在试验进行的中后期，水压变动缓慢并趋于稳定，定时观测即可。
步骤二：建立研究区域的数值模型，模型离散化后的网格尺寸不大于相关长度的1/2。具体过程为：设定研究区域为非均质饱和含水层，充分考虑边界对结果的影响，确定最终的反演区域。根据反演区域的规模及对反演精度的需求，模型离散化后的网格尺寸不大于相关长度的1/2。
步骤三：设定含水层参数的随机分布函数值。具体过程为：先验地质信息的来源包括地质填图、钻孔岩芯取样及抽水试验结果，根据已知的地质信息预估水文地质参数的均值，将参数的初始估计值设定为均质分布。也可选择在均值基础上，为每个离散元的参数附加符合正态分布的扰动值，参数的初始估计值设定为非均质分布。
步骤四：选择适用的协方差函数，通过方差和相关长度描述含水层参数的变异程度。具体过程为：研究区在无经验协方差函数的情况下，根据含水层类型，方差和相关长度的赋值参考如下：对于孔隙含水层，方差可设0.5-1，相关长度可设20-50m；对于裂隙含水层，方差可设1-2，相关长度可设50-200m；对于岩溶含水层，方差可设1-2，相关长度可设100-500m。
步骤五：利用放水层析法，迭代计算含水层参数的连续分布。具体过程为：依据现场实际情况确定时间步长，将测得的水压数据处理为相等时间间隔的数据集，并进行平滑降噪，提取出数据的明显变化规律。确定钻孔在反演区域中的相对位置。平滑处理过后的水压值按照时间顺序依次赋值于对应钻孔结点处。边界包括定水头边界和定流量边界。依据现场实际情况确定边界为定水头边界或零流量边界。在试验过程中，定水头边界保证含水层始终维持饱和。设定试验开始时间为反演计算的初始时间，且整个区域内部结点的水头值与定水头边界的水头值数值相等。根据放水试验的顺序逐次反演并迭代计算，动态显示连续参数的刻画过程。迭代停止的条件包括两方面，一方面是限制迭代次数，另一方面是结果位于误差允许范围之内。
以导水系数为例，初始估计结果为
式中：
为x0处的导水系数扰动值的估计值；H*为实际水压测量值；He为通过地下水流控制方程求得的水压模拟值；Nf为参数测量点的总数；Np为放水试验的总次数；Nh(k)为第k次放水试验中水压测量点的总数；Nt(k,j)为第k次放水试验中第j个水压测量点随时间变化的水压测量数据的总数；λ0i为协同克里格权重，即第i个参数测量点的参数值f*对x0处的参数估计值的权重；μ0kjl为实际水压测量值H*(k,xj,tl)对参数估计值的权重。第一次参数估计结果，通过
变换代入地下水流控制方程进行下一步迭代计算。
迭代计算公式为
式中：
为第r次迭代时，第k次放水试验中第j个水压测量点xj在tl时刻，实际水压测量值H*(k,xj,tl)和水压模拟值H(r)(k,xj,tl)的差值对x0处的参数估计值的权重。同样地，将参数估计值代入地下水流控制方程，计算在水压测量点不同时刻的水压模拟值，判断模拟值与实际测量值之间的差异，若差异没有达到指定精度，上述计算过程会重复迭代进行。
下面通过一个实施例来进一步说明本发明。
本实施例以某煤矿的S1210工作面井下叠加放水试验为例，步骤如下：
步骤一：S1210工作面井下施工钻孔数量较多，自切眼往北约2000m范围内工作面顶板有超过60个钻孔，故可利用这些井下钻孔进行采前大流量放水试验。根据放水前对井下所有钻孔涌水量的排查情况，共选取工作面顶板上30个水量较大的钻孔作为放水孔。分区域选取28个井下钻孔(其中工作面顶板钻孔17个，巷道外侧钻孔11个)和6个地面钻孔作为水压观测孔。放水试验共持续9天，放水次序为自南向北，具体采用叠加放水的方式，观测孔全程测量水压。试验正式开始时，开启第一阶段的6个放水孔，位于研究区南部区域，第一阶段放水持续2天；随后开启第二阶段的7个放水孔，此时共有13个钻孔联合放水，第二阶段放水持续2天；继而开启第三阶段的8个放水孔，此时共有21个钻孔联合放水，第三阶段放水持续3天；最后开启第四阶段的9个放水孔，位于研究区最北部，此时共有30个钻孔联合放水，第四阶段放水持续2天。放水孔流量、地面观测孔水压、井下观测孔水压同步观测。观测频率：每阶段放水开始后3小时内为1次/5min，之后为1次/30min。
步骤二：建立研究区域的二维数值模型为1800m×3200m的矩形区域，将反演区域划分为90×80个单元格，每个单元格的规格为20m×40m。
步骤三：目标含水层岩性为直罗组砂岩，设定导水系数对数表示的均值为1.3，储水系数对数表示的均值为-5.0。
步骤四：对于直罗组砂岩含水层，设定导水系数对数表示的方差为1.0，储水系数对数表示的方差为0.5，x和y方向的相关长度均为200m。
步骤五：依据现场实际情况确定时间步长为1h，将测得的水压数据处理为相等时间间隔1h的数据集，并进行平滑降噪，提取出数据的明显变化规律。确定钻孔在反演区域中的相对位置。平滑处理过后的水压数据按照时间顺序依次赋值于对应钻孔结点处。设置四周为1000m的定水头边界，初始水头为1000m。由反演结果可推知，S1210工作面南北富水性差异大，从涌水区往北，导水系数及储水系数整体上逐渐减小。南部导水系数和储水系数均较大，属于富水异常区，可能位于古河流冲刷带，裂隙发育，富水性强。
本次使用的方法是一种基于传统放水试验的改进分析方法，通过一系列分阶段分区域叠加放水试验，设定随机分布的含水层参数，分析试验的水压数据，以获取含水层水文地质参数的连续分布，为判断富水异常区提供详细的参数信息。该方法操作简单，适用性广，获取的水文地质参数详细具体，可准确刻画含水层的非均质性，精确定位富水异常区的位置及范围，有效降低矿井突水灾害风险，为矿区安全生产提供综合技术保障。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于叠加放水试验获取含水层连续参数的方法，其特征在于：包括如下步骤：
在研究区域实施叠加放水试验；
建立研究区域的数值模型并划分网格离散模型；
设定含水层参数的随机分布函数值；
选择适用的协方差函数描述含水层参数的变异程度；
利用放水层析法迭代计算含水层参数的连续分布。
2.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：实施叠加放水试验时，潜在富水区的钻孔加倍布设。
3.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：在潜在富水区的钻孔密度加密2-3倍。
4.根据权利要求2所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：在研究区域分阶段分区域依次实施放水，所有观测孔同步测量水压，每次试验的前期加密观测水压，中后期定时观测；
进一步的，叠加放水试验的次数大于或等于3次。
5.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：网格离散模型中的网格尺寸小于或等于相关长度的1/2；
进一步的，设定研究区域为非均质饱和含水层，根据研究区域的规模及对反演精度的需求，确定离散网格的规格。
6.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：设定含水层参数的随机分布函数值时，根据已知的地质信息预估水文地质参数的均值，将参数的初始估计值设定为均质分布；
或，在均值基础上，为每个离散元的参数附加符合正态分布的扰动值，参数的初始估计值设定为非均质分布。
7.根据权利要求6所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：地质信息包括地质填图、钻孔岩芯取样及抽水试验结果。
8.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：选择适用的协方差函数描述含水层参数的变异程度时，通过方差和相关长度描述含水层参数的变异程度。
9.根据权利要求8所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：研究区域在无经验协方差函数时，根据含水层类型，方差和相关长度的赋值如下：对于孔隙含水层，方差设为0.5-1，相关长度设为20-50m；对于裂隙含水层，方差设为1-2，相关长度设为50-200m；对于岩溶含水层，方差设为1-2，相关长度设为100-500m。
10.根据权利要求1所述的基于叠加放水试验获取含水层连续参数的方法，其特征在于：利用放水层析法迭代计算含水层参数的连续分布的具体步骤如下：
将水压数据整理为相等时间间隔的数据集，并进行平滑降噪，提取数据的明显变化规律；
确定钻孔在反演区域中的相对位置，平滑处理过后的水压值按照时间顺序依次赋值于对应钻孔结点处，边界包括定水头边界和定流量边界，定水头边界保证含水层始终维持饱和；
设定试验开始时间为反演计算的初始时间，初始时间区域内部结点的水头值与定水头边界的水头值数值相等；
根据放水试验的顺序逐次反演并迭代计算，获取含水层连续参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010161020.2A CN111401716B (zh) | 2020-03-10 | 2020-03-10 | 一种基于叠加放水试验获取含水层连续参数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010161020.2A CN111401716B (zh) | 2020-03-10 | 2020-03-10 | 一种基于叠加放水试验获取含水层连续参数的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111401716A true CN111401716A (zh) | 2020-07-10 |
| CN111401716B CN111401716B (zh) | 2022-03-15 |
Family
ID=71430759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010161020.2A Active CN111401716B (zh) | 2020-03-10 | 2020-03-10 | 一种基于叠加放水试验获取含水层连续参数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111401716B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111767670A (zh) * | 2020-07-17 | 2020-10-13 | 中煤科工集团西安研究院有限公司 | 基于井下仰斜单孔放水试验的含水层地质参数反演方法 |
| CN111881974A (zh) * | 2020-07-27 | 2020-11-03 | 河南理工大学 | 一种基于Piper-PCA-FCL判别模型的突水水源识别方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104216023A (zh) * | 2014-09-01 | 2014-12-17 | 山东科技大学 | 矿井掘进巷道高密度三维直流勘探方法 |
| CN109869187A (zh) * | 2019-03-11 | 2019-06-11 | 山东科技大学 | 基于水文探查及放水试验的含水层疏放可行性分析方法 |
-
2020
- 2020-03-10 CN CN202010161020.2A patent/CN111401716B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104216023A (zh) * | 2014-09-01 | 2014-12-17 | 山东科技大学 | 矿井掘进巷道高密度三维直流勘探方法 |
| CN109869187A (zh) * | 2019-03-11 | 2019-06-11 | 山东科技大学 | 基于水文探查及放水试验的含水层疏放可行性分析方法 |
Non-Patent Citations (5)
| Title |
|---|
| YONG-LI GUO 等: "The Risk Assessment of Groundwater Pollution in the Dawu Water Source", 《2012 INTERNATIONAL CONFERENCE ON BIOMEDICAL ENGINEERING AND BIOTECHNOLOGY》 * |
| 董艳辉 等: "应用水力层析法刻画含水层非均质性", 《工程勘察》 * |
| 蒋立群 等: "水力层析法与克立金法估算非均质含水层渗透系数场比较", 《地球科学》 * |
| 赵瑞珏 等: "基于水力层析法的某煤矿承压含水层叠加放水试验分析", 《水文地质工程地质》 * |
| 邵红旗 等: "一种放水试验分析方法及其应用", 《水文地质工程地质》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111767670A (zh) * | 2020-07-17 | 2020-10-13 | 中煤科工集团西安研究院有限公司 | 基于井下仰斜单孔放水试验的含水层地质参数反演方法 |
| CN111767670B (zh) * | 2020-07-17 | 2024-04-09 | 中煤科工集团西安研究院有限公司 | 基于井下仰斜单孔放水试验的含水层地质参数反演方法 |
| CN111881974A (zh) * | 2020-07-27 | 2020-11-03 | 河南理工大学 | 一种基于Piper-PCA-FCL判别模型的突水水源识别方法 |
| CN111881974B (zh) * | 2020-07-27 | 2023-06-09 | 河南理工大学 | 一种基于Piper-PCA-FCL判别模型的突水水源识别方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111401716B (zh) | 2022-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110795878B (zh) | 一种隧道涌水量预测方法 | |
| US9988900B2 (en) | Method of geometric evaluation of hydraulic fractures by using pressure changes | |
| US20160376885A1 (en) | Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs | |
| US6886632B2 (en) | Estimating formation properties in inter-well regions by monitoring saturation and salinity front arrivals | |
| CN105223612B (zh) | 一种基于地震信息的煤矿水害预测评价方法 | |
| GB2528384A (en) | Completion design based on logging while drilling (LWD) data | |
| CN109541685B (zh) | 一种河道砂体识别方法 | |
| CN111401716B (zh) | 一种基于叠加放水试验获取含水层连续参数的方法 | |
| Aslanyan et al. | Waterflood Study of High Viscosity Saturated Reservoir with Multiwell Retrospective Testing and Cross-Well Pressure Pulse-Code Testing | |
| CN103088803B (zh) | 鉴别岩溶塌陷致塌因素的方法 | |
| CN107895092B (zh) | 一种基于复杂非线性注采建模的井间连通定量评价方法 | |
| CN113742995B (zh) | 一种基于煤矿大数据的矿井涌水量预测方法及系统 | |
| CN105467438A (zh) | 一种基于三模量的泥页岩地应力三维地震表征方法 | |
| Comerford et al. | Controls on geothermal heat recovery from a hot sedimentary aquifer in Guardbridge, Scotland: Field measurements, modelling and long term sustainability | |
| Landa et al. | Reservoir characterization constrained to well test data: a field example | |
| CN111734394B (zh) | 一种确定致密油藏压裂井不定常流井底压力的方法 | |
| CN111415038A (zh) | 一种多工作面采空区涌水量随采精细化预测方法 | |
| Lee et al. | Aquifer analysis and modeling in a fractured, heterogeneous medium | |
| CN110320574B (zh) | 基于缓坡三角洲薄层砂体刻画的方法 | |
| Aronovici et al. | Soil‐permeability as a criterion for drainage‐design | |
| Winkler et al. | Scale dependent hydraulic investigations of faulted crystalline rocks—Examples from the Eastern Alps, Austria | |
| CN110687603A (zh) | 一种海上油田储层内部渗流屏障的地质建模方法 | |
| CN111236934B (zh) | 水淹级别确定方法和装置 | |
| Moran | Interpretation of long-term Grande Ronde aquifer testing in the Palouse Basin of Idaho and Washington | |
| CN113094864B (zh) | 强非均质碳酸盐岩水平井分段设计方法 |
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 |