解决方案对交通建模过程废水基于算子分裂技术充电gydF4y2Ba

文摘gydF4y2Ba

地下水运动的耦合和被动运输期间地下水补给与废水导致一个复杂的数学模型,涉及的术语来描述对流扩散、吸附/解吸和/或生物降解,等等。已经发现很难解决这样一个耦合模型分析或数值。本研究采用算子分裂技术,耦合模型分解为两个子不同的内在特征。通过应用逆风有限差分格式的有限体积积分通量对流项,一个隐式解决方案过程推导出求解对流占优方程。色散项是一个标准的中心差分离散方案而dispersion-dominant方程是解决预先处理使用雅可比共轭梯度(PJCG)方法或基于local-one-dimensional托马斯方法方案。方法在这项研究中提出的解决方案应用于地下水补给的示范项目二次废水在高碑店污水处理厂(STP)成功。gydF4y2Ba

1。介绍gydF4y2Ba

通常,地下水系统可以self-maintain液压稳态和化学平衡之前一些外部的刺激,如降雨、刺激、或污水补给,介绍了。然而,地下水系统将经历这样的刺激后很长一段时间不稳定过程,达到新的平衡状态。不稳定的持续时间由外部刺激引起的,特别是污水充电、地下系统经历了复杂的物理/化学/生物过程溶质:污染物对流扩散、吸附/解吸、溶解/沉淀,geo-chemical反应,和/或生物降解。论文的每一个过程可以通过合适的数学模型描述。上述过程的耦合称为活性污染物运输。数学模型描述上述部分或全部过程的耦合现象是极为复杂的。此外,获得解决方案的模型已经被证明是一个艰难的任务。近年来许多调查被执行调查现象hydrophysical过程有或没有反应在各种运输多孔介质([gydF4y2Ba1gydF4y2Ba- - - - - -gydF4y2Ba16gydF4y2Ba])。因为几个过程参与废水的地下水补给,描述耦合系统的数学模型显示强烈的非线性,只能通过数值方法解决。gydF4y2Ba

数值解的方法反应运输方程的套件可以分成两类:一步法方法和两步方法。同时,标准的伽辽金有限元法(FEM)或其修改版本采用流动和溶质运移方程的离散化,大多数现有的数值模型。一步的方法,整个方程系统描述地下水流动的耦合模型和无功传输同时解决,而在两步方法/传输流和反应部分分别对待。一步的方法已经被证明是低效的,很少被用于解决实际问题。大多数的调查集中在两步或多步方法。其中一些采用两步方法基于顺序迭代方法(SIA)来解决无功传输模型。d·谢弗et al。gydF4y2Ba5gydF4y2Ba)建立了一个3 d模型传输、生物化学和化学(TBC)处理无功传输在饱和地下水流动。顺序的解决方案过程使用的advection-dispersion运输方程解决独立地为每一个物种或水组件通过添加源/汇项代表交流阶段从之前的时间步长和化学/生物化学反应。化学/生物化学反应方程解决与运输的浓度变化明确的源/汇项。Tebes-Stevens et al。gydF4y2Ba8gydF4y2Ba)开发了一种无功传输模型FEREACT检查地下水二维稳态流的耦合影响,平衡水物种形成反应,活动控制界面反应。新加坡航空的一个改进的两步解算法是用来把地球化学和微生物反应过程的影响在地下溶质运移的控制方程。gydF4y2Ba

