非均质含水层参数识别和分数阶对流弥散方程数值解法研究

来源 :南京大学 | 被引量 : 0次 | 上传用户:superficalness
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
我国淡水资源相对匮乏,而地下水资源作为淡水资源的重要组成部分,在经济和生活中占有重要地位。随着我国地下水资源开发利用程度及工业活动的增加,地下水污染问题越来越受到关注。地下水数值模拟作为一种通用手段,在地下水资源评价、地下水污染物迁移预测、地下水修复技术指导方面有广泛应用。自然界的含水层,受形成条件及地质作用影响,多具有强烈的非均质性,导致污染物在其中的迁移,表现出反常特征,如污染物穿透曲线的拖尾现象和早到达现象。如何刻画含水层的非均质性,开发出能精确模拟非均质含水层中污染物迁移转化规律的方法,并应用到实际案例为地下水开发利用及修复提供指导,成为当前地下水溶质运移模拟研究的一个热点。当前对非均质介质中地下水溶质运移模拟,主要分为两种方法:(1)在典型单元体尺度上,对含水层进行细化,精确赋值模型每个网格的参数,然后用基于Fick定律的传统对流弥散方程模拟;(2)修改传统的对流弥散方程,得到新的可用于描述溶质反常运移性质的分数阶对流弥散方程进行模拟。针对方法(1),局域化集合卡尔曼滤波(Ensemble Kalman Filter,EnKF)是一种有效的参数识别方法,其通过同化含水层观测数据,校正模型参数的估计,达到对含水层性质精确刻画。本文修改课题组已有模型,将局域化EnKF同时同化地下水水位和溶质运移数据,识别对数渗透系数场,探究相对于同化单一种类的观测数据,同时同化两种观测数据对对数渗透系数场识别结果的改进,并探讨了溶质运移范围对渗透系数识别的影响;然后尝试同化地下水溶质运移观测数据,探究仅有对流弥散作用下,局域化EnKF在识别对数弥散度场的作用,并探讨了实现数目、初始猜想、观测点数目、观测频率、观测点空间分布和观测误差对弥散度场识别结果的影响;最后,探究局域化EnKF在同化溶质运移数据,同时识别对数弥散度场和对数吸附参数场的效果。得到如下结论:①局域化EnKF同时同化水位和溶质运移数据识别对数渗透系数场时,参数估计效果比同化单一种类数据好,说明应用多种类数据进行参数识别的有效性;若溶质运移范围有限,仅同化只对溶质运移范围内参数识别有效。②利用局域化EnKF同化溶质运移数据识别对数弥散度场,可以有效地对弥散度场进行识别,不同因素对参数识别结果影响如下:模型实现数目有一定要求,过少实现数不能对参数进行有效识别,过多实现数会增大计算负担;观测点数目并非越多越好:过少的观测点不能提供足够的参数识别信息,过多的观测点会由于观测误差的累积导致参数识别结果变差;初始猜想越接近真实值,利用相同观测数据,参数识别结果越好;观测误差越大,为局域化EnKF系统提供的误导信息越多,弥散度场识别效果越差;观测井位置、数目及观测频率对参数识别结果有较大影响,特定的点位和时间观测,能够使参数识别结果最好。③利用局域化EnKF同化含水层溶质运移数据(同时存在对流弥散和吸附作用),可以同时识别对数弥散度场和吸附参数场,展示了利用局域化EnKF同时识别两种参数的有效性。实现数目对同时识别多种参数有影响:某一个实现数目,可能只对单一参数识别结果较好,选择实现数目时,除了考虑计算消耗,参数类别与实现数目对结果的综合作用也应考虑;观测点位的空间分布及观测频率对参数估计结果也有较大影响:整体趋势上,观测数据越多参数识别结果越好;对于文中模型,最优观测间隔是1.0天。对于方法(2),虽然分数阶对流弥散方程理论研究有很大进展,但是其求解精度和计算效率一直限制其在实际算例中的推广,因此高效的求解方法一直是当前研究的一个热点。本文在粒子追踪法求解时间分数阶对流弥散方程的基础上,提出了变权重粒子追踪方法,对每一个粒子赋予可以变化的行走概率和权重,粒子每次跳跃后,根据到达点位概率与粒子原权重的比值,相应调整粒子的数量、权重及行走概率,以求提高计算结果精度。变权重粒子追踪法求解时间分数阶对流弥散方程的一维和二维算例结果表明:相对于传统粒子追踪方法,变权重粒子追踪法可以有效提高计算结果精度;在提高计算结果精度的同时,变权重粒子追踪法不会增加计算量,造成额外的计算负担。用变权重粒子追踪方法求解空间分数阶对流弥散方程,并没有明显提高粒子追踪方法的精确度,效率较低的原因可能与粒子跳跃步长变化较大,两次跳跃之间的位置缺乏关联性。
其他文献
第一章血清miR146a与miR-155可作为强直性脊柱炎早期诊断分子标记物的筛选研究目的:比较强直性脊柱炎患者(ankylosing spondylitis,AS)与正常对照组血清microRNAs(miRNAs)表达谱的差异,并评估miRNAs作为早期诊断AS生物标志物的潜在价值。方法:选择10例AS患者(病例组)和10例健康成人(正常对照组),AS的诊断均符合1984修订的纽约标准。首先,采
限制性三体问题模型下的平动点在空间探测任务中具有非常重要的作用。平动点既是对空间天体观测的理想区域,又可以作为人类进行太阳系行星际空间探测任务的中转站。学者们围绕平动点周期轨道和到达平动点的转移轨道设计进行了大量工作。本文以此为背景,针对空间任务中涉及到的相关动力学问题进行了研究,本文的主要研究内容和创新点如下:构造了圆参考轨道相对运动模型下任意平动点附近的周期构型。首先将平动点(主星)附近的周期
脂肪肝(fatty liver disease,FLD)是各种病因导致的肝脏中脂类沉积过多,肝脏细胞有明显的脂滴存在。非酒精性脂肪肝(non-alcoholic fatty liver disease,NAFLD)是正在迅速蔓延的全球性健康问题,在欧美和澳洲等发达国家的发病率为20~40%,在中国NAFLD是仅次于病毒肝炎的第二大肝病。严重的NAFLD可促进肝纤维化、肝癌和2型糖尿病等,因此,预防
如何高效地调控声波的传播是声学研究的一个重要问题。近年来,声学超构材料(acousticmetamaterials)和声子晶体(phononic crystals)作为新型的声学功能材料,展现出了传统声学材料所不具备的奇异特性,为声波的操纵带来了全新的自由度。声子晶体和声学超构材料之间没有严格的区分,都属于声学人工结构材料。一般来说,声子晶体主要研究的是声波在周期性结构中的传播问题,声子晶体所展现
白云岩成因和白云岩储层演化过程中相关流体的来源和性质对解决相关问题至关重要。从白云岩成因角度来看,涉及到的模式包括咸化海水模式、海水与大气水混合模式及卤水渗透回流模式等;从储层的形成来看,既有同生期的大气降水淋滤,也有成岩期的溶蚀,还有深部热流体的改造等。因此,要把白云岩及其相关储层的形成过程和机理研究清楚,流体来源的确定是关键环节。但从目前的研究现状来看,这一问题还没有得到很好解决。针对这一科学
自从2004年石墨烯被成功制备以来,二维材料优异的物理化学性能震惊了科研工作者,其在各行各业的应用越来越广泛,深刻地改变着人们的生活。对新型二维材料的探索和研究作为前沿科学之一,越发吸引众多国家纷纷从战略层面部署发展计划。具有类似石墨烯结构的二维过渡金属硫化物(Transitional metal dichalcogenides,TMDC)具有优异的光学、电学、量子物理等特性,在各领域取得了诸多研
功能材料是二十一世纪高新技术发展最重要的先导性和基础性材料之一,在国际上,对于新型功能材料的开发和研究一直都很活跃。它的进步推动着诸如信息、电子、机械、化工、航天和通讯等尖端的高科技产业的飞速发展。稀土正铁氧体材料作为一种新型电子材料,目前已在信息产业,电子传感等领域被广泛应用,因此具有较大的研究价值和应用前景,一直是凝聚态物理和材料科学领域探索和开发功能体系的热点之一。力学谱是研究固体材料中的相
借助水文地质学和同位素水文地球化学理论研究油气成藏理论是学科交叉的一个热点。我国天然气资源储量丰富,但常规气开采成本较高,因此需要开展非常规油气资源的勘探。天然气水合物是高效清洁的新能源,近年来吸引各国的关注。海洋陆坡区是天然气水合物的主要储库,它只能在低温高压条件下存在。天然气水合物的一个重要研究方向是甲烷及其气源。甲烷是有机质的分解产物之一,通过热力学计算可以得到甲烷在海洋沉积物孔隙水中的溶解
砷既是人体必须的微量元素,也是危害非常严重的污染元素。不同赋存形式的砷,毒性不同,通常无机砷的毒性要高于有机砷,As(Ⅲ)的毒性要高于As(V)。矿山开发往往是导致地区性砷污染的主要原因,因此揭示矿山环境中次生砷矿物的组成及其可溶性对评估砷的环境危害、治理砷污染具有重要意义。本文首先系统研究了安徽铜陵铜金矿区和湖南石门雄黄(AsS)-雌黄(As2S3)矿山中砷的赋存形式和主要的次生砷矿物类型,从矿
邹平火山岩盆地位于山东省邹平县南部,盆地内广泛分布有早白垩世侵入岩和火山岩,其中王家庄石英二长岩、碑楼二长斑岩和南洞子二长岩与Cu(-Mo-Au)矿化有关。本文在详细的野外地质工作基础上,通过盆地内侵入岩和火山岩锆石U-Pb定年、全岩元素地球化学、全岩Sr-Nd-Pb-Hf同位素、锆石原位Hf-O同位素,王家庄Cu(-Mo)矿床辉钼矿Re-Os定年、矿石硫化物S、Pb同位素、含矿石英H-O同位素和