Numerical Modeling and Spatial Moment Analysis of Solute Mobility and Spreading in a Coupled Fracture-Skin-Matrix System

This study deals with transport of solutes through a saturated sub-surface rock formation with well-defined horizontal parallel fractures. For this purpose, a simplified conceptual model consisting of a single fracture and its associated rock-matrix is considered in the presence of a fracture-skin in order to study the mobility and mixing of solutes along the fracture. In this paper, a coupled fracture-skin-matrix system is modeled numerically using finite difference method in a pseudo two-dimensional domain with a constant continuous source at fracture inlet. Flow and transport processes are considered parallel to the fracture axis, while the transport processes in fracture-skin as well as in rock-matrix are considered perpendicular to the fracture axis. Having obtained the concentration distribution along the fracture, method of spatial moments is employed to study the mobility and spreading of solutes. Sensitivity analyses have been done to understand the effect of various fracture-skin parameters like porosity, thickness, and diffusion coefficient. Further, the influence of non-linear sorption and radioactive decaying of solutes are carried out for different sorption intensities and decay constants. Results suggest that the presence of fracture-skin significantly influences the mobility and spreading of solutes along the fracture in comparison with a coupled fracture-matrix system without fracture-skin.     

深部含水层储热系统的数值模拟研究

深部含水层储热系统, 是一种以深部含水层介质为载体的“地热+”多能互补储/供能系统。该系统可将各种形式的能量储存于地下并按需求取出加以利用, 能够弥补能源供需的时空分布的不平衡, 是综合利用多种可再生能源, 实现节能减排的有效途径。深部含水层储热是一种水热型地热资源开发利用的新方式。传统的水热型地热资源开采, 受限于地热尾水回灌导致的热突破, 热田寿命有限。而利用深部含水层储能, 不但可以增强热储能力, 而且能够延长热田寿命。本研究通过数值模拟研究表明, 与传统的水热型地热系统相比, 通过采取深部含水层储热技术, 在相同开采流量下, 可以提升单井供热能力20%。其储热效率可以达到60%~90%, 能够保证更高的地热能提取率; 同时, 与传统水热型地热井寿命受限于热突破时间的情况不同, 深部含水层储热系统可以延长地热井寿命, 实现可持续开采。通过参数敏感性分析的方法进一步评估了不同参数对深部含水层储热性能的影响。通过参数研究发现, 储层渗透率和储热温度等是影响储热能力和效率的关键参数。本研究为后续开展深部含水层储热系统工程的设计与优化提供科学参考。     

非饱和土热湿迁移的数值模拟

本文基于菲利普-迪弗瑞斯(Philip-De Vries)热湿耦合迁移模型,采用一维积分有限差分程序HM1,模拟了一维轴对称热湿耦合迁移问题。模拟中,有意识地调整水分热扩散率D_T,水分等温扩散率D_θ,热导率λ,得到了有关物性参数变化对土壤热湿迁移影响的规律;同时,又有意识地调整某些初始条件(初始含水率)及边界条件(输入热功率),取得了初始含水率及输入热功率变化对土壤热湿迁移影响的规律。这些规律性成果对埋地电缆敷设工程设计具有指导意义。     

Modeling of a Space Heating System Coupled with Underground Energy Storage

Mathematical Geosciences - Heat pump systems and radiant floor heating systems are extensively employed to adjust indoor temperatures. Both types of system can reduce energy consumption and...     

季节冻土区融雪冻土水热耦合模型研究

为描述季节性冻土地区的积融雪及冻土水热过程,采用能量平衡及水量平衡法建立融雪模型,采用有限体积法离散热传导方程和非饱和土壤水运动方程,并对其交界面处进行耦合处理,从而建立季节冻土区融雪冻土耦合数学模型。结合2005年曲玛莱实测站点的观测资料,从不同角度描述积雪变化,冻土活动层内土壤水热特征,检验模型方法的有效性,为揭示该地区的水文运动规律提供有力的支持。     