在许多其他现有的调查,采用两步方法基于算子分裂技术解决无功传输模型。在这些研究之前,算子分裂技术应用于解决对流扩散方程,提出了其优势。Zaidel和利未(gydF4y2Ba17gydF4y2Ba)提出了一个基于该方法准确二阶显式有限差分格式获得对流方程的数值解。Russo et al。gydF4y2Ba18gydF4y2Ba)使用一个标准的中心差分方案解决色散方程。近年来,研究活动一直在集中无功传输模型的数值解的方法利用算子分裂技术。巴里et al。gydF4y2Ba19gydF4y2Ba)开发了一个替代算子分裂方法求解化学反应地下水运输模型。程和叶gydF4y2Ba4gydF4y2Ba)提出了一个三维数值模型DHYDROGEOCHEM地下流动,传热和反应化学运输一套完整的地球化学反应是考虑。强耦合弱耦合时用于稳态模拟用于瞬态模拟来节省计算时间。流动方程的离散采用伽辽金有限元法虽然反应的化学运输和传热方程解决混合场或传统的欧拉有限元法。皮卡德的方法被用来处理非线性反应的化学运输和不饱和地下流,而牛顿迭代法来计算化学平衡。巴里et al。gydF4y2Ba11gydF4y2Ba)解决耦合算子分裂方法相比,化学非平衡反应/地下水传输模型。Cirpka et al。gydF4y2Ba9gydF4y2Ba,gydF4y2Ba10gydF4y2Ba和叶等。gydF4y2Ba16gydF4y2Ba)使用新加坡航空两步方法和算子分裂两步方法的比较。gydF4y2Ba

虽然这样的被动运输过程的许多方面进行了调查实验和数值,现有的研究大多是跟踪到一个特定的问题和我们理解的耦合过程远非令人满意。实际问题引起的工程要求,地下水补给与废水及其可能的应用领域规模激励我们洞察机制/基本反应动力学过程和它们之间的交互进行进一步调查被动运输的现象,并建立数值模型有效地解决耦合问题。在目前的研究中,一个新的解决方案开发过程基于算子分裂技术解决耦合的对流扩散模型和被动运输。部分gydF4y2Ba2gydF4y2Ba州数学模型描述饱和地下水流动和污染物传输的细节。处理不同的水位,基于体积积分的有限体积方法介绍了偏微分方程(pde)。部分gydF4y2Ba3gydF4y2Ba描述算子分裂技术的应用无功传输模型。建立一个隐式逆风有限差分格式求解对流占优方程。数值计算解决方案以显式的方式通过优势的潜在的流场。扩散通量项是一个标准的中心差分离散方案。托马斯PJCG方法或方法基于局部一维方式适合解决dispersion-dominant方程。应用数值计算方法在本研究开发的地下水补给的示范项目在高碑店二级污水STP节中描述gydF4y2Ba4gydF4y2Ba。这个示范项目的细节被称为(gydF4y2Ba20.gydF4y2Ba]。gydF4y2Ba

2。数学模型gydF4y2Ba

地下水补给的物理模型如图gydF4y2Ba1gydF4y2Ba。在多孔介质地下水流动的基本假设是(a)达西定律,impl,流速与压力梯度成正比;(b) nondeformable多孔介质骨架;(c) Bossinesq地下水的密度近似;(d)等温和局部平衡的化学反应,生物降解;(e)固定稳固阶段,Fickian分散混合和理想。gydF4y2Ba

2.1。的偏微分Equations-PDEsgydF4y2Ba

2.1.1。控制方程的饱和地下水流动gydF4y2Ba

饱和地下水流模型的控制方程可以表示为gydF4y2Ba 在哪里gydF4y2Ba是特定存储系数,gydF4y2Ba和gydF4y2Ba地下水的压缩系数和压缩系数的骨架,分别;在笛卡儿坐标系统gydF4y2Ba是水头;gydF4y2BakgydF4y2Ba相关的渗透系数张量吗gydF4y2Ba通过gydF4y2Ba,在那里gydF4y2Ba是内在的渗透率张量,gydF4y2BaμgydF4y2Ba地下水的粘度,gydF4y2BaggydF4y2Ba重力加速度,gydF4y2BazgydF4y2Ba基准面的垂直坐标朝上低于地下水位。gydF4y2Ba是协会的源项与充电或生产。gydF4y2Ba

通常情况下,通过选择坐标系统配合的主要渗流方向,可以简化为渗透系数张量gydF4y2Ba

