专利名称：一种基于优化控制点确定水文地质参数的自动配线方法
技术领域：
本发明涉及一种确定承压含水层水文地质参数的方法，尤其是涉及一种基于优化控制点实现自动配线确定承压含水层导水系数、贮水系数等水文地质参数的方法。
背景技术：
近年来，随着我国人口的迅速增长和经济的快速发展，城市建设规模和速度前所未有，城市建设用地日益紧张，地下空间的开发和利用也日益紧迫，北京和沿海一些大城市地下轨道交通和车站、大型地下商场、地下停车场等大型深基坑工程层出不穷，这些深基坑的施工成败与能否较好的控制地下水关系密切，为降低深基坑工程的设计与施工风险，有必要通过抽水试验确定各含水层，特别是承压含水层的水文地质参数。此外，抽水试验也是煤炭资源等地质勘探的重要手段，其目的是研究含水层重要水文地质特征，取得含水层水文地质参数，评价含水层的富水性，并为预计矿井涌水量与对地下水综合利用的评价提供资料。因此，承压含水层水文地质参数的确定对揭示抽取地下水引起地下水变化和运动的基本规律有着重要的意义，而且为开采地下水和石油以及地下工程建设等生产活动提供理论依据和技术保障。非稳定流承压完整井求取水文地质参数的方法主要有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在横、纵坐标方向有绝对距离Axi4和Ag ,
所述绝对距离即控制点至理论点的实际距离，与坐标无关，将控制点和所有实测数据点按
平移量(M〗，Ay).移动至标准曲线附近，计算每次的离差平方和fA(k) (k=l、2…m)，确定fA
的最小值，所述k表示第k次平移，同时也表示第k个理论点，k的取值范围为[l，m]；B4、取辅助控制点B(xb，yB),xB=2/3X (实测数据最大值_实测数据最小值)+实测数据最小值，然后在xB左右各取0.25，计算该范围内实测点纵坐标的均值yB，得到控制点B(xb, yB)，用与步骤B3同样的方法寻找离差平方和fB的最小值，取fA、fB中较小值，按对应的平移量移动实测数据，即可完成配线。上述步骤B3中，标准曲线(理论点组成)和实测曲线(实测数据点组成)的双对数坐标值可能不同。所述步骤B3中，第k次平移的离差平方和fA(k)，即每次移动后，计算第i个实测数据点与标准曲线上m个理论点的距离式(k=l、2-m)，所述4为第i个数据点和第k个理
论点的绝对距离，与坐标无关，找出4 (k=l、2…m)中的最小值，该值是第i个实测数据点至
理论点的最短斜距，与求第i个实测数据点到标准曲线的最小距离等效，因为求离差平方和的目的在于找到实测曲线距离标准曲线最近的位置，所以该等效基本不影响自动配线的精度，而且不需要求点到曲线的距离，便于实现。所有实测数据点用同样的方法确定到标准曲线的绝对距离，平方求和就可以得到离差平方和fA(k);步骤B4中，同样的方法，得到离差平方和fB(k)。为了提高精度，上述离差平方和最小值的确定包括如下步骤:l)k=l,即控制点A带动实测数据第一次平移，所述k同时是第k个理论点，k的取值范围为[1，m]；2)控制点A及n个实测数据点均平移所述为控制点A至第k个理论点的绝对距离，即控制点至理论点的实际距离，与坐标无关，作第k次平移；3) =minjy丨](k=l、2…m; i=l、2*“n)，所述是第i个数据点和m个理论
点距离的最小值(k=U2-m),即第i个实测数据点至理论点的最短斜距，与求第i
个实测数据点到标准曲线的最小距离等效，其中4力第i个数据点和第k个理论点的绝对距离，与坐标无关；4)/400 = 5_[(<皿)2](1=1、2-11)，所述4(10为第k次平移后n个实测数据点与标准曲线的离差平方和；5)当k〈m时，k=k+l,即平移次数小于m次时,均由控制点A带动实测数据向下一个理论点平移，求得相应的离差平方和，直到m次平移完成后，终止循环；每一次平移后求得的离差平方和fA(k) (k=l、2…m)均存储在数列中；6)循环终止,fA=min[fA(k) ] (1^=1、2夂111),得到控制点A带动实测数据平移m次的
最小尚差平方和。上述步骤B4中按对应的平移量移动实测数据，即可完成配线，所述相应的平移量，如果fA彡fB,则相应的平移量为(Ax14^i)，其中k为fA所对应的第k次平移；如果
fA>fB，则相应的平移量为，其中k为fB所对应的第k次平移，根据平移量得到最终配线位置。上述步骤C中的水文地质参数一般包括导水系数和贮水系数。导水系数和贮水系数的按如下方法计算:
权利要求
1.一种基于优化控制点确定水文地质参数的自动配线方法，其特征在于:包括如下步骤: A、进行承压含水层非稳定流抽水试验，得试验资料； B、基于非稳定流承压含水层抽水试验资料，绘制实测数据曲线与配线法标准曲线，给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合； C、选取匹配点计算水文地质参数。
2.如权利要求1所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤B为:采用控制点和辅助控制点带动所有实测数据点沿标准曲线移动，寻找离差平方和最小值，实现自动配线；所述步骤C为:根据配线结果计算水文地质参数。
3.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤B包括如下步骤: B1、绘制配线法标准曲线，所述标准曲线由m个理论点连线而成； B2、录入实测数据，所述实测数据为η个，由程序直接读取Excel表格录入； B3、确定控制点坐标A(xa，yA)，所述控制点A取η个实测数据的均值作为横坐标χΑ，然后在xA左右各取0.25，计算该范围内实测点纵坐标的均值yA，得到控制点A (xA, yA)；由控制点A带动实测数据平移，每个理论点与控制点A在横、纵坐标方向有绝对距离Δχ:和,所述绝对距离即控制点至理论点的实际距离，与坐标无关，将控制点和所有实测数据点按平 移量移动至标准曲线附近，计算每次的离差平方和fA(k) α=1、2...πι)，确定fA的最小值，所述k表示第k次平移，同时也表示第k个理论点，k的取值范围为[l，m]； B4、取辅助控制点B(xB，yB)，xb=2/3X (实测数据最大值-实测数据最小值)+实测数据最小值，然后在xB左右各取0.25，计算该范围内实测点纵坐标的均值yB，得到控制点B(xb, yB)，用与步骤B3同样的方法寻找离差平方和fB的最小值，取fA、fB中较小值，按对应的平移量移动实测数据，即可完成配线。
4.如权利要求3所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤B3中，第k次平移的离差平方和fA (k)，即每次移动后，计算第i个实测数据点与标准曲线上m个理论点的距离c/i(k=l、2…m)，所述<为第i个数据点和第k个理论点的绝对距离，与坐标无关，找出^ U=l、2…m)中的最小值，该值是第i个实测数据点至理论点的最短斜距，与求第i个实测数据点到标准曲线的最小距离等效，所有实测数据点用同样的方法确定到标准曲线的绝对距离，平方求和就可以得到离差平方和fA(k);步骤B4中，同样的方法，得到离差平方和fB(k)。
5.如权利要求4所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:离差平方和最小值的确定包括如下步骤: 1)k=l,即控制点A带动实测数据第一次平移，所述k同时是第k个理论点，k的取值范围为[I，m]； 2)控制点A及η个实测数据点均平移，所述为控制点A至第k个理论点的绝对距离，即控制点至理论点的实际距离，与坐标无关，作第k次平移；3)d'mm =min[^] (k=l>2---m; i=l>2---n),所述4是第i个数据点和m个理论点距离的最小值(k=U2-m),即第i个实测数据点至理论点的最短斜距，与求第i个实测数据点到标准曲线的最小距离等效，其中4为第i个数据点和第k个理论点的绝对距离，与坐标无关； 4)臟(i=l、2…n)，所述fA(k)为第k次平移后n个实测数据点与标准曲线的离差平方和； 5)当k〈m时，k=k+l,即平移次数小于m次时,均由控制点A带动实测数据向下一个理论点平移，求得相应的离差平方和，直到m次平移完成后，终止循环；每一次平移后求得的离差平方和fA(k) (k=l、2…m)均存储在数列中； 6)循环终止，fA=min[fA(k)](1^=1、2夂111),得到控制点A带动实测数据平移m次的最小尚差平方和。
6.如权利要求3所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:步骤B4中按对应的平移量移动实测数据，即可完成配线，所述相应的平移量，如果fA ( fB,贝U相应的平移量为其中k为fA所对应的第k次平移；如果fA>fB，则相应的平移量为(A4，，其中k为fB所对应的第k次平移，根据平移量得到最终配线位置。
7.如权利要求1所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤C中的水文地质参数包括导水系数和贮水系数。
8.根据权利要求7所述一种基于优化控制点确定水文地质参数的自动配线方法，其特征在于:导水系数和贮水系数的按如下方法计算: S=-Q-W{u)⑴ AnT 2 * r a,、 (J) ATt在双对数坐标纸下，Ku)，-1/U和S，-t/r2的曲线形态基本一致，当配线完成后，在曲线上或曲线任取一匹配点，，记下匹配点的对应坐标值:W(u)，1/u, s和t/r2或t，代入(I)、(2)，分别计算有关参数: ()* t(4)jl—口y I Ir" _ Ji —式中:W(u) —井函数；Q—抽水井流量，单位为m3/h ;T—导水系数，单位为m2/d —含水层的贮水系数；s—抽水影响范围内，任意一点任一时刻的水位降深，单位为m;t—自抽水开始到计算时刻的时间，单位为min ；r一计算点到抽水井距离,单位为m。
9.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤A抽水试验为在抽水孔周围配置若干个观测孔，抽水过程中，保持抽水量固定而观测地下水位变化的抽水试验。
10.如权利要求1或2所述的基于优化控制点确定水文地质参数的自动配线方法，其特征在于:所述步骤A抽水试验资料包括抽水井流量、观测孔与抽水井的距离及观测孔降深随时间 的变化。
全文摘要
本发明公开了一种基于优化控制点确定水文地质参数的自动配线方法，包括如下步骤进行承压含水层非稳定流抽水试验，得试验资料；基于非稳定流承压含水层抽水试验资料，绘制实测数据曲线与配线法标准曲线，给出通过基于优化控制点的自动配线法实现实测曲线与标准曲线的拟合；选取匹配点计算水文地质参数。本发明提出的计算机自动配线的新方法，人机交互简单，接近传统配线法的基本操作模式；通过控制点带动所有实测点，可以大大减小沿标准曲线移动的次数，相当于寻找了一条最优路径移动实测数据点，在保证精度的同时，缩短了后台操作时间，提高了工作效率。
文档编号G01V9/00GK103149600SQ20131004834
公开日2013年6月12日 申请日期2013年2月6日 优先权日2013年2月6日
发明者赵燕容, 朱旭芬, 张静, 李志利, 李伟, 武立军, 王璐佳 申请人:河海大学