CN111077587A - 一种地下岩溶结构精细刻画方法及装置 - Google Patents
一种地下岩溶结构精细刻画方法及装置 Download PDFInfo
- Publication number
- CN111077587A CN111077587A CN201911391156.6A CN201911391156A CN111077587A CN 111077587 A CN111077587 A CN 111077587A CN 201911391156 A CN201911391156 A CN 201911391156A CN 111077587 A CN111077587 A CN 111077587A
- Authority
- CN
- China
- Prior art keywords
- water level
- value
- well
- pumping
- level change
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 137
- 238000005086 pumping Methods 0.000 claims abstract description 99
- 230000008859 change Effects 0.000 claims abstract description 80
- 230000006870 function Effects 0.000 claims description 70
- 239000002349 well water Substances 0.000 claims description 28
- 235000020681 well water Nutrition 0.000 claims description 28
- 238000012512 characterization method Methods 0.000 claims description 24
- 230000000737 periodic effect Effects 0.000 claims description 10
- 238000012886 linear function Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 238000012360 testing method Methods 0.000 description 19
- 239000000243 solution Substances 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种地下岩溶结构精细刻画方法及装置。包括：将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。本发明通过将抽水井抽水流速设置为以时间为自变量的余弦函数，对观测井内水位变化值进行分析，结合抽水井内水位变化值(相位和振幅)，能够准确快速对地下岩溶结构进行精细刻画。
Description
技术领域
本发明涉及地下岩溶结构刻画技术领域，尤其涉及一种地下岩溶结构精细刻画方法及装置。
背景技术
许多年来，国内外学者致力于对地层岩溶结构刻画技术的研究，此技术在水资源开发与保护，石油开采，地热能开发等领域都有广泛应用。井间抽水测试方法就是地层参数刻画技术的其中一种。所谓的井间抽水测试技术，是指针对目标区域，从抽水井抽取一定水量而在某距离外的观测井测定不同时间井内水位的变化，并利用各种地下水流理论式或图解法分析抽水试验的结果从而获得地层参数的方法。
但是，传统的井间抽水试验，在抽水井内的抽水速度或者井内水位是恒定的。这种方法存在着以下的缺点。第一，传统抽水测试会受到一系列的外在因素影响，导致结果不准，例如：地下水流速、河流、蒸发潮汐等因素。第二，传统抽水试验往往采用是log坐标表示时间，测试周期长，测试时间往往需要几天甚至一周以上。所以，亟需一种周期短、抗干扰性强的抽水试验技术来实现地下岩溶结构的精准刻画。
上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。
发明内容
有鉴于此，本发明提出了一种地下岩溶结构精细刻画方法及装置，旨在解决现有技术无法提供一种周期短、抗干扰性强的抽水试验技术来实现地下岩溶结构的精准刻画的技术问题。
本发明的技术方案是这样实现的：
一方面，本发明提供了一种地下岩溶结构精细刻画方法，所述地下岩溶结构精细刻画方法包括以下步骤：
S1，将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；
S2，设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
在以上技术方案的基础上，优选的，步骤S1中，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，还包括以下步骤，所述待计算数据包括：流速周期变化振幅、流速平均值以及注水周期，所述流速平均值根据历史流速通过平均算法获得。
在以上技术方案的基础上，优选的，还包括以下步骤，所述时间流速算法为：
Q(t)＝-QAcos(ωt)+Qm；
其中，QA代表流速周期变化振幅，ω代表频率，
Qm代表流速平均值，t代表时间，Q(t)代表流速，T代表注水周期。
在以上技术方案的基础上，优选的，步骤S1中，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据，还包括以下步骤，根据该水位变化数据获取水位变化曲线，将观测井水位变化曲线转化为曲线图形，判断该曲线图形是否存在周期性变化，若存在，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值；若不存在，则重新选择抽水井。
在以上技术方案的基础上，优选的，步骤S2中，设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画之前，还包括以下步骤，获取曲线函数格式，根据该曲线函数格式将水位变化数据转化为水位变化曲线函数，通过傅里叶变换对所述水位变化曲线函数进行分解，所述水位变化曲线函数分解为：
h(x,y,t)＝hosc(x,y,t)+hlin(x,y,t)；
其中，h(x,y,t)代表水位变化函数，hosc(x,y,t)代表水位变化三角函数，hlin(x,y,t)代表水位变化线性函数，x代表水位的宽度，y代表变化水位的高度，t代表时间。
在以上技术方案的基础上，优选的，步骤S2中，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，还包括以下步骤，获取观测井水位变化三角函数以及抽水井水位变化三角函数，将观测井水位变化三角函数与抽水井水位变化三角函数进行比较，获取观测井水位变化三角函数的振幅占抽水井水位变化三角函数的振幅的百分比数值以及观测井水位变化三角函数相对于抽水井水位变化三角函数的相位差值。
在以上技术方案的基础上，优选的，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画，还包括以下步骤，设定比较数值阈值以及对应的判断结果，所述判断结果分为三类，对振幅的百分比数值进行判断，当振幅的百分比数值为零时，获取第一类判断结果，并对地下岩溶结构进行刻画；当当振幅的百分比数值不为零时，将相位差值与比较数值阈值进行比较，当相位差值小于比较数值阈值时，获取第二类判断结果，并对地下岩溶结构进行刻画；当相位差值大于比较数值阈值时，获取第三类判断结果，并对地下岩溶结构进行刻画。
更进一步优选的，所述地下岩溶结构刻画装置包括：
记录模块，用于将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；
刻画模块，用于设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
第二方面，所述地下岩溶结构精细刻画方法还包括一种设备，所述设备包括：存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的地下岩溶结构精细刻画方法程序，所述地下岩溶结构精细刻画方法程序配置为实现如上文所述的地下岩溶结构精细刻画方法的步骤。
第三方面，所述地下岩溶结构精细刻画方法还包括一种介质，所述介质为计算机介质，所述计算机介质上存储有地下岩溶结构精细刻画方法程序，所述地下岩溶结构精细刻画方法程序被处理器执行时实现如上文所述的地下岩溶结构精细刻画方法的步骤。
本发明的一种地下岩溶结构精细刻画方法相对于现有技术具有以下有益效果：
(1)通过将抽水井抽水流速设置为以时间为自变量的余弦函数，不仅可以提供更多的有效信息，而且抗干扰因素比较强，减少了测试的时间，提高了整个地下岩溶结构刻画过程的效率；
(2)通过利用傅里叶变换将曲线函数分解为一个三角函数和一个线性函数的叠加，根据三角函数来对岩溶结构进行判断，能够精准且快速获取岩溶结构，提高了整个刻画过程的效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例方案涉及的硬件运行环境的设备的结构示意图；
图2为本发明地下岩溶结构精细刻画方法第一实施例的流程示意图；
图3为本发明地下岩溶结构精细刻画方法第一实施例的功能模块示意图。
具体实施方式
下面将结合本发明实施方式，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本发明一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
如图1所示，该设备可以包括：处理器1001，例如中央处理器(Central ProcessingUnit，CPU)，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(WIreless-FIdelity，WI-FI)接口)。存储器1005可以是高速的随机存取存储器(Random Access Memory，RAM)存储器，也可以是稳定的非易失性存储器(Non-Volatile Memory，NVM)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的结构并不构成对设备的限定，在实际应用中设备可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种介质的存储器1005中可以包括操作系统、网络通信模块、用户接口模块以及地下岩溶结构精细刻画方法程序。
在图1所示的设备中，网络接口1004主要用于建立设备与存储地下岩溶结构精细刻画方法系统中所需的所有数据的服务器的通信连接；用户接口1003主要用于与用户进行数据交互；本发明地下岩溶结构精细刻画方法设备中的处理器1001、存储器1005可以设置在地下岩溶结构精细刻画方法设备中，所述地下岩溶结构精细刻画方法设备通过处理器1001调用存储器1005中存储的地下岩溶结构精细刻画方法程序，并执行本发明实施提供的地下岩溶结构精细刻画方法。
结合图2，图2为本发明地下岩溶结构精细刻画方法第一实施例的流程示意图。
本实施例中，所述地下岩溶结构精细刻画方法包括以下步骤：
S10：将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据。
应当理解的是，将水井分为观测井以及抽水井两组之间，会人为选定要测试底层参数的目标区域，根据要研究的目标来开钻探测井网，同时会设计2-4组抽水试验，来确保刻画的精确性，每一组试验会选定一口井作为抽水井，其它的井作为观测井，为了方便计算，系统会给所有的井分配编号，编号从1号开始，根据井的数量来进行编号，然后会确定待计算数据，所述待计算数据包括：流速周期变化振幅、流速平均值以及注水周期，所述流速平均值根据历史流速通过平均算法获得，历史流速是指系统内部存储的抽水井的水流速度。
应当理解的是，之后会进行周期性的抽水测试，建立时间流速算法，根据这个时间流速算法来进行抽水，然后在观测井内观测并记录抽水井以及观测井内的水位变化数值。
应当理解的是，所述时间流速算法为：
Q(t)＝-QAcos(ωt)+Qm；
其中，QA代表流速周期变化振幅，ω代表频率，
Qm代表流速平均值，t代表时间，Q(t)代表流速，T代表注水周期。
应当理解的是，根据实际情况的不同，每组抽水试验至少采用两个不同的周期T进行实验。一个短周期抽水试验足以找出那些恰好位于岩溶结构上的观测井，但是增加一个长周期的抽水试验可以准确地找出那些不在岩溶结构上但是就在岩溶结构附近的观测井。根据现场操作经验，选用1min和5min作为周期效果较好。
应当理解的是，系统还会获取曲线函数格式，根据该曲线函数格式将水位变化数据转化为水位变化曲线函数，然后观测井水位变化曲线都没有出现明显的周期性波动，则此组抽水试验无效，则需要重新选定抽水井以及观测井进行试验；如果出现明显周期性波动，说明抽水井的位置成功打在了地下岩溶结构上。此组抽水试验可以为我们提供有效信息。
S20：设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
应当理解的是，为了后续方便进行比较，系统会将所有井内的水位变化数据转化为水位变化曲线函数，然后通过傅里叶变换对所述水位变化曲线函数进行分解，将水位变化函数分解为一个三角函数和一个线性函数的叠加，所述水位变化曲线函数分解为：
h(x,y,t)＝hosc(x,y,t)+hlin(x,y,t)；
其中，h(x,y,t)代表水位变化函数，hosc(x,y,t)代表水位变化三角函数，hlin(x,y,t)代表水位变化线性函数，x代表水位的宽度，y代表变化水位的高度，t代表时间。
应当理解的是，傅里叶变换是能将满足一定条件的某个函数表示成三角函数(正弦和/或余弦函数)或者它们的积分的线性组合。在不同的研究领域，傅立叶变换具有多种不同的变体形式，如连续傅立叶变换和离散傅立叶变换。
应当理解的是，系统会获取观测井水位变化三角函数以及抽水井水位变化三角函数，将观测井水位变化三角函数与抽水井水位变化三角函数进行比较，获取观测井水位变化三角函数的振幅占抽水井水位变化三角函数的振幅的百分比数值以及观测井水位变化三角函数相对于抽水井水位变化三角函数的相位差值，如表1所示。
表1幅度相对大小以及相位差值表
应当理解的是，系统会设定相对应的判断条件以及判断结果，判断条件如表2所示，设定比较阈值时判断条件之一，这个阈值由用户自己设定，对振幅的百分比数值进行判断，当振幅的百分比数值为零时，获取第一类判断结果，并对地下岩溶结构进行刻画；当当振幅的百分比数值不为零时，将相位差值与比较数值阈值进行比较，当相位差值小于比较数值阈值时，获取第二类判断结果，并对地下岩溶结构进行刻画；当相位差值大于比较数值阈值时，获取第三类判断结果，并对地下岩溶结构进行刻画。
表2判断条件
应当理解的是，对于表2的判断条件，可以获得如下结论，提高抽水周期导致幅度增加，相位变化下降，观测井距离岩溶结构越远，则幅度越小，相位变化值越高。
应当理解的是，通过表2，以表1为例，对表1中的观测井进行判断，如果1号井为抽水井，则2号观测井以及3号观测井都属于第二类；如果2号井为抽水井，则1号观测井以及3号观测井都属于第三类，由此，在结合当前的实际地理情况，就可以刻画出地下岩溶结构。
需要说明的是，以上仅为举例说明，并不对本申请的技术方案构成任何限定。
通过上述描述不难发现，本实施例通过将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。本实施例通过将抽水井抽水流速设置为以时间为自变量的余弦函数，对观测井内水位变化值进行分析，结合抽水井内水位变化值(相位和振幅)，能够准确快速对地下岩溶结构进行精细刻画。
此外，本发明实施例还提出一种地下岩溶结构刻画装置。如图3所示，该地下岩溶结构刻画装置包括：记录模块10、刻画模块20。
记录模块10，用于设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画；
刻画模块20，用于设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
此外，需要说明的是，以上所描述的装置实施例仅仅是示意性的，并不对本发明的保护范围构成限定，在实际应用中，本领域的技术人员可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的，此处不做限制。
另外，未在本实施例中详尽描述的技术细节，可参见本发明任意实施例所提供的地下岩溶结构精细刻画方法，此处不再赘述。
此外，本发明实施例还提出一种介质，所述介质为计算机介质，所述计算机介质上存储有地下岩溶结构精细刻画方法程序，所述岩地下岩溶结构精细刻画方法程序被处理器执行时实现如下操作：
S1，将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；
S2，设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
所述待计算数据包括：流速周期变化振幅、流速平均值以及注水周期，所述流速平均值根据历史流速通过平均算法获得。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
所述时间流速算法为：
Q(t)＝-QAcos(ωt)+Qm；
其中，QA代表流速周期变化振幅，ω代表频率，
Qm代表流速平均值，t代表时间，Q(t)代表流速，T代表注水周期。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
根据该水位变化数据获取水位变化曲线，将观测井水位变化曲线转化为曲线图形，判断该曲线图形是否存在周期性变化，若存在，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值；若不存在，则重新选择抽水井。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
获取曲线函数格式，根据该曲线函数格式将水位变化数据转化为水位变化曲线函数，通过傅里叶变换对所述水位变化曲线函数进行分解，所述水位变化曲线函数分解为：
h(x,y,t)＝hosc(x,y,t)+hlin(x,y,t)；
其中，h(x,y,t)代表水位变化函数，hosc(x,y,t)代表水位变化三角函数，hlin(x,y,t)代表水位变化线性函数，x代表水位的宽度，y代表变化水位的高度，t代表时间。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
获取观测井水位变化三角函数以及抽水井水位变化三角函数，将观测井水位变化三角函数与抽水井水位变化三角函数进行比较，获取观测井水位变化三角函数的振幅占抽水井水位变化三角函数的振幅的百分比数值以及观测井水位变化三角函数相对于抽水井水位变化三角函数的相位差值。
进一步地，所述地下岩溶结构精细刻画方法程序被处理器执行时还实现如下操作：
设定比较数值阈值以及对应的判断结果，所述判断结果分为三类，对振幅的百分比数值进行判断，当振幅的百分比数值为零时，获取第一类判断结果，并对地下岩溶结构进行刻画；当当振幅的百分比数值不为零时，将相位差值与比较数值阈值进行比较，当相位差值小于比较数值阈值时，获取第二类判断结果，并对地下岩溶结构进行刻画；当相位差值大于比较数值阈值时，获取第三类判断结果，并对地下岩溶结构进行刻画。
以上所述仅为本发明的较佳实施方式而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种地下岩溶结构精细刻画方法，其特征在于：包括以下步骤；
S1，将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；
S2，设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
2.如权利要求1所述的地下岩溶结构精细刻画方法，其特征在于：步骤S1中，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，还包括以下步骤，所述待计算数据包括：流速周期变化振幅、流速平均值以及注水周期，所述流速平均值根据历史流速通过平均算法获得。
3.如权利要求2所述的地下岩溶结构精细刻画方法，其特征在于：还包括以下步骤，所述时间流速算法为：
Q(t)＝-QAcos(ωt)+Qm；
其中，QA代表流速周期变化振幅，ω代表频率，
Qm代表流速平均值，t代表时间，Q(t)代表流速，T代表注水周期。
4.如权利要求3所述的地下岩溶结构精细刻画方法，其特征在于：步骤S1中，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据，还包括以下步骤，根据该水位变化数据获取水位变化曲线，将观测井水位变化曲线转化为曲线图形，判断该曲线图形是否存在周期性变化，若存在，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值；若不存在，则重新选择抽水井。
5.如权利要求4所述的地下岩溶结构精细刻画方法，其特征在于：步骤S2中，设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画之前，还包括以下步骤，获取曲线函数格式，根据该曲线函数格式将水位变化数据转化为水位变化曲线函数，通过傅里叶变换对所述水位变化曲线函数进行分解，所述水位变化曲线函数分解为：
h(x,y,t)＝hosc(x,y,t)+hlin(x,y,t)；
其中，h(x,y,t)代表水位变化函数，hosc(x,y,t)代表水位变化三角函数，hlin(x,y,t)代表水位变化线性函数，x代表水位的宽度，y代表变化水位的高度，t代表时间。
6.如权利要求5所述的地下岩溶结构精细刻画方法，其特征在于：步骤S2中，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，还包括以下步骤，获取观测井水位变化三角函数以及抽水井水位变化三角函数，将观测井水位变化三角函数与抽水井水位变化三角函数进行比较，获取观测井水位变化三角函数的振幅占抽水井水位变化三角函数的振幅的百分比数值以及观测井水位变化三角函数相对于抽水井水位变化三角函数的相位差值。
7.如权利要求6所述的地下岩溶结构精细刻画方法，其特征在于：将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画，还包括以下步骤，设定比较数值阈值以及对应的判断结果，所述判断结果分为三类，对振幅的百分比数值进行判断，当振幅的百分比数值为零时，获取第一类判断结果，并对地下岩溶结构进行刻画；当当振幅的百分比数值不为零时，将相位差值与比较数值阈值进行比较，当相位差值小于比较数值阈值时，获取第二类判断结果，并对地下岩溶结构进行刻画；当相位差值大于比较数值阈值时，获取第三类判断结果，并对地下岩溶结构进行刻画。
8.一种地下岩溶结构刻画装置，其特征在于，所述地下岩溶结构刻画装置包括：
记录模块，用于将水井分为观测井以及抽水井两组，建立时间流速算法，获取待计算数据，根据时间流速算法进行计算，获取流速数值，根据该流速数值对抽水井进行抽水，记录所有水井内水位变化数据；
刻画模块，用于设定比较数值阈值，根据观测井水位变化数值以及抽水井水位变化数值确定获取待比较数值，将比较数值阈值与待比较数值进行比较，根据比较结果对地下岩溶结构进行刻画。
9.一种设备，其特征在于，所述设备包括：存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的岩溶地下岩溶结构精细刻画方法程序，所述岩溶地下岩溶结构精细刻画方法程序配置为实现如权利要求1至7任一项所述的岩溶地下岩溶结构精细刻画方法的步骤。
10.一种介质，其特征在于，所述介质为计算机介质，所述计算机介质上存储有岩溶地下岩溶结构精细刻画方法程序，所述岩溶地下岩溶结构精细刻画方法程序被处理器执行时实现如权利要求1至7任一项所述的岩溶地下岩溶结构精细刻画方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911391156.6A CN111077587B (zh) | 2019-12-30 | 2019-12-30 | 一种地下岩溶结构精细刻画方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911391156.6A CN111077587B (zh) | 2019-12-30 | 2019-12-30 | 一种地下岩溶结构精细刻画方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111077587A true CN111077587A (zh) | 2020-04-28 |
| CN111077587B CN111077587B (zh) | 2022-04-15 |
Family
ID=70319842
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911391156.6A Active CN111077587B (zh) | 2019-12-30 | 2019-12-30 | 一种地下岩溶结构精细刻画方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111077587B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08311852A (ja) * | 1995-05-17 | 1996-11-26 | Shimizu Corp | 揚水試験方法 |
| JP2006161340A (ja) * | 2004-12-03 | 2006-06-22 | Tone Geo Tech Co Ltd | 揚水試験解析方法及び揚水試験解析用記憶媒体 |
| JP2008169622A (ja) * | 2007-01-11 | 2008-07-24 | Toho Chisui Kk | 揚水試験方法と揚水試験装置 |
| CN102213775A (zh) * | 2011-04-12 | 2011-10-12 | 河海大学 | 获取分布式水文模型下垫面条件的地下水连通性探测方法 |
| CN102435543A (zh) * | 2011-12-05 | 2012-05-02 | 湖南科技大学 | 在线全孔连续检测的稳定流抽水试验设备及其检测方法 |
| KR101205445B1 (ko) * | 2012-05-29 | 2012-11-27 | 주식회사 지오그린이십일 | 지하수 양수시험 장치 및 방법 |
| CN103303986A (zh) * | 2013-07-05 | 2013-09-18 | 华北电力大学 | 基于周期性原位抽灌一体化修复地下水污染的装置及方法 |
| CN103410175A (zh) * | 2013-08-05 | 2013-11-27 | 中铁十二局集团有限公司 | 一种岩溶富水地质地下结构防水施工渗水控制施工方法 |
| CN105297753A (zh) * | 2015-10-16 | 2016-02-03 | 上海交通大学 | 一种在止水帷幕作用下水文地质参数的确定方法 |
| TW201723532A (zh) * | 2015-12-28 | 2017-07-01 | 逸奇科技股份有限公司 | 一種地下水位之評估方法 |
| CN107289997A (zh) * | 2017-05-05 | 2017-10-24 | 济南轨道交通集团有限公司 | 一种岩溶裂隙水探测系统及方法 |
| JP2019127701A (ja) * | 2018-01-22 | 2019-08-01 | 株式会社竹中工務店 | 揚水制御装置及び揚水制御モデル学習装置 |
-
2019
- 2019-12-30 CN CN201911391156.6A patent/CN111077587B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08311852A (ja) * | 1995-05-17 | 1996-11-26 | Shimizu Corp | 揚水試験方法 |
| JP2006161340A (ja) * | 2004-12-03 | 2006-06-22 | Tone Geo Tech Co Ltd | 揚水試験解析方法及び揚水試験解析用記憶媒体 |
| JP2008169622A (ja) * | 2007-01-11 | 2008-07-24 | Toho Chisui Kk | 揚水試験方法と揚水試験装置 |
| CN102213775A (zh) * | 2011-04-12 | 2011-10-12 | 河海大学 | 获取分布式水文模型下垫面条件的地下水连通性探测方法 |
| CN102435543A (zh) * | 2011-12-05 | 2012-05-02 | 湖南科技大学 | 在线全孔连续检测的稳定流抽水试验设备及其检测方法 |
| KR101205445B1 (ko) * | 2012-05-29 | 2012-11-27 | 주식회사 지오그린이십일 | 지하수 양수시험 장치 및 방법 |
| CN103303986A (zh) * | 2013-07-05 | 2013-09-18 | 华北电力大学 | 基于周期性原位抽灌一体化修复地下水污染的装置及方法 |
| CN103410175A (zh) * | 2013-08-05 | 2013-11-27 | 中铁十二局集团有限公司 | 一种岩溶富水地质地下结构防水施工渗水控制施工方法 |
| CN105297753A (zh) * | 2015-10-16 | 2016-02-03 | 上海交通大学 | 一种在止水帷幕作用下水文地质参数的确定方法 |
| TW201723532A (zh) * | 2015-12-28 | 2017-07-01 | 逸奇科技股份有限公司 | 一種地下水位之評估方法 |
| CN107289997A (zh) * | 2017-05-05 | 2017-10-24 | 济南轨道交通集团有限公司 | 一种岩溶裂隙水探测系统及方法 |
| JP2019127701A (ja) * | 2018-01-22 | 2019-08-01 | 株式会社竹中工務店 | 揚水制御装置及び揚水制御モデル学習装置 |
Non-Patent Citations (3)
| Title |
|---|
| 吕全标等: "基于钻孔抽水试验和示踪试验的岩溶地区含水层结构研究", 《中国岩溶》 * |
| 朱昊: "岩溶含水介质刻画 ——以桂林毛村流域为例", 《中国优秀硕士学位论文全文数据库（基础科学辑）》 * |
| 郭续磊等: "新型地下水流速流向测量技术及其在岩溶区调查中的应用", 《地质科技情报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111077587B (zh) | 2022-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hagen et al. | Two‐dimensional, unstructured mesh generation for tidal models | |
| Tee | Tide-induced residual current, a 2-D nonlinear numerical tidal model | |
| US8862450B2 (en) | Method of developing a petroleum reservoir from history matching of production data and seismic data | |
| Holman et al. | A parametric model for barred equilibrium beach profiles | |
| CN109543356A (zh) | 考虑空间非平稳性的海洋内部温盐结构遥感反演方法 | |
| Dickinson et al. | A screening tool for delineating subregions of steady recharge within groundwater models | |
| CN115203945B (zh) | 一种水文地质参数估计方法、装置、设备和存储介质 | |
| Limber et al. | Modeling large-scale shoreline change caused by complex bathymetry in low-angle wave climates | |
| CN106991279A (zh) | 一种致密砂岩气成藏概率的确定方法和装置 | |
| CN111475920A (zh) | 一种深水盆地古水深的获取方法、系统、电子设备及存储介质 | |
| Elken et al. | Circulation patterns in the Gulf of Finland derived from the EOF analysis of model results | |
| Xiong et al. | On estimation of coastal wave parameters and wave‐induced shear stresses | |
| CN103616496B (zh) | 一种测定滩涂湿地早晚潮间期土壤co2通量的方法 | |
| CN111077587B (zh) | 一种地下岩溶结构精细刻画方法及装置 | |
| Cael et al. | A theory for the relationship between lake surface area and maximum depth | |
| Ibrahim et al. | Comparison between inverse distance weighted (IDW) and Kriging | |
| CN111985747B (zh) | 油藏开发方法及装置 | |
| Kambesis et al. | Morphometric analysis of cave patterns using fractal indices | |
| CN104063614B (zh) | 一种基于odp的海底沙波特征自动识别方法 | |
| Ding et al. | Cross-shore sediment transport for modeling long-term shoreline evolution | |
| Alghalandis et al. | Similarity analysis of discrete fracture networks | |
| CN108510134B (zh) | 一种确定地质变量分位值的方法和系统 | |
| CN112001055B (zh) | 一种基于微构造的低幅稀油油藏含水率预测方法 | |
| Crenshaw et al. | Water table profiles in the vicinity of horizontal drains | |
| Lecointre et al. | Depth dependence of westward-propagating North Atlantic features diagnosed from altimetry and a numerical 1/6 model |
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 |