与水的解决方案(gydF4y2Ba1gydF4y2Ba),速度场可以用达西定律计算gydF4y2Ba 在哪里gydF4y2Ba。gydF4y2Ba

2.1.2。控制方程的被动运输gydF4y2Ba

溶质的运输过程污染物同时发生与地下水运动,包括对流、土壤吸附/解吸,和/或生物降解。为组件gydF4y2Ba米gydF4y2Ba可溶性物质,无功传输模型可以表示由以下方程:gydF4y2Ba 在上面的方程中,参数gydF4y2Ba是岩石的孔隙度/土壤;gydF4y2Ba是溶质的浓度有关系吗gydF4y2Ba在注入浓度/生产。gydF4y2Ba代表组件的质量流量gydF4y2Ba米gydF4y2Ba由于分散,可以表示为gydF4y2Ba 在这gydF4y2Ba弥散张量(包括分子扩散),可以计算吗gydF4y2Ba 在哪里gydF4y2Ba是指分子扩散;gydF4y2Ba是曲折的因素;gydF4y2Ba和gydF4y2Ba纵向/切向扩散系数;gydF4y2Ba(gydF4y2Ba)或= 0 (gydF4y2Ba)(gydF4y2Ba)。通常,分子扩散是小到可以被安全地排除在分析([gydF4y2Ba21gydF4y2Ba,gydF4y2Ba22gydF4y2Ba])。gydF4y2Ba

这个词gydF4y2Ba右边的gydF4y2Ba4gydF4y2Ba)代表浓度变化由于土壤吸附,可表达的水浓度gydF4y2Ba 在哪里gydF4y2Ba土壤干容重;gydF4y2Ba是固相吸附的污染物浓度粒子表面;和gydF4y2Ba是线性分布系数gydF4y2Ba。gydF4y2Ba

生物降解的术语可以写成gydF4y2Ba 在哪里gydF4y2Ba降解系数,标1、2、3代表有氧状态,缺氧模式,分别或厌氧模式。gydF4y2Ba被称为undegradable浓度。有机污染物的生物降解在特定时间可以确定为三种降解模式之一。因此,上标gydF4y2Ba我gydF4y2Ba将被忽视的和实际值将分配给退化系数gydF4y2Ba根据生物降解过程中,好氧、缺氧、厌氧。gydF4y2Ba

用(gydF4y2Ba7gydF4y2Ba)和(gydF4y2Ba8gydF4y2Ba)(gydF4y2Ba4gydF4y2Ba),我们有gydF4y2Ba 在哪里gydF4y2Ba被称为发育迟滞系数。非降解性组件,我们有gydF4y2Ba而对于nonadsorbable组件,gydF4y2Ba。gydF4y2Ba新源项是由于充电/生产和生物降解。gydF4y2Ba

2.1.3。初始和边界条件gydF4y2Ba

将遇到的特定问题的充电站点,初始水头梯度可以被忽视而溶质污染物平均初始分布gydF4y2Ba 在哪里gydF4y2Ba。gydF4y2Ba

有限差分方程(gydF4y2Ba29日一gydF4y2Ba),(gydF4y2Ba29个bgydF4y2Ba)和(gydF4y2Ba29度gydF4y2Ba)可以表示为gydF4y2Ba 在哪里gydF4y2Ba 解决方案(gydF4y2Ba30度gydF4y2Ba)的最终解决无功传输模型:PDE (gydF4y2Ba4gydF4y2Ba)或有限体积法的形式(gydF4y2Ba13 bgydF4y2Ba)。证明了本方法的有效性,我们还计算的显式解(gydF4y2Ba13 bgydF4y2Ba)。250的目标字段gydF4y2Ba110米gydF4y2Ba20米(新西兰= 10)被选中来模拟做的分布集中在180天。明确的解决方案一个很小的时间步长为0.1天可以作为近似的解决方案。对比目前的隐式解决方案,明确方案展示一个好的协议。gydF4y2Ba