大陆俯冲带热结构的数值模拟

俯冲带热结构是控制俯冲板块演化的最主要因素之一.前人通过建立解析模型和数值模型对大洋俯冲带热结构进行了一系列研究,发现俯冲板块年龄和俯冲速度是影响俯冲带热结构的关键因素.为了认识大陆俯冲带热结构,特别是理解数值模型结果与岩石学结果之间的差异,我们建立了二维大陆俯冲带运动学和动力学数值模型研究其热结构演化.模型结果显示,如果大陆俯冲板块的俯冲速度与角度和大洋板块一致的话,较低的大陆俯冲带初始温度导致其板块温度比大洋俯冲带低.但是,当大陆俯冲板块的初始温度较高,俯冲速度超慢并且考虑大陆地壳中的放射性元素生热时,模型得到的大陆俯冲带热结构能够解释通过高压和超高压变质岩得到的较热的俯冲板块温度.另一方面,如果俯冲板块与上覆板块存在动力学解耦作用,也能够得到较热的俯冲温压数据.      

Design and Analysis for Modeling and Predicting Spatial Contamination

Sampling and prediction strategies relevant at the planning stage of the cleanup of environmental hazards are discussed. Sampling designs and models are compared using an extensive set of data on dioxin contamination at Piazza Road, Missouri. To meet the assumptions of the statistical model, such data are often transformed by taking logarithms. Predicted values may be required on the untransformed scale, however, and several predictors are also compared. Fairly small designs turn out to be sufficient for model fitting and for predicting. For fitting, taking replicates ensures a positive measurement error variance and smooths the predictor. This is strongly advised for standard predictors. Alternatively, we propose a predictor linear in the untransformed data, with coefficients derived from a model fitted to the logarithms of the data. It performs well on the Piazza Road data, even with no replication.     

人地耦合系统脆弱性研究进展

当前在以人地耦合系统为中心的气候变化、生态风险和可持续管理等研究中,脆弱性理论及其评价已成为地理学及很多相关学科的研究热点,其涵义不断泛化,内容日渐扩展。不同学科背景的学者对脆弱性解释存在很大分歧,这阻碍了研究的科学化,也影响研究成果的实践意义,迫切需要一个可行的框架来系统整合其理论和实践的研究。从审视脆弱性的理论沿革出发,分析其在国内外可持续发展和全球变化领域的应用现状,总结出分歧的主要原因是对相关概念关系的不理解以及研究时空尺度和学科视角的不明确。基于对脆弱性评价的整合框架的研究现状,选取VSD模型作为统领具体的脆弱性评价的依据。VSD模型的优势在于:①有明确的定义,将脆弱性分解为暴露程度、敏感性和适应潜力3个维度;②用方面层—指标层—参数层逐级递进、细化的方式组织评价数据;③有规范评价流程的8个步骤。最后以我国典型的干旱区为例,运用VSD模型构建了示例区脆弱性评价的指标和参数。     

基于多模式降水集成的陆气耦合洪水预报

针对基于单模式陆气耦合模型在洪水预报中存在的不确定性问题,利用三种数值天气预报模式MC2、GEM与T213开展基于多模式降水集成的陆气耦合洪水预报研究。通过淮河流域2007～2009年3场暴雨洪水验证表明:基于多模式降水集成的陆气耦合洪水预报可以有效地减小不确定性、提高精度和稳定性;但基于多模式降水集成的陆气耦合洪水预报对小尺度局地强降水仍存在较大的不确定性,今后需对此展开进一步的研究,以给出更加稳定、可靠的洪水预报结果。     

Coupled Numerical Evaluations of the Geomechanical Interactions Between a Hydraulic Fracture Stimulation and a Natural Fracture System in Shale Formations

