智能优化算法应用：基于麻雀搜索算法的水文地质参数优化 -附代码

智能优化算法应用：基于麻雀搜索算法的水文地质参数优化

文章目录

• 智能优化算法应用：基于麻雀搜索算法的水文地质参数优化

• 1.麻雀搜索算法
• 2.目标函数建模
• 3.实验测试
• 4.参考文献
• 5.Matlab代码

摘要：为提高水文地质参数求解精度，拓展水文地质参数优化方法，本文基于泰斯(Theis)基本公式，利用麻雀搜索算法优化求解 Theis 公式导水系数和储水系数，旨在验证麻雀算法用于水文地质参数优化求解的可行性和有效性。

1.麻雀搜索算法

 

2.目标函数建模

在抽水试验中，对于符合泰斯(Theis)假设的含水层系统，观测井水位降深可表示为:
$s = \frac{QW(u)}{4\pi T} \tag{1}$

$W(u) = ln\frac{2.25Tt}{r^2S} \tag{2}$

式中 $s$ 为水位降深（m）; Q为抽水井流量（$m^3/d$）; $T$为导水系数($m^2/d$); $W(u)$ 为泰斯井函数，$u$为参变量，在满足 $u ≤ 0. 01$ 或$u ≤ 0. 05$ 时，$W(u) = ln\frac{2.25Tt}{r^2S}$成立; t为自抽水开始到计算时刻的时间，$d$; $r$为观测孔距抽水井的距离(m); $S$

基于式(1)、(2)构建目标适应度函数:
$\begin{cases} minf(T,S)=\sum_{i=1}^N(s-s$
式中$s$为实测水位降深，m; $s$为优化求解水位降深，m; i 为抽水序号或观测井号，i =1，2，…，N ; T 为待优化参数导水系数;S为储水系数。

3.实验测试

实例 1:在某一地区承压含水层进行多孔抽水试验，抽水井稳定流量为 60 m 3 /h，14 号井为完整抽水井，观 2 孔、观 15 孔、观 16 孔、观 10 孔及观 9 孔为观测孔。观孔 15 距抽水井距离 r 为 125 m，抽水过程中观测井中不同时刻的降深见表 1。本文以观
孔 15 抽水实验观测资料为例，利用 Gold － SA 算法优化求解导水系数 T 和储水系数 S。

表1 观孔 15 抽水试验数据

序号	时间t/min	降深s/m
1	10	0.16
2	20	0.48
3	30	0.54
4	40	0.65
5	60	0.75
6	80	1
7	100	1.12
8	120	1.22
9	150	1.36
10	210	1.55
11	270	1.7
12	330	1.83
13	400	1.89
14	450	1.98
15	645	2.17
16	870	2.38
17	990	2.46
18	1185	2.54

参数设置。设置导水系数 T 和储水系数 S的搜索范围均为［0，1000］;

麻雀搜索实验结果如下：

从上图曲线可以看出，优化得到的T，S模拟的降深与实际降深非常接近，表明该种方法具有一定的可行性。

SSA 获得的最佳T，S : 9.1354 0.013582
SSA 获得的最佳适应度值 : 0.0061864

4.参考文献

[1] 周有荣, 李娜, 周发辉. 黄金正弦算法在水文地质参数优化中的应用[J]. 人民珠江, 2020, 041(006):117-120,128.

5.Matlab代码

[1] 周有荣, 李娜, 周发辉. 黄金正弦算法在水文地质参数优化中的应用[J]. 人民珠江, 2020, 041(006):117-120,128.
个人资料介绍