4所示。应用程序gydF4y2Ba

目前的解决方案过程数值模拟工作的一部分,为示范工程进行了地下水补给的二级污水在高碑店STP (gydF4y2Ba20.gydF4y2Ba]。目标站点在北京冲积平原,那里有一个紧密的粘土层表面16 m以下,可以设置为基准。地下水位是9.7米。三个充电池是为替代操作(图gydF4y2Ba6gydF4y2Ba)。一个生产井和三个观测井用以监测地下水运动和污染物运输。充电实现显示充电率很低(20 - 30吨/天)的充电模式充电池。这样的充电速率,没有什么可以被监控观测井,从仿真结果也不能区分。渗透系数测试值的解释这一现象:0.12 /天(0 - 4米深);0.48 /天(4 - 8米深度);1.92 /天(9 - 12米的深度);3.84 /天(13 - 16米深度)。gydF4y2Ba

这些实践的基础上,我们建议注入井的充电模式充电水直接进入层水力传导率高,给仿真结果。建立了数值模型网格系统gydF4y2Ba模拟域的1000米gydF4y2Ba600米gydF4y2Ba20米。初始条件是地下水位,深度为9.7米;DOC浓度为2.4 mg / L;氧浓度、0.01 mg / L。undegradable浓度被认为是1.2 mg / L。根据实验室检测,吸附分配系数为0.00012(公斤/ l)。有氧环境退化系数(> 0.9 mg / L),缺氧环境(0.9 mg / L > > 0.01 mg / L),和厌氧环境(0.01 mg / L >)是0.139,0.07,和0.02,分别。液压和被动运输过程模拟的一些参数估计,包括hydraulic-steady-state;hydraulic-penetration;maximum-recharge; reused-water-quality. The results are partly presented in this paper.

地下水质量的决定性因素可能重用的二级污水通过地下充电。在这里,我们选择医生作为代表污染物对地下水质量进行评估。图gydF4y2Ba7gydF4y2Ba显示了医生分布在180天的补给率产量200吨/天,200吨/天。gydF4y2Ba

可以看出,DOC浓度升高随着时间的增加,尤其是对一个接近排水好的网站。没有足够的生物降解污染物达到生产井。井300米之外,DOC浓度的增加可以被忽视,即使720天充电。生产井的水质可以提高通过扩大其排水距离网站和/或通过增加浓度提高生物降解。做水的质量浓度的影响在生产井150远离充电站点表所示gydF4y2Ba1gydF4y2Ba。gydF4y2Ba

5。结论和讨论gydF4y2Ba

重用的废水通过地下充电是一种有效的水资源管理模式。流体由于废水总是伴随着活性溶质运输补给污染物。为了避免外来污染地下水的补给,必须进行数值模拟来评估对地下水水质的影响和充电前的地下环境中实现。地下水运动的耦合和被动运输导致复杂的非线性模型,已被证明是很难获得数值解。基于算子分裂技术,本研究发展一个新的解决方案过程来解决无功传输模型的饱和流与废水在浅层地下水补给。处理不同水位、有限体积的积分偏微分方程(pde)的介绍。采用算子分裂技术,无功传输模型分解成一个对流占优方程加上dispersion-dominant方程。一个隐式逆风有限差分格式用于求解对流占优的方程。数值计算解决方案以显式的方式通过优势的潜在的流场。一个标准的中心差分方案应用于dispersion-dominant方程。 Either the PJCG method or the Thomas method based on local one-dimensional is suitable for solving the dispersion-dominant equation. The explicit solution with a very small time step, which can be taken as the approximate of the real solution, is also computed to demonstrate the validity of the present method.