Due to the low permeability of many shale reservoirs, multi-stage hydraulic fracturing in horizontal wells is used to increase the productive, stimulated reservoir volume. However, each created hydraulic fracture alters the stress field around it, and subsequent fractures are affected by the stress field from previous fractures. The results of a numerical evaluation of the effect of stress field changes (stress shadowing), as a function of natural fracture and geomechanical properties, are presented, including a detailed evaluation of natural fracture shear failure (and, by analogy, the generated microseismicity) due to a created hydraulic fracture. The numerical simulations were performed using continuum and discrete element modeling approaches in both mechanical-only and fully coupled, hydro-mechanical modes. The results show the critical impacts that the stress field changes from a created hydraulic fracture have on the shear of the natural fracture system, which in-turn, significantly affects the success of the hydraulic fracture stimulation. Furthermore, the results provide important insight into: the role of completion design (stage spacing) and operational parameters (rate, viscosity, etc.) on the possibility of enhancing the stimulation of the natural fracture network (‘complexity’); the mechanisms that generate the microseismicity that occurs during a hydraulic fracture stimulation; and the interpretation of the generated microseismicity in relation to the volume of stimulated reservoir formation.     

油藏数值模拟中地质模型的建模流程与方法

围绕油藏数值模拟过程中三维地质模型的建模技术进行了分析与探讨,详细介绍了基于角点网格模型的建模方法,给出了相应的实现步骤。其主要流程是:首先根据断层数据构造断层模型,在断层模型的基础上构建骨架模型;然后在骨架模型约束下采用地层恢复技术实现含断层的地层模型;最后基于结构模型插值物性参数完成属性模型。以塔河油田缝洞型油藏为例,对建模流程和技术的可行性进行验证,结果表明,其建模结果与专业地质建模软件Petrel相符。     

Progress and Challenges in Coupled Hydrodynamic-Ecological Estuarine Modeling

Numerical modeling has emerged over the last several decades as a widely accepted tool for investigations in environmental sciences. In estuarine research, hydrodynamic and ecological models have moved along parallel tracks with regard to complexity, refinement, computational power, and incorporation of uncertainty. Coupled hydrodynamic-ecological models have been used to assess ecosystem processes and interactions, simulate future scenarios, and evaluate remedial actions in response to eutrophication, habitat loss, and freshwater diversion. The need to couple hydrodynamic and ecological models to address research and management questions is clear because dynamic feedbacks between biotic and physical processes are critical interactions within ecosystems. In this review, we present historical and modern perspectives on estuarine hydrodynamic and ecological modeling, consider model limitations, and address aspects of model linkage, skill assessment, and complexity. We discuss the balance between spatial and temporal resolution and present examples using different spatiotemporal scales. Finally, we recommend future lines of inquiry, approaches to balance complexity and uncertainty, and model transparency and utility. It is idealistic to think we can pursue a “theory of everything” for estuarine models, but recent advances suggest that models for both scientific investigations and management applications will continue to improve in terms of realism, precision, and accuracy.     

三峡库区陆气耦合研究及应用

以延长洪水预见期、提高预报精度为目标,研究气象水文耦合机制,利用数值天气预报模式WRF(Weather Research and Forecasting)驱动分布式VIC(Variable Infiltration Capacity)水文模型,构建三峡库区陆气耦合洪水预报系统,并对2007~2008年期间四场暴雨洪水进行日滚动预报试验。结果表明,WRF模式在三峡库区内有着良好的短期降水预报精度,基于数值天气预报模式和分布式水文模型的陆气耦合洪水预报系统能有效延长三峡入库洪水预见期、提高洪水预报精度,具有较大的应用潜力。     

矩、概率权重矩与线性矩的关系分析

本文叙述了矩、概率权重矩与线性矩之间的关系，并以Г分布为例进行分析。线性矩的计算式为概率权重矩的线性组合，两者的计算结果完全相同。概率权重矩和线性矩均与指定的频率分布型式和作为权重的概率有关，结果的敏感性较差。详细分析了计算参数的近似公式及其精度。对水文应用而言，它们的计算结果仅是估计的初值，需经过合理性分析才能取用。     

Numerical Modeling of the Coupled Mechanical and Hydrological Processes during Deformation and Mineralization in the Mount Isa Block, Australia

