CN103149600B - 一种基于优化控制点确定水文地质参数的自动配线方法 - Google Patents
一种基于优化控制点确定水文地质参数的自动配线方法 Download PDFInfo
- Publication number
- CN103149600B CN103149600B CN201310048344.5A CN201310048344A CN103149600B CN 103149600 B CN103149600 B CN 103149600B CN 201310048344 A CN201310048344 A CN 201310048344A CN 103149600 B CN103149600 B CN 103149600B
- Authority
- CN
- China
- Prior art keywords
- point
- measured data
- reference mark
- wiring method
- automatic wiring
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于优化控制点确定水文地质参数的自动配线方法,包括如下步骤:进行承压含水层非稳定流抽水试验,得试验资料;基于非稳定流承压含水层抽水试验资料,绘制实测数据曲线与配线法标准曲线,给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合;选取匹配点计算水文地质参数。本发明提出的计算机自动配线的新方法,人机交互简单,接近传统配线法的基本操作模式;通过控制点带动所有实测点,可以大大减小沿标准曲线移动的次数,相当于寻找了一条最优路径移动实测数据点,在保证精度的同时,缩短了后台操作时间,提高了工作效率。
Description
技术领域
本发明涉及一种确定承压含水层水文地质参数的方法,尤其是涉及一种基于优化控制点实现自动配线确定承压含水层导水系数、贮水系数等水文地质参数的方法。
背景技术
近年来,随着我国人口的迅速增长和经济的快速发展,城市建设规模和速度前所未有,城市建设用地日益紧张,地下空间的开发和利用也日益紧迫,北京和沿海一些大城市地下轨道交通和车站、大型地下商场、地下停车场等大型深基坑工程层出不穷,这些深基坑的施工成败与能否较好的控制地下水关系密切,为降低深基坑工程的设计与施工风险,有必要通过抽水试验确定各含水层,特别是承压含水层的水文地质参数。此外,抽水试验也是煤炭资源等地质勘探的重要手段,其目的是研究含水层重要水文地质特征,取得含水层水文地质参数,评价含水层的富水性,并为预计矿井涌水量与对地下水综合利用的评价提供资料。因此,承压含水层水文地质参数的确定对揭示抽取地下水引起地下水变化和运动的基本规律有着重要的意义,而且为开采地下水和石油以及地下工程建设等生产活动提供理论依据和技术保障。
非稳定流承压完整井求取水文地质参数的方法主要有Jacob直线法、配线法和水位恢复法,其中配线法是最为常用的方法。传统方法采用人工配线,步骤相对繁琐,工作量较大,而且利用配线观察者的目测进行配线,存在较大的随意性,从而引入不必要的误差。随着计算机技术的发展,计算机配线早年主要停留在手动配线的阶段,通过程序实现了数据的自动录入,但曲线的拟合仍然依据主观的判断;近年来兴起了智能算法,如基于配线法的BP神经网络求解Theis模型(2010)等,但是这些智能算法已失去了配线法实现过程的直观性。
发明内容
发明目的:本发明基于配线法标准曲线与实测曲线基本重合,采用控制点带动所有实测点沿标准曲线移动,寻找离差平方和最小值的原理,提出一种实现自动配线确定承压含水层水文地质参数的方法,并结合抽水试验进行了验证。这一方法不仅精度高、易操作,而且具有计算快速,配线过程直观可见等优点。
技术方案:为了解决上述技术问题,采用如下技术方案:
一种基于优化控制点确定水文地质参数的自动配线方法,包括如下步骤:
A、进行承压含水层非稳定流抽水试验,得试验资料;
B、基于非稳定流承压含水层抽水试验资料,绘制实测数据曲线与配线法标准曲线,给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合;
C、选取匹配点计算水文地质参数。
为了提高配线的精度和易操作性,所述步骤B为:采用控制点和辅助控制点带动所有实测数据点沿标准曲线移动,寻找离差平方和最小值,实现自动配线;所述步骤C为:根据配线结果计算水文地质参数。
为了进一步提高配线的精度、易操作性及直观性,所述步骤B包括如下步骤:
B1、绘制配线法标准曲线,所述标准曲线由m个理论点连线而成;
B2、录入实测数据,所述实测数据为n个,由程序直接读取Excel表格录入;
B3、确定控制点坐标A(xA,yA),所述控制点A取n个实测数据的均值作为横坐标xA,然后在xA左右各取0.25,计算该范围内实测点纵坐标的均值yA,得到控制点A(xA,yA);由控制点A带动实测数据平移,每个理论点与控制点A在横、纵坐标方向有绝对距离和所述绝对距离即控制点至理论点的实际距离,与坐标无关,将控制点和所有实测数据点按平移量移动至标准曲线附近,计算每次的离差平方和fA(k)(k=1、2…m),确定fA的最小值,所述k表示第k次平移,同时也表示第k个理论点,k的取值范围为[1,m];
B4、取辅助控制点B(xB,yB),xB=2/3×(实测数据最大值-实测数据最小值)+实测数据最小值,然后在xB左右各取0.25,计算该范围内实测点纵坐标的均值yB,得到控制点B(xB,yB),用与步骤B3同样的方法寻找离差平方和fB的最小值,取fA、fB中较小值,按对应的平移量移动实测数据,即可完成配线。
上述步骤B3中,标准曲线(理论点组成)和实测曲线(实测数据点组成)的双对数坐标值可能不同。
所述步骤B3中,第k次平移的离差平方和fA(k),即每次移动后,计算第i个实测数据点与标准曲线上m个理论点的距离(k=1、2…m),所述为第i个数据点和第k个理论点的绝对距离,与坐标无关,找出(k=1、2…m)中的最小值,该值是第i个实测数据点至理论点的最短斜距,与求第i个实测数据点到标准曲线的最小距离等效,因为求离差平方和的目的在于找到实测曲线距离标准曲线最近的位置,所以该等效基本不影响自动配线的精度,而且不需要求点到曲线的距离,便于实现。
所有实测数据点用同样的方法确定到标准曲线的绝对距离,平方求和就可以得到离差平方和fA(k);步骤B4中,同样的方法,得到离差平方和fB(k)。
为了提高精度,上述离差平方和最小值的确定包括如下步骤:
1)k=1,即控制点A带动实测数据第一次平移,所述k同时是第k个理论点,k的取值范围为[1,m];
2)控制点A及n个实测数据点均平移所述为控制点A至第k个理论点的绝对距离,即控制点至理论点的实际距离,与坐标无关,作第k次平移;
3)(k=1、2…m;i=1、2…n),所述是第i个数据点和m个理论点距离的最小值(k=1、2…m),即第i个实测数据点至理论点的最短斜距,与求第i个实测数据点到标准曲线的最小距离等效,其中为第i个数据点和第k个理论点的绝对距离,与坐标无关;
4)(i=1、2…n),所述fA(k)为第k次平移后n个实测数据点与标准曲线的离差平方和;
5)当k<m时,k=k+1,即平移次数小于m次时,均由控制点A带动实测数据向下一个理论点平移,求得相应的离差平方和,直到m次平移完成后,终止循环;每一次平移后求得的离差平方和fA(k)(k=1、2…m)均存储在数列中;
6)循环终止,fA=min[fA(k)](k=1、2…m),得到控制点A带动实测数据平移m次的最小离差平方和。
上述步骤B4中按对应的平移量移动实测数据,即可完成配线,所述相应的平移量,如果fA≤fB,则相应的平移量为其中k为fA所对应的第k次平移;如果fA>fB,则相应的平移量为其中k为fB所对应的第k次平移,根据平移量得到最终配线位置。
上述步骤C中的水文地质参数一般包括导水系数和贮水系数。
导水系数和贮水系数的按如下方法计算:
在双对数坐标纸下,W(u),-1/u和s,-t/r2的曲线形态基本一致,当配线完成后,在曲线上或曲线任取一匹配点,,记下匹配点的对应坐标值:W(u),1/u,s和t/r2或t,代入(1)、(2),分别计算有关参数:
式中:W(u)—井函数;Q—抽水井流量,单位为m3/h;T—导水系数,单位为m2/d;μ*—含水层的贮水系数;s—抽水影响范围内,任意一点任一时刻的水位降深,单位为m;t—自抽水开始到计算时刻的时间,单位为min;r—计算点到抽水井距离,单位为m。
为了方便检测,上述步骤A抽水试验为在抽水孔周围配置若干个观测孔,抽水过程中,保持抽水量固定而观测地下水位变化的抽水试验;所述步骤A抽水试验资料包括抽水井流量、观测孔与抽水井的距离及观测孔降深随时间的变化。
有益效果:与现有技术相比,本发明所提供的一种基于优化控制点确定水文地质参数的自动配线方法具有以下优点:
①提出的计算机自动配线的新方法,人机交互简单,接近传统配线法的基本操作模式;
②相比将区域分成点的阵列,越加密,其精度越高,但后台操作次数剧增,比如等间距加密,每次加密后台操作次数增加(2n-1)2-n2,式中n为第n次加密时单行(或单列)的点数,在W(u)数值表中1/u有288个点,如果要达到该精度,则至少平移328329次,每次平移要计算每个实测点到标准曲线的距离,求离差平方和,这样运算量是很大的,甚至会导致死机。通过控制点带动所有实测点,可以大大减小沿标准曲线移动的次数,相当于寻找了一条最优路径移动实测数据点,在保证精度的同时,缩短了后台操作时间,提高了工作效率。
附图说明
图1为自动配线法的流程图;
图2(a)、2(b)、2(c)为控制点定位示意图,以控制点A为例;
图3为确定离差平方和最小值的流程图,以控制点A带动实测数据平移为例;
图4为配线法标准曲线;
图5为本发明实施例自动配线结果。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
确定承压含水层水文地质参数,首先进行承压含水层非稳定流抽水试验,给出抽水井流量、观测孔与抽水井距离及观测孔降深随时间变化的试验资料,采用Excel表格以规定格式记录试验数据;然后给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合,最终选取匹配点计算导水系数、贮水系数等水文地质参数的方法。
所述承压含水层非稳定流抽水试验,即在抽水孔周围配置若干个观测孔,抽水过程中,保持抽水量固定而观测地下水位变化的抽水试验,所述抽水孔的进水部分揭穿整个承压含水层。
所述通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合,最终选取匹配点计算导水系数、贮水系数等水文地质参数。如图1所示,其步骤为:
(1)绘制配线法标准曲线,所述标准曲线由m个理论值连线而成;
(2)录入实测数据,所述实测数据一共n个,由程序直接读取Excel表格录入;
(3)确定控制点坐标A(xA,yA)和辅助控制点B(xB,yB),所述控制点A取n个实测数据的均值作为横坐标xA,如图2(a)所示,取坐标均值可以避免选取某一异常数据点作为控制点;然后在xA左右各取0.25,如图2(b)所示,因为对数坐标轴的单位为1,相当于采用单位刻度的1/2,该范围不能取太小,否则可能数据点太少,容易受异常值的影响,该范围也不能取太大,否则控制点在曲线拐弯处,y坐标将发生偏离,根据标准曲线判断单位刻度的1/2范围内线段较平直,所以在x左右各取0.25;计算该范围内实测点(如图2(b)中的①②③)纵坐标的均值y,如图2(b)中所示;得到控制点A(xA,yA),如图2(c)中所示。所述辅助控制点B(xB,yB)的选取有利于提高计算精度,因为配线法一般后半段拟合较好,坐标xB=2/3×(实测数据最大值-实测数据最小值)+实测数据最小值,然后在xB左右各取0.25,计算该范围内实测点纵坐标的均值yB,得到控制点B(xB,yB)。
(4)由控制点A(辅助控制点B)带动实测数据点平移m次,得到离差平方和最小值fA(fB),以控制点A带动实测数据平移为例,如图3所示,其过程为:
1)k=1,即控制点A带动实测数据第一次平移,所述k同时是第k个理论点,k的取值范围为[1,m];
2)控制点A及n个实测数据点均平移所述为控制点A至第k个理论点的绝对距离,即控制点至理论点的实际距离,与坐标无关,作第k次平移;
3)(k=1、2…m;i=1、2…n),所述是第i个数据点和m个理论点距离的最小值(k=1、2…m),即第i个实测数据点至理论点的最短斜距,与求第i个实测数据点到标准曲线的最小距离等效,因为求离差平方和的目的在于找到实测曲线距离标准曲线最近的位置,所以该等效基本不影响自动配线的精度,而且不需要求点到曲线的距离,便于实现。其中为第i个数据点和第k个理论点的绝对距离,与坐标无关;
4)(i=1、2…n),所述fA(k)为第k次平移后n个实测数据点与标准曲线的离差平方和;
5)当k<m时,k=k+1,即平移次数小于m次时,均由控制点A带动实测数据向下一个理论点平移,求得相应的离差平方和,直到m次平移完成后,终止循环;每一次平移后求得的离差平方和fA(k)(k=1、2…m)均存储在数列中;
6)循环终止,fA=min[fA(k)](k=1、2…m),得到控制点A带动实测数据平移m次的最小离差平方和。
(5)min(fA,fB),即比较fA和fB的大小,取较小值;
(6)确定相应的平移量,得到最终配线位置,所述相应的平移量,如果fA≤fB,则相应的平移量为其中k为fA所对应的第k次平移;同理,如果fA>fB,则相应的平移量为其中k为fB所对应的第k次平移,根据平移量得到最终配线位置;
(7)选取匹配点,计算水文地质参数,承压含水层完整井定流量非稳定流:
式中:W(u)—井函数;Q—抽水井流量(m3/h);T—导水系数(m2/d);μ*—含水层的贮水系数;s—抽水影响范围内,任意一点任一时刻的水位降深(m);t—自抽水开始到计算时刻的时间(min);r—计算点到抽水井距离(m)。
在双对数坐标纸下,W(u)-1/u和s-t/r2的曲线形态基本一致,当配线完成后,任取一匹配点(在曲线上或曲线外均可),记下匹配点的对应坐标值:W(u),1/u,s和t/r2(或t),代入(1)、(2),分别计算有关参数:
本实施例中一种基于优化控制点确定水文地质参数的自动配线法的具体步骤如下:
(1)承压含水层非稳定流抽水试验
已知某一承压含水层抽水试验,抽水井号14,抽水井稳定流量为60m3/h,其中四个观测孔,观测孔2距离抽水井43m,观测孔15距离抽水井140m,观测孔10距离抽水井510m,观测孔1距离抽水井780m,累计抽水1185min,给出各个观测孔降深随时间变化的试验资料,如表1所示:
表1观测孔试验数据
(2)采用Excel表格以规定格式记录试验数据,即抽水时间除以观测孔与抽水井距离的平方(min/㎡)记录在表格的第一列,观测孔的降深(m)记录在表格第二列。
(3)通过一种基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合,最终选取匹配点计算导水系数、贮水系数等水文地质参数:
本实施例中采用VB程序实现自动配线法,首先,输入抽水井稳定流量60m3/h并读取Excel表格中的实测数据,进入绘图,自动生成配线法标准曲线,本实施例中标准曲线采用157个点理论点,因为157个理论点之后,标准曲线在横坐标达到108,早已趋于平稳,其后的标准曲线予以省略,不影响计算结果,同时可以减小后台计算量;然后单击“配线”,完成自动配线,实现标准曲线与实测数据曲线的拟合,如图5所示;选取任意一个拟合结果较好的匹配点,计算得到导水系数T=227.2122m2/d,贮水系数μ*=1.4251×10-4,计算速度很快,仅需2-3s,若采用阵列的方法,则至少需要1min才能完成配线。
(4)采用人工配线法计算得到的导水系数为212.3m2/d,贮水系数为1.47×10-4,计算结果与其比较精度满足工程要求。
Claims (9)
1.一种基于优化控制点确定水文地质参数的自动配线方法,其特征在于:包括如下步骤:
A、进行承压含水层非稳定流抽水试验,得试验资料;
B、基于非稳定流承压含水层抽水试验资料,绘制实测数据曲线与配线法标准曲线,给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合;
C、选取匹配点计算水文地质参数;
所述步骤B包括如下步骤:
B1、绘制配线法标准曲线,所述标准曲线由m个理论点连线而成;
B2、录入实测数据,所述实测数据为n个,由程序直接读取Excel表格录入;
B3、确定控制点坐标A(xA,yA),所述控制点A取n个实测数据的均值作为横坐标xA,然后在xA左右各取0.25,计算该范围内实测点纵坐标的均值yA,得到控制点A(xA,yA);由控制点A带动实测数据平移,每个理论点与控制点A在横、纵坐标方向有绝对距离和所述绝对距离即控制点至理论点的实际距离,与坐标无关,将控制点和所有实测数据点按平移量移动至标准曲线附近,计算每次的离差平方和fA(k)(k=1、2…m),确定fA的最小值,所述k表示第k次平移,同时也表示第k个理论点,k的取值范围为[1,m];
B4、取辅助控制点B(xB,yB),xB=2/3×(实测数据最大值-实测数据最小值)+实测数据最小值,然后在xB左右各取0.25,计算该范围内实测点纵坐标的均值yB,得到辅助控制点B(xB,yB),用与步骤B3同样的方法寻找离差平方和fB的最小值,取fA、fB中较小值,按对应的平移量移动实测数据,即可完成配线。
2.如权利要求1所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:所述步骤B为:采用控制点和辅助控制点带动所有实测数据点沿标准曲线移动,寻找离差平方和最小值,实现自动配线;所述步骤C为:根据配线结果计算水文地质参数。
3.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:所述步骤B3中,第k次平移的离差平方和fA(k),即每次移动后,计算第i个 实测数据点与标准曲线上m个理论点的距离(k=1、2…m),所述为第i个数据点和第k个理论点的绝对距离,与坐标无关,找出(k=1、2…m)中的最小值,该值是第i个实测数据点至理论点的最短斜距,与求第i个实测数据点到标准曲线的最小距离等效,所有实测数据点用同样的方法确定到标准曲线的绝对距离,平方求和就可以得到离差平方和fA(k);步骤B4中,同样的方法,得到离差平方和fB(k)。
4.如权利要求3所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:离差平方和最小值的确定包括如下步骤:
1)k=1,即控制点A带动实测数据第一次平移,所述k同时是第k个理论点,k的取值范围为[1,m];
2)控制点A及n个实测数据点均平移所述为控制点A至第k个理论点的绝对距离,即控制点至理论点的实际距离,与坐标无关,作第k次平移;
3)(k=1、2…m;i=1、2…n),所述是第i个数据点和m个理论点距离的最小值(k=1、2…m),即第i个实测数据点至理论点的最短斜距,与求第i个实测数据点到标准曲线的最小距离等效,其中为第i个数据点和第k个理论点的绝对距离,与坐标无关;
4)(i=1、2…n),所述fA(k)为第k次平移后n个实测数据点与标准曲线的离差平方和;
5)当k<m时,k=k+1,即平移次数小于m次时,均由控制点A带动实测数据向下一个理论点平移,求得相应的离差平方和,直到m次平移完成后,终止循环;每一次平移后求得的离差平方和fA(k)(k=1、2…m)均存储在数列中;
6)循环终止,fA=min[fA(k)](k=1、2…m),得到控制点A带动实测数据平移m次的最小离差平方和。
5.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:步骤B4中按对应的平移量移动实测数据,即可完成配线,所述对应的平移量,如果fA≤fB,则对应的平移量为其中k为fA所对应的第k次平移;如果fA>fB,则对应的平移量为其中k为fB所对应的第k次平移,根据平移量得 到最终配线位置。
6.如权利要求1所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:所述步骤C中的水文地质参数包括导水系数和贮水系数。
7.根据权利要求6所述一种基于优化控制点确定水文地质参数的自动配线方法,其特征在于:导水系数和贮水系数的按如下方法计算:
在双对数坐标纸下,W(u),-1/u和s,-t/r2的曲线形态基本一致,当配线完成后,在曲线上任取一匹配点,记下匹配点的对应坐标值:W(u),1/u,s和t/r2或t,代入(1)、(2),分别计算有关参数:
式中:W(u)—井函数;Q—抽水井流量,单位为m3/h;T—导水系数,单位为m2/d;μ*—含水层的贮水系数;s—抽水影响范围内,任意一点任一时刻的水位降深,单位为m;t—自抽水开始到计算时刻的时间,单位为min;r—计算点到抽水井距离,单位为m。
8.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:所述步骤A抽水试验为在抽水孔周围配置若干个观测孔,抽水过程中,保持抽水量固定而观测地下水位变化的抽水试验。
9.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法,其特征在于:所述步骤A抽水试验资料包括抽水井流量、观测孔与抽水井的距离及观测孔降深随时间的变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310048344.5A CN103149600B (zh) | 2013-02-06 | 2013-02-06 | 一种基于优化控制点确定水文地质参数的自动配线方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310048344.5A CN103149600B (zh) | 2013-02-06 | 2013-02-06 | 一种基于优化控制点确定水文地质参数的自动配线方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103149600A CN103149600A (zh) | 2013-06-12 |
CN103149600B true CN103149600B (zh) | 2015-06-17 |
Family
ID=48547781
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310048344.5A Expired - Fee Related CN103149600B (zh) | 2013-02-06 | 2013-02-06 | 一种基于优化控制点确定水文地质参数的自动配线方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103149600B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104933117B (zh) * | 2015-06-05 | 2017-02-22 | 中国地质大学(武汉) | 一种抽水试验数据的处理方法 |
CN106501157B (zh) * | 2016-12-13 | 2019-07-12 | 河海大学 | 现场确定外管弱透水层水文地质参数的内管降深双管法 |
CN107369191B (zh) * | 2017-08-15 | 2021-01-15 | 国网湖南省电力有限公司 | 电网气象灾害预测色斑图修正方法、系统及装置 |
CN108596429A (zh) * | 2018-03-14 | 2018-09-28 | 中国石油天然气股份有限公司 | 证实已开发储量年度修正值的计算方法和装置 |
CN108959762A (zh) * | 2018-06-29 | 2018-12-07 | 江铃汽车股份有限公司 | 一种高效衬套拟合方法 |
CN110865008A (zh) * | 2019-10-28 | 2020-03-06 | 河海大学 | 基于有限尺度下圆形定水头边界非稳定流抽水的水文地质参数测定方法 |
CN112698002A (zh) * | 2020-12-29 | 2021-04-23 | 张勇 | 非稳定流抽水试验水文地质参数智能计算方法 |
CN113050190A (zh) * | 2021-03-03 | 2021-06-29 | 河北益坤岩土工程新技术有限公司 | 直线边界非稳定流抽水试验水文地质参数智能计算方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809642A (zh) * | 2012-08-03 | 2012-12-05 | 河海大学 | 一种确定弱透水层水文地质参数的方法 |
-
2013
- 2013-02-06 CN CN201310048344.5A patent/CN103149600B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809642A (zh) * | 2012-08-03 | 2012-12-05 | 河海大学 | 一种确定弱透水层水文地质参数的方法 |
Non-Patent Citations (3)
Title |
---|
《Theis解析模型的数值解及其在MatLab中的实现》;李林子等;《江苏地质》;20070331;第31卷(第3期);242-245 * |
《利用配线法求水文地质参数计算机程序简介》;朱国祥等;《工程勘查》;19990331(第3期);39-41 * |
《求水文地质参数的计算机配线法》;邹正盛等;《工程勘查》;20010630(第6期);31-32 * |
Also Published As
Publication number | Publication date |
---|---|
CN103149600A (zh) | 2013-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103149600B (zh) | 一种基于优化控制点确定水文地质参数的自动配线方法 | |
CN112417573B (zh) | 基于ga-lssvm与nsga-ⅱ盾构下穿既有隧道施工多目标优化的方法 | |
CN105389467B (zh) | 一种获得井间连通关系的方法及装置 | |
CN107423518B (zh) | 利用Civil3D与Revit结合的综合管廊全线BIM协同设计方法 | |
CN105239976A (zh) | 考虑水驱油藏渗流参数动态变化的数值模拟方法 | |
CN109948272A (zh) | 基于井间连通性的调堵动态预测方法和系统 | |
CN103886135B (zh) | 基于二维非恒定流数值模型的电力工程选址方法 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN111625922A (zh) | 一种基于机器学习代理模型的大规模油藏注采优化方法 | |
CN105631529A (zh) | 一种边水气藏见水时间预测方法 | |
CN104314547A (zh) | Hafd水平井压裂酸化分段优化方法及系统 | |
CN104598996A (zh) | 基于最小二乘支持向量机的施工引起地表变形的预测方法 | |
CN103020709B (zh) | 基于人工蜂群和量子粒子群算法的一维水质模型参数率定方法 | |
CN109839493A (zh) | 地下工程岩体质量评价方法、装置、存储介质及电子设备 | |
CN103530700B (zh) | 城区配电网饱和负荷综合预测方法 | |
CN104732091A (zh) | 基于自然选择蚁群算法的元胞自动机河床演变预测方法 | |
CN102692491B (zh) | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 | |
CN115860197A (zh) | 一种基于数据驱动的煤层气产量预测方法及系统 | |
CN106650028A (zh) | 一种基于敏捷卫星设计参数的优化方法及系统 | |
CN104727789A (zh) | 中高渗砂岩油藏水驱波及系数及过水倍数动态描述方法 | |
CN105740501A (zh) | 一种基于定常数值计算的泵装置水力稳定性分析方法 | |
CN103886129A (zh) | 将测井数据离散到储层网格模型的方法和装置 | |
CN106846481B (zh) | 一种地质剖面图的生成方法 | |
CN116187001A (zh) | 一种地浸采铀有效溶浸范围的自动提取方法 | |
Lian et al. | Evaluation and applicability study on prediction methods of water inflow in mines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 Termination date: 20190206 |
|
CF01 | Termination of patent right due to non-payment of annual fee |