模拟二维对流占优被动运输过程,Cirpka et al。gydF4y2Ba9gydF4y2Ba)开发了一种新的基于cell-centered streamline-oriented网格有限体积方法方案。在他们的方法中,运输方案是一个适应的主要方向,因此,横向通量减少。计算后的水力模型基于矩形网格,包含四边形streamline-oriented网格细胞应该是生成的,而空间梯度沿流线和正交的方向是必需的。由此产生的模型是解决一个算子分裂方法,对流运输的解决非线性显式方法分散运输时隐式地解决。虽然这种方法显示其优势(gydF4y2Ba11个gydF4y2Ba),(gydF4y2Ba11 bgydF4y2Ba),它只能用于二维问题,因为使用流线。对于瞬态问题,每个时间步的一代的面向数据流网格将耗费时间。此外,只有明确的解决方案可以获得对流运输,这始终是被动运输的主要过程。而Cirpka等的方法。gydF4y2Ba9gydF4y2Ba),目前的解决方法并不局限于二维问题。它可以自由地应用于三维瞬变流动的问题。gydF4y2Ba

应该提到,在这项研究中所描述的解决方案过程效率高是被动运输功能由于二级污水补给下溶质污染物的浓度相对较低。浓度变化由于吸附和生物降解可以近似莫诺类型。更常见的情况下的耦合非线性被动运输问题,辛普森和乡下人gydF4y2Ba23gydF4y2Ba]开发了分裂算子方法应用于无功传输和莫诺动力学分析了删除错误的边界和内部错误。gydF4y2Ba

引用gydF4y2Ba