Mount Isa is a major Australian and world Pb‐Zn‐Ag mineral province. The wide varieties of mineralization in the province are believed to be closely related to the geodynamic processes of Isan Orogeny, which occurred between ca 1500 and 1620 Ma. In order to understand the geodynamic processes associated with the Isan Orogeny and the giant mineralization systems in the Mount Isa district, a series of numerical models has been constructed to simulate coupled mechanical–hydrological processes, using Fast Lagrangian Analysis of Continua (FLAC), a finite difference computer code. The numerical modeling results have demonstrated that the most probable far‐field stress orientation during the Isan Orogeny is the asymmetrical E–W shortening, which led to greater easternward tectonic movement at the west boundary of the district in comparison with westward movement at the east boundary. During the initial and early stage of the Isan Orogeny, the mechanical and hydrological conditions in the Leichardt Fault Trough of the West Fold Belt are much more favorable for fluid accumulation and mineralization than in the East Fold Belt. The Mount Isan fault zone developed as a high dilation shear zone where the fluids were focused. As the asymmetrical shortening progressed, shortening deformation and shear strain localization became intensified in the eastern part of the orogenic district. The eastern region therefore became a more favorable locality for hydrothermal mineralization. This structural development feature seems to explain why mineralization in the East Fold Belt is generally later than in the West Fold Belt. Fluid production from the Williams–Naraku granites could result in fluid over‐pressuring, and this probably contributed to the extensive brecciation and related mineralization in the East Fold Belt.     

Coupled Reactive Transport Modeling of Redox Processes in a Nitrate-Polluted Sandy Aquifer

The reactive transport modeling of a complicated suite of reactions apparent in the aquifer during the application of N-containing fertilizers is reported. The unconfined sandy aquifer can be subdivided into an oxic zone which contains groundwater with oxygen and nitrate and an anoxic zone characterized by elevated iron and sulfate concentrations in groundwater. Oxygen and nitrate are being reduced by pyrite and organic matter that commonly apparent in the aquifer. The oxidation of pyrite is modeled using the local equilibrium approach, whereas decomposition of organic matter, with the adoption of kinetic approach. The system is buffered by dissolution of aluminum and iron oxides. The modeling process is a two-step procedure. First, the processes are modeled in the one-dimensional (1D) column using PHREEQC code. Subsequently, the calibrated and verified data were copied and used in two-dimensional (2D) PHAST model. Prior to the performance of reactive transport modeling operations with PHAST, a reliable flow model was executed. Finally, predictions are made for the distribution of water chemistry for the year 2008. Model predicts that sulfate derived from the ongoing pyrite oxidation is reduced by the dissolved organic carbon at the higher depth and forms pyrite by the reaction with iron. The results of this study highlight the importance of understanding the interplay between the transport and chemical reactions that occur during the input of nitrate to the aquifer. Reactive transport modeling incorporating the use of a newly developed code PHAST have proved to be a powerful tool for analyzing and quantifying such interactions.     

季节冻土区风积沙土路基冻结过程中水热迁移数值分析

基于Harlan模型和Darcy定律,并考虑温度梯度对水分迁移影响、温度和含水量对水热参数影响以及各种环境气候因素的影响,建立了完全依赖气象资料和水热参数的风积沙土路基冻结过程中水热耦合迁移数学模型,采用全隐式有限差分格式和TDMA迭代法对内蒙古锡林浩特地区沙漠公路207国道K135+000处冻结期间路基水热迁移规律进行了数值模拟.结果表明:该地区道路冻结深度随时间近似线性变化,冻结速度达到2~3 cm·d^-1,最大冻深为3 m左右,冻融时间约为180 d;水分迁移主要发生在冻结锋面附近,从未冻区向冻结区迁移,且随着冻结锋面前移,迁移量逐渐增大;整个冻融期间最大冻深底部层位含水量变化较大,路面下0~50 cm范围内温度变化比较剧烈.