基于Matlab求解稳定流的水文地质参数.pdf

《基于Matlab求解稳定流的水文地质参数.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、江西煤炭科技2O年第3期JIANGXICOALSCIENCE＆TECHNOLOGYNO．32014基于Matlab求解稳定流的水文地质参数吴家杰，史丽娜(江西省煤田地质局一九五地质队，江西南昌330001)摘要：相比传统的迭代法，通过利用MATLAB软件fzero函数，求解抽水试验水文地质参数更为快捷简便。关键词：稳定流；水文地质参数；MATLAB；渗透系数中图分类号：P641．73文献标识码：B文章编号：1006—2572(2014)03—0124—03SolutiontoHydrogeologicalParametersofSt

2、eadyFlowBasedonMatlabWuJiajie．ShiLina(No．195GeologicalTeam，JiangxiCoalFieldGeologyBureau，Nanchang，Jiangxi330001)Abstract：Comparedwithtraditionaliterativemethod，itismorequickandeasytOsolvehydrogeologyparametersinpumpingtestbyusingMATLABsoftfzerofunction．Keywords：steadyf

3、low；hydrogeologyparameters；MATLAB；permeabilitycoefficient1概述梅山矿区(南区)范围内分布的主要地层由老至新为：二叠系中统茅口组(Pm)、三叠系上统安源组紫家冲段导水系数(T)和渗透系数(K)是进行地下资源论证、(Ta)、三叠系上统安源组龙王寨段(T。a)、三叠系上统安地下水资源评价和地下水数值模拟的水文地质参数。抽源组高家段(Ta。)、侏罗系下统水北组(Jsh)、白垩系南雄水试验法是确定含水层参数的主要方法。组(K。n)、第四系(Q)。主要含水层为白垩系砾岩裂隙岩MATLAB

4、是矩阵实验室(MatrixLaboratory)的简称，溶含水层、安源组紫家冲段底部砾岩含水层，茅口组灰岩是美国MathWorks公司出品的商业数学软件，可用于算岩溶含水层，其它可视为隔水层或相对隔水层。法开发、数据可视化、数据分析以及数值计算。采用2．2抽水试验情况MATLAB软件，不仅可以快速进行大量数据的科学计算，梅山矿区(南区)在详查阶段共进行了二孔二层次的还可以克服人工配线过程中产生的人工误差。与其它软稳定流抽水试验工作，两个水文地质钻孔施工过程中均严件相比，MATLAB操作比较容易，编程相对较简单。本文格按照水文地质钻探

5、施工要求进行施工，对含水层采用清依据梅山矿区抽水试验资料，采用MATLAB程序来计算水钻进，并采用套管隔水，止水效果好，钻探质量等级为甲承压含水层导水系数T和渗透系数K，取得了较好的结级；抽水试验前统一进行含水层盐化测井，测井成果准确果。可靠，测井质量等级为甲级；抽水试验严格遵照《煤矿床水2抽水试验文地质、工程地质及环境地质勘查评价标准》进行，抽水质2．1水文地质概况量均达合格以上，钻孔质量等级综合评定均为甲级。·124·MATLAB求解渗透系数后，依据公式R一104结语SwK，求出各孔各次降深对应的影响半径R，水1孔的3以往传统求

6、解抽水试验水文地质参数(K、R)，一般用次降深的影响半径R1、R2、R3分别为537m、720m、870迭代法，在实际过程应用偏繁琐，而通过装布依自流完整井m；水2孔的3次降深的R1、R2、R3分别为3722m、3065公式和影响半径经验公式，化为一元方程，利用MATLABm、2653m。软件fzero函数求解，优点在于计算过程简单，产生误差的2个抽水钻孔各个水文地质参数见表3。概率小，精度高。而传统迭代法先任意取值，一直迭代，不表3某矿区抽水试验水文地质参数一览断重复计算，直至数值接近目标值，实际生产过程容易出差，且难以发现。因而

7、利用MATLAB函数求解水文地质参数在实际生产过程运用比传统迭代法更为方便、可靠。参考文献：(1]李俊庭，王愈吉．地下水动力学(M]．北京：地质出版社，1986．3．3与传统迭代法对比[2]石志远，许模．抽水试验资料的遗传算法反演水文地质参MATLAB函数求解水文地质参数与传统迭代法求解数[J]．四川地质学报，2011，31(4)：466—469，480．水文地质参数结果见表4，从表中可以看出，水1孔两种方[33刘卫国．MATLAB升序设计与应用(第二版)[M]．北京：高法求解水文地质参数数值结果误差不大，而水2孔两种方等教育出版社

8、，2006．法求解水文地质参数数值结果产生一定误差，可能导致估[4]徐东艳，孟晓刚．MATLAB函数库查询辞典[M]．北京：中国算涌水量产生较大偏差。铁道出版社，2005．表4MATLAB函数与传统迭代法求解结果／(1MATLAB传统