1. j . e .容易和m . Jekel”模拟地下水补给与先进的废水处理,”gydF4y2Ba水科学与技术gydF4y2Ba,33卷,不。年级,409 - 418年,1996页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
2. r . Therrien和e·a·Sudicky”三维分析variably-saturated discretely-fractured多孔介质中流动和溶质运移,”gydF4y2Ba《污染物水文gydF4y2Ba,23卷,不。1 - 2,在美国,1996页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
3. a . Habbar和o . Kolditz反应运输过程的数值模拟在裂隙多孔介质,”gydF4y2Ba学报》国际会议上地下水质量(《GQ》98)gydF4y2BaTruebingen,页21 - 25日,德国,1998年9月。gydF4y2Ba视图:gydF4y2Ba谷歌学术搜索gydF4y2Ba
4. H.-P。程和G.-T。叶”,开发和示范应用地下流动的三维数值模型,传热,和反应性化学运输:3 dhydrogeochem,”gydF4y2Ba《污染物水文gydF4y2Ba,34卷,不。1 - 2,页47 - 83,1998。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
5. d·谢弗w·谢弗和w . Kinzelbach模拟含水层1中有关生物降解反应过程。三维结构被动运输模式”,gydF4y2Ba《污染物水文gydF4y2Ba没有,卷。31日。1 - 2、167 - 186年,1998页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
6. c . i Steefel和p . c . Lichtner”,多组分反应传输在离散fractures-I:控制反应几何面前,“gydF4y2Ba《水文gydF4y2Ba,卷209,不。1 - 4、186 - 199年,1998页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
7. c . i Steefel和p . c . Lichtner”,多组分反应传输在离散fractures-II: Maqarin超碱性地下水的渗透,约旦,自然模拟网站”gydF4y2Ba《水文gydF4y2Ba,卷209,不。1 - 4、200 - 224年,1998页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
8. c . Tebes-Stevens a . j . Valocchi j . m . VanBriesen b . e . Rittmann,“多组分传输与耦合的地球化学和微生物反应:模型描述和实例模拟,”gydF4y2Ba《水文gydF4y2Ba,卷209,不。1 - 4,8-26,1998页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
9. o . a . Cirpka e·o·Frind, r . Helmig”Streamline-oriented交通模型的网格生成二维域包括井,“gydF4y2Ba水资源的进步gydF4y2Ba,22卷,不。7,697 - 710年,1999页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
10. o . a . Cirpka e·o·Frind, r . Helmig”数值方法对矩形和streamline-oriented电网无功传输,”gydF4y2Ba水资源的进步gydF4y2Ba,22卷,不。7,711 - 728年,1999页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
11. d . a . Barry m·克拉普k Bajracharya表示,h .舞会和c·j·坎宁安,“比较split-operator方法求解耦合化学非平衡反应/地下水传输模型,”gydF4y2Ba数学和计算机模拟gydF4y2Ba,53卷,不。1 - 2、113 - 127年,2000页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
12. h .舞会、d·a·巴里和g·b·戴维斯”数值模拟设计和评价地下水的修复方案,“gydF4y2Ba生态模型gydF4y2Ba,卷128,不。2 - 3、181 - 195年,2000页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
13. k·t·b·三联和e·a·Sudicky“多组分模拟wastewater-derived氮和碳的浅非承压含水层。即模型配方和性能。”gydF4y2Ba《污染物水文gydF4y2Ba卷,47号1,53 - 84年,2001页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
14. k·t·b·麦克卡利·e·a·Sudicky和w·d·罗伯逊”数值模拟取消化粪池的细粒度的脱氮层从浅层地下水硝酸盐,”gydF4y2Ba《污染物水文gydF4y2Ba,52卷,不。1 - 4,29-55,2001页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
15. 美国大肠Serrano”下溶质传输非线性吸附和衰变”,gydF4y2Ba水的研究gydF4y2Ba,35卷,不。6,1525 - 1533年,2001页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
16. G.-T。叶,m·d·西格尔,M.-H。李”,数值模拟耦合不定地饱和流体流动和无功传输与快速和缓慢的化学反应,”gydF4y2Ba《污染物水文gydF4y2Ba卷,47号2 - 4、379 - 390年,2001页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
17. j . Zaidel和b·李维”数值方案阶段的解决方案,三个组件流在多孔介质和传输方程,”gydF4y2Ba连续介质力学的数值方法gydF4y2Ba11卷,第89 - 75页,1980年。gydF4y2Ba视图:gydF4y2Ba谷歌学术搜索gydF4y2Ba
18. j·d·Russo Zaidel, a . Laufer”部分饱和非均匀土壤中溶质运移的随机分析,“gydF4y2Ba水资源研究gydF4y2Ba,30卷,不。3、769 - 779年,1994页。gydF4y2Ba视图:gydF4y2Ba谷歌学术搜索gydF4y2Ba
19. d·a·巴里和c·t·米勒,k . Bajracharya表示,“替代split-operator方法求解化学反应/地下水传输模型,”gydF4y2Ba水资源的进步gydF4y2Ba,19卷,不。5,261 - 275年,1996页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
20. 朱,g .赵k .彭,y,“饱和造型污水补给及其对地下水水质的影响,”gydF4y2Ba废水再利用和地下水质量gydF4y2BaIAHS出版285年,页100 - 109,IAHS出版社,瓦林福德,英国,2004年。gydF4y2Ba视图:gydF4y2Ba谷歌学术搜索gydF4y2Ba
21. j .熊,gydF4y2Ba在多孔介质流体动力学gydF4y2Ba爱思唯尔,纽约,纽约,美国,1972年。gydF4y2Ba
22. j·d·Russo Zaidel, a . Laufer”流动数值分析和运输在结合多相渗流zone-groundwater系统中,“gydF4y2Ba水资源的进步gydF4y2Ba,24卷,不。1,49 - 62年,2000页。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba
23. m·j·辛普森和k·a·兰德曼”分裂算子的分析方法应用于无功与莫诺动力学运输,”gydF4y2Ba水资源的进步gydF4y2Ba,30卷,不。9日,第2033 - 2026页,2007年。gydF4y2Ba视图:gydF4y2Ba出版商的网站gydF4y2Ba|gydF4y2Ba谷歌学术搜索gydF4y2Ba

版权gydF4y2Ba

版权©2008 Shutang朱镕基等。这是一个开放分布式下文章gydF4y2Ba知识共享归属许可gydF4y2Ba,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。gydF4y2Ba