用非线性回归方法研究地下水动态规律.pdfVIP

用非线性回归方法研究地下水动态规律 张子贤 (徐州建筑职业技术学院，江苏徐州 221008) 摘要采用非线性回归方法研究地下水动态规律，给出了应用最小二乘法原理和双．最割 线迭代法求解一类纯非线性回归模型的方法。实例表明，使用本文方法简便易行，效果良 好。时实际资料分析计算表明，阿雏扬诺夫潜水蒸发公式中的指数n存在小于1的情况。 关键词 地下水动态规律；纯非线性模型；最小二乘法；双点割线迭代法}拟合精度 对地下水动态规律的分析与研究，是地下水资源研究与计算的重要组成部分。由于地 下水动态的影响因素比较复杂，且一般具有随机性，因此目前对地下水动态规律的研究广 泛采用数理统计法，涉及的数学模型也常常是纯非线性的。对于纯非线性回归模型，可用 将反映地下水动态的非线性规律简化为线性的，然后通过线性回归求解。因而一定程度上 降低了地下水动态资料的拟合效果和所求得的地下水含水层参数的精度，从而影响水资源 计算与评价的成果。本文结合地下水动态分析研究的有关问题，针对仅一个参数以非线性 形式出现的纯非线性回归模型，采用一种简便易行的计算方法求解模型的回归系数。实例 表明，效果良好。 1基本原理 以三个未知参数口，b，c(当然可以多于3个)为例，设纯非线性模型为 ，一f(x，口，btc) (1) 设式(1)中仅参数f以非线性形式出现。利用m组观察值(z，，，．)．i一1～m，根 据最小二乘法，求n，6，c使残差平方和 Q(a，6，c)一∑[弘一f(x朋，b，c)]2 (2) |2】 为最小。则 一 2 [，，一m一，缸)](一誓) (3) 一 2 b—m刖，¨)](一菇) (4) 硷瓦啦丽硷瓦 一 2 ，∑㈦。∑㈦。∑㈨ b-f(如柚∽](一誓) (5) 彳 52 由于参数n，6以线性形式出现，联解式(3)、(4)，可求得口，b，但它们均依赖于 程。利用求解非线性方程根的数值解法，可求得c，进而可求得a，b。 2应用实傍0 2．1潜水位上升和次降雨量关系的确定 针对开封惠北试区一地下永动态观测点[“，在地下水非开采条件下，雨前地下水埋深 为3m以内的降雨和潜水位上升关系的观测资料，经分析次降雨量与潜水位上升值之间有 如下关系式： △H=Ⅱ’(p+丸一y—A)6 (6) 式中：△H为次降雨引起的潜水位上升高度。mm；P为次降雨量，mill；丸为前期影响 雨量，mmY为次降雨地表径流量ram!。7，b为待定常数；p。为能引起地下水位上升 的最小雨量，也称临界雨量，mm，待求。 Y=口+bln(x—c) (7) 可见式(7)为纯非线性模型，但只有参数f以非线性形式出现，采用本文方法求式(7) 中的拟合参数，使 Q(a，b，c)一∑[，，～n—bln(x：一f)]2 l一1 为最小。则 筹一2耋b—n-6In(置刊]：o 鬻一2塞(y，一n_6ln(”c))l。(Xi--C)：o 鬻一z耋‰一一瞰一m(暑)c川一。 引入如下记号并化简，得式(8)、(9)、(10) L，一∑弘，L。=∑ln(x，一c)，L。一∑yiln(x；一c) L3一蚤k2(x,--chLt 2蚤毒b， i‘…。， ^ L