应用并行PEST算法优化地下水模型参数
更新时间:2023-06-01 07:27:01 阅读量: 实用文档 文档下载
- 图像处理并行算法与应用推荐度:
- 相关推荐
工程地质学报1004—9665/2010/18(1)一0140—05
应用并行PEST算法优化地下水模型参数
董艳辉①
李国敏①郭永海②
徐海珍①
(①中国科学院地质与地球物理研究所北京100029)(②核7-_业北京地质研究院北京100029)
摘要基于列文伯格一马夸尔特(Levenberg—Marquardt)算法的PEsT参数优化程序具有寻优速度快、健壮性好的优点,在地下水模型参数优化研究中有许多成功的应用实例。但是,对于大尺度、高精度和高复杂性的大规模地下水模拟,使用PEST进行参数优化需要大量的计算时间,优化效率较低。本文应用OpenMP并行编程方法对PEST算法进行了并行化,使之可以在共享存储并行计算机上进行参数优化的并行计算。并将此方法应用于甘肃北山区域地下水模型的参数优化中,并行实验表明,使用并行化的PEST可以将地下水模型参数优化效率提高3.7倍。关键词PEST
MODFLOW参数优化甘肃北山
文献标识码:A
并行计算
中图分类号:P641.73
oPTD江IZATIONOF
MODELPARAMETERS
FOR
GRoUNDWATERFLoW
USINGPARALLELIZEDPEST
METHOD
DONGYanhui①LIGuomin①
GUOYonghai②XUHaizhen①
Academy
((I)lnstituteofGeologyandGeophysics,ChineseofSciences,Beijing100029)
100029)
(②既彬昭Research
Abstract
Institute
ofUraniumGeology,Beijing
Duringthelastdecades,anumberofmethodshavebeendevelopedfortheoptimizationofhydrologic
a
modelparameters.OnefrequentlyusedmethodisthePESTmethodwhichiszationalgorithm.As
a
Levenberg—Marquardt—basedoptimi—
nonlinearparameterestimator,thePESTmethod
can
call
existindependentlyofanyparticular
models,canestimateparametersand/orexcitations,andtypes
cover
a
carryout
variouspredictiveanalysistasks.Itsmodel
widerange.Itispreferabletocalibrate
a
a
MODFLOWmodelusingthePESTmethodinmanysitua—
tions.However,for
large—scale,high—resolution
or
high—complexitygroundwaterflowmodeling,theparameterop—
timizationrequiresmassivecomputingtime.Inthisstudy,theOpenMPprogrammingparadigmisusedtoachievethemodestparallelism
to
ona
shared memorycomputerforthePESTmethod.Application
can
oftheparallelPEST
method
Beishanarea,Gansuprovinceshowsthatthecomputingspeed
beincreasedupto3.7times.Thecomputing
speedisdefinedastheratiooftheserialcomputingtimeandtheparallelcomputingtime.TheparallelPESTpa-rameteroptimization
method
can
be
an
usefulandeffectivetoolforlargemodels.
KeywordsPEST,MODFLOW,Parameteroptimization,BeishanMountain,Parallelcomputing
收稿日期:2009—02—27;收到修改稿日|胡:2009-03—31.
基金项目:“中国科学院知识创新工程重要方向项目(Kzcx2一yw一126—04,Kzcx2一yw一116),水利部公益性行业科研专项
(200801020),环保公益性行业科研专项(200709055).
第一作者简介:董艳辉,水文地质专业.Email:lemondyh@mail.iggcas.∽.cn
万方数据
18(1)董艳辉等:应用并行PEST算法优化地下水模型参数
1
引言
为使得所建立的地下水模型可以较真实地反映研究区实际水文地质条件,确定模型中的水文地质参数是重要的研究工作。实际工作中,传统的方法
一般采取试估一校正法(试算法)来进行参数的确
定。这种人工迭代方法缺乏收敛判别标准,难以达到参数最优化且工作量很大。因此,将水文地质参数确定归结为求极值问题,使得误差评价函数达到最小的最优化方法得到广泛地应用,如最速下降法、二次规划、共轭梯度方法、高斯一牛顿法等¨.2J。而针对基于梯度搜索优化方法的缺点,蚁群算法、人工神经网络、模拟退火算法及遗传算法等也被引入参数优化研究中,并取得了较好地效果¨。J。
随着地下水研究的发展,研究尺度的扩大、精度的提高以及复杂性的增加,使得地下水模型的规模、水文地质参数分区数目也相应增大,导致自动优化的求解效率和计算精度的下降。例如大尺度的地下水模型,如果进行模型规模和参数分区时的简化,进行水文地质参数的优化时计算时间可以满足要求,但是精度会大大减低;而使用较大规模的模型和分区时,需要耗费计算时间,效率很低。
PESTL61参数优化算法是一种比较成熟、应用较多的地下水模型参数优化方法,该方法在美国DeathValley等区域地下水模型参数优化中取得了很好效果的ET,S]。但是对于大规模的地下水模型,同样存在优化计算时间过长的问题。本文研究针对此问题,使用OpenMP【91并行编程方法对PEST算法进行了并行化,使该算法可以在共享存储并行计算机上进行参数优化的并行计算。应用并行化的PEST算法,对甘肃北山区域地下水模型进行了水文地质参
数并行优化,并对参数优化结果和并行计算效果进
行了分析。2
PEST参数优化算法
PEST算法的核心是使用列文伯格一马夸尔特
算法对目标函数求解最小值,目标函数则是模型参
数计算值与实际观测值的差异函数,而模型的计算
值主要依赖于模型的参数∞J。假设模型参数存储
在矩阵菇中,而观测值存储在矩阵Y中,则菇与Y的
关系可以用下式表示:
Y=M(菇)
(1)
万方数据
141
戈为有n个(参数的个数)元素的矩阵,Y为n个
(观测值的个数)元素的矩阵。M是一个非线性函数,它将模型参数与观测值关联在一起。将式(1)
线性化可以得到下式:
多=Y+H(茹一髫)
(2)
茹和多分别是参数与模型结果的矩阵。H为m行
n列的雅可比偏导数矩阵,计算公式如下:
H[i,j]=aM(圣)[i]/ax[j]
(3)
i和_『分别是H的行列数。PEST在计算时对矩
阵曼循环更新,每次循环加入增量矩阵/2,:
『MI=(H^ToH女)。H。To(y一靠)
…
I;:=未i+配。
‘4,
式中,k为循环次数,矩阵上标“一”代表更新矩阵前的计算值,而“+”代表更新矩阵后的值,r为转置符号。0为m行/71,列的观测值权重矩阵,假设有两
类观测值时(如渗透系数和降雨入渗系数),矩阵0
的对角线元素为每一个观测在整个模型误差中的相对权重,非对角线元素均为零。由式(4)可以得到目标函数如下:
咖=[多一Y—H(金一茗)]70[多一Y—H(未一菇)]
(5)
算法计算过程概括如下:首先,选定初始参数矩
阵拓,然后进入整个模型的校正期;模型模拟结果存人矩阵菇中,相对应的观测值存入矩阵Y中;使用一阶泰勒展开式对凰进行数值求解;应用式(4)计算参数增量矩阵‰及更新的矩阵量÷;将;÷的值存
入参数矩阵;i中;循环整个算法直到达到收敛。
在进行参数估计计算的同时,可以进行参数敏感性分析。参数敏感性分析可以用来评价已知数据对于待优化的参数是否充足。研究中使用综合敏感
性来分析,综合敏感性代表了对于一个参数的所有敏感性。对于参数Z,其综合敏感性计算如下:
sf=(H70H)n1/2/,,m
(6)
H为雅可比矩阵,0为观测值权重矩阵,m为观
测值的个数。s。比较大的参数对于整个优化过程是比较容易的,而s。比较小时,说明该参数比较难或者
不可能被正确的优化。
3
PEST参数优化算法的并行方法
MODFLOW—ASP是使用PEST算法来进行
MODFLOWⅢ3地下水模型参数优化的程序。对于待
估的水文地质参数,PEST在其优化空间内循环调用
142
MODFLOW的计算结果,从而进行最优化组合的求解(图1a)。
PEST算法中对目标函数的极小值的计算是比
较快的,然而由于在每次求解模拟值时都要调用MODFLOW进行计算,使整个优化时间大大增加。
在整个优化时间中,MODFLOW计算所花的时间往往占整个时间的99%以上,并且随着模型的增大,所占比例也大大增加。因此将PEST调用MODF-Low计算的部分进行并行化,并在高性能计算机上运算将会大大缩短优化时间。
OpenMP是用于共享存储编程的一套可移植、可扩展的应用程序接口(Application
ProgramInter-
face),通过一系列编译指示、运行时函数库和环境变量来描述程序的并行特性,支持C、C++和Fortran(77,90和95)语言。OpenMP是共享存储平台的并行编程标准,为其建立了一套简单的编译指示。有时仅仅使用3、4种编译指示就可以得到显著的并行效果。OpenMP支持增量并行,对于一个串行程序,可以每次只对一部分代码并行化,最终将整个程序
并行化,这也是共享存储模型相对于消息传递模型
的优点。同时,不管是粗粒度还是细粒度,OpenMP都提供了很好的支持。
在研究中使用OpenMP对MODFLOW的PCG¨u算法进行了并行化,并与PEST进行连接(图lb)。经过并行化的MODFLOW—ASP可以在共享存储并行计算机或多核电脑上进行并行参数优化计算。
初始参数Il初始参数
...........[
坚茎至茎1
0DFLOw计算
ODFLOw计算
二二工二二
得到模拟值
得到模拟值
优化结果优化结果
b
图1
MODFLOW—ASP参数优化流程图
Fig.1
Flowchartof
MODFLOW—ASPprocesses
a.MODELOW—ASP流程;b.优化MODELOW—ASP流程
4
应用实例
北山地区位于甘肃省西部,北纬40。30’一42。
10’,东经96。30’一99。50’,总面积约20000km2。北万方数据
JournalofEngineeringGeology工程地质学报2010
山地区是我国高放废物地质处置库预选场地,对该地区的区域地下水流动系统的认识和准确刻画是一项重要的研究内容,它将为处置库的最终确定提供参考依据‘12,13]。4.1地下水流动模型
由于研究区面积很大,水文地质边界较难确定,因此首先使用SRTMDEM数据进行水文空间分析,及数字河网的提取及流域划分,以此从流域整体对边界进行把握。从分析结果(图2)可以看出,北山地区大部分位于黑河流域的西北角,南部一部分地区属于疏勒河流域。
图2水文空间分析结果
Fig.2
Resultsofhydrologicspatialanalysis
地下水模拟区域选取相对完整的北山子流域。
区域东北部边界为北山地区对黑河流域的补给边
界,其他边界则是流域间的分水岭。由于空间分析获得的是地表水分水岭,补给及地形的不对称会引
起地下水与地表水的分水岭不一致,但是对于大区域地下水的模拟来说影响是可以接受的。
北山地区地下水主要有松散岩类孔隙水、碎屑
岩类孔隙一裂隙水和基岩裂隙水3种主要类型¨4I。
研究中使用MODFLOW建立地下水流动模型,并假设其能够代表区域地下水流动。在裂隙发育强烈的区域,模型中将其渗透系数进行显著的放大或缩小来满足要求。模型在平面上离散为180行和275
列,垂向上剖分为三层,模型有效计算单元共计75891个。研究中使用SRTM3DEM数据来构建模型,精细的地表高程刻画有助于提高潜水蒸发的计
算精度㈣。
18(1)董艳辉等:应用并行PEST算法优化地下水模型参数
从研究区内的岩性来看,大面积分布的花岗岩由于存在裂隙,透水性要比变质泥岩好,砂砾岩、碎屑岩透水性更好。为了从一个较简单的模型开始校正,模型中只考虑占主导地位的水文地质特征信息。对于渗透系数,根据不同的岩性分区将其简化为高(K1)、中(1(2)、低(1(3)、极低(K4)4个渗透系数分区。渗透系数参数分区是不连续的,每个分区包含模型中的若干单元(图3)。
图3模型渗透系数分区
Fig.3
Hydrauliceonducfivity
zone
ofmodel
Kl—K4渗透系数
北山地区高差较大,降水量表现出随地势升高
而递增的规律性,因此考虑使用地表高程信息对降
水补给量的空间分布进行分析,以确定不同的降水补给区。对于降雨入渗的补给量划分为高(RCHl)、中(RCH2)、低(RCH3)3个分区。
使用较少的参数分区进行模型校正,可以对整个地表水文地质信息进行清晰的评价。同时,使用
综合敏感陛来研究渗透参数和补给量,进行深一步
的分析。
4.2参数并行优化及分析
待优化参数包括渗透系数及补给量共7个,其中渗透系数在优化中使用其对数值,而补给量为真实值(表1)。选取研究区的17口水并的水位数据
作为观测值。
整个参数优化需要调用176次MODFLOW模
型的进行迭代运算。研究中分别使用原始串行PEST和并行PEST程序进行参数优化。并行计算
测试在中国科学院地质与地球物理研究所地下水数值模拟高性能实验室进行,测试机器拥有8个处理器以及16G内存。使用串行PEST以及使用8个处
万方数据
143
理器并行PEST进行(表2)。
表1
参数定义表
Table1
Parameterdefinitions
表2优化计算耗时表
Table2
ExecutiontimeofsequentialPEsrrandparallelPEST
使用原始的PEST程序进行参数优化需要近6h的时间,由于模型需要多次调试和多个情景分析,因
此计算耗时对研究影响很大。当使用8个处理器并
行计算时,同样的问题只需约1.5h就可完成,效率提高了3.7倍。同时,由于使用OpenMP进行并行化并不需要对算法核心进行修改,因此并行计算和串行程序计算的结果是一致的。
为对参数优化结果进行评价,对每一个参数进行综合敏感度的计算(表3,图4)。由敏感性分析结果可以看出,参数Kl和RCH3的综合敏感度非常低,参数优化过程不能对其进行优化。综合敏感度低是因为在模型的分区中没有观测井信息,即参数Kl和RCH3基本不会对目标函数产生影响,因而不能进行参数估计。在下一步的工作中需尽量收集更多的观测信息来满足优化的需要。其他参数的
综合敏感度较大,说明其优化结果是比较可信的。
表3
参数优化结果
Table3
Resultsofparameterestimation
参数
优化值
K1/m.d一10.31357K2/m.d一15.18E—03K3/m d’‘8.24E一04K4/m.d一15.03E—05RCHl/mm a一13.73E一02RCI}12/mm.8’1
2.92E一01
RCH3mm a一1
1.00E—05
144
50.0
u,
魁40 O
5.’058
_一
罄3m
o
18.787
20.0
11.35I
10.338
r]
10.O。
广1
0.0
0.000I
r]lI
0.013
KI
K2
K3
K4
RCHl
RCH2
RCH3
参数
图4优化参数的综合敏感度
Fig.4
Compositesensitivitiesforestimatedparameters
5结论
本文研究针对大尺度、高精度和高复杂性的大规模地下水模型参数优化数据量大、计算时间冗长
的问题,将PEST参数优化算法进行了并行化,使该方法可以在并行计算机上进行参数优化的并行计算。通过甘肃北山实例应用可以看到,并行化的PEST算法不仅可以通过水位限制条件进行地下水模型参数优化,得到相对最优的参数组合使得模型
最大程度的反映真实的地下水系统。同时,由于使
用了高性能计算机进行并行计算,大大减少了优化时间,提高了效率。
PEST参数优化算法的并行计算是在共享存储并行计算机实现的,在多核电脑上也可以很好地应用,在多核电脑飞速发展的今天,这种并行环境是非常易于实现的。并行PEST参数优化算法为大规模地下水模型的参数优化提供了有力的工具,有着广阔的应用前景。
参考文献
[1]
姚磊华.用改进的遗传算法和高斯牛顿法联合反演三维地下
水流模型参数[J].计算物理,2005,22(4):31l一318.
Yao
Leihun.Parameteridentification
ina
3一Dgroundwaterflow
numericalmodel:animprovedgenetic
algorithmandtheGauss—
Newton
method.ChineseJournalofComputationalPhysics,2005,
22(4):311~318.[2]
刘立才,陈鸿汉,张达政.梯度法在水文地质参数估值中的应用[J].水文地质工程地质,2003,(3):39—41.
IjuUcai,Chen
Honghan,Zh∞gDazheng.Apphcationofgmdi
ent
method
to
calculationofhydrogeologicalparameters.Hydrnge-
ologyand
Engineering
Geology,2003,(3):39-41.
[3]李守巨,上官子昌,刘迎曦等.地下水渗流模型参数识别的模
拟退火算法[J].岩石力学与工程学报,2005,24(s1):5031—
5036.Li
Shouju,ShangguanZichang,LiuYingxi,eta1..Parameter
i-
dentifieationprocedureforgroundwater
flow
modelwithsimulated
万方数据
JournalofEngineeringGeology工程地质学报2010
annealing.ChineseJournalofRockMechanicsandEngineering,
2005,24(s1):5031—5036.
[4]江思珉,朱国荣,胡西嘉.单纯形一模拟退火混合算法反求水文
地质参数及其并行求解[J].地质评论,2007,53(1):92—9r7.
JiangSimin,ZhuGuorong,Hu
Xijia.Hybndsimplex—simulated
annealingmethodforsolfinghydrogeologic
parameters
anditspar-
allel
implementation.GeologicalRenew,2007,53(1):92~97.
[5]刘增荣,崔伟华,王鑫.一种土的非线性弹性本构模型参数的
反演方法[J].工程地质学报,2008,16(3):338—341.
Liu
Zengrong,CuiWei,hun,WangXin.YaoLeihua.Back—
analysys
methodfor
parameters
ofsoilnonlinearelasticconstitutive
models.Journalof
EngineeringGeology,2008,16(3):338—341.
[6]Doherty,J..PEST:Model~Independent
ParameterEstimation
User
Manual(5thed.)[R],2002.
[7]Zyvoloski
G.,KwieklisE.,Eddebbarh
A.A.,et
a1.Thesite—
scalesaturated
zone
flowmodelforYuccaMountain:calibrationof
differentconceptual
modelsandtheir
impact
on
flow
paths[J].
JournalofContaminant
Hydrology,2003,62(63):731—750.
[8]Goegebeu
M.,PauwelsVRN.Improvementof
thePESTparameter
estimation
algorithmthrough
Extended
Kalman
Filtering[J].
Journalof
Hydrology,2007,(337):436—451.
[9]OpenMPArchitectureRenewBoard.OpenMPApplication
Pro-
gram
Interface(Version2.5)[R],2005.
[10]HarhaughA.W.,Banta
E.R.,HillM.C.,eta1..MODF-
LOW-2000,the
U.S.geological
survey
modularground—water
model—-userguide
to
modularization
concepts
andtheground—
waterflow
process[R].USGSOpen—FileReport
00—92,
2000.
[11]HillM.C..Preconditionedconjugate一伊adient2(PCG2),a
computerpmgram
for
solvingground—waterflowequations[R].
USGSWater—ResourcesInvestigationsReport90—4048,2003.
[12]郭永海,刘淑芬,王驹,等.高放废物处置库选址中的水文地
质特性评价[J].世界核地质科学,2007,24(4):233—237.
Guo
Yonghai,LiuShufen,WangJu,eta1..Hydmlogicalper-
fonnanceassessment
on
siting
the
high
levelradioactivewastere—
pository.WorldNuclearGeoscience,2007,24(4):233—237.
[13]郭永海,刘淑芬,吕川河.高放废物地质处置库选址中的水文
地质调查[J].铀矿地质,2005,21(5):296~299.
GuoYonghai,LiuShufen,LuZuanhe,et
a1..Hydrological
in-
vestigation
for
sitting
disposalrepositoryforhighlevelradioactive
waste.Uranium
Geology,2005,21(5):296—299.
[14]刘淑芬,郭永海,王驹等.高放废物地质处置库北山预选区地
下水的形成和分布[J].铀矿地质,2007,23(6):356~362.
LiuShufen,GuoYonghai,WangJu,eta1..Groundwaterforma-
tionanddistributionin
thepreselectedBeishan
al'ea
of
repository
for
Ilighlevelradioactivewaste.UraniumGeology,2005,23(5):
296—299.
[15]董艳辉,李国敏,郭永海.SRTM3DEM在区域地下水数值模
拟中的应用[J].工程勘察,2008,(11):41—43.
Dong
Yanhui,LiGuomin,GuoYonghai.ApplicationofSRTM3
DEM
in
numericalmodekofregional
groundwaterflow.Geoteeh
nical
Investigation&Surveying,2008,(11):41-43.
正在阅读:
应用并行PEST算法优化地下水模型参数06-01
二层框架结构施工组织设计完整版本06-07
叠眉妆的美文11-03
文山县社区文化活动中心规划建设方案05-31
2014年攀枝花市政府工作报告09-12
小鸟你为什么哭泣作文600字07-14
二年级作文猫鼠会谈400字06-20
温暖作文500字07-16
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 地下水
- 并行
- 算法
- 模型
- 优化
- 参数
- 应用
- PEST
- 毕业生大学生活自我评价例文
- 现代物流概论--第1章 概述
- 建筑节能概论 选修论文
- 风险管理【外文翻译】
- 中小学教育培训机构的组织架构、岗位职能要求
- 无心8.0+实验室MOD排序
- 见习报告_范文.doc1
- 初中物理实验操作评分标准及实验记录单
- 托盘使用管理办法0505
- LED照明产品销售渠道建设的7大策略
- 3D打印增材数字化制造技术的发展现状和趋势
- Unit 5 I&39;m watching TV 第一课时教学设计
- 汽车保险发展历程与现状浅析
- 东风汽车集团股份有限公司2009年度第一期短期融资券跟踪评级报告
- 2012广东高考语文试题与答案详解
- 最新2011年职称英语考试大纲词汇表
- 大学生党员挂职街道社区的思考
- 浙大远程计算机网络基础离线作业
- 辽宁瓦房店金刚石矿集区区域成矿控制条件及资源潜力预测
- 2014二级黑龙江省建造师机电工程考试答题技巧