Numerical modeling study of mineralization induced by methane cold seep at the sea bottom     

《Marine and Petroleum Geology》2016

In order to investigate the response of authigenic minerals to gas hydrate geo-systems, the biogeochemical processes and its induced mineralization were predicted by employing the comprehensive reactive transport modeling approach. Based on the available data extracted from the northern continental slope area of the South China Sea, a 1-D vertical column model was developed. Three cases with different upward methane flux rates and three cases with different mineral compositions, i.e., a total of six cases were designed to investigate the effects of variations in the depth of sulfate methane transition zone (SMTZ) and in the mineral composition on the formation of authigenic minerals. The simulation results indicate that the SMTZ depth influenced by both the upward methane flux rate and the initial composition played an important role in the formation of authigenic minerals. The AOM reaction is intensive at the interface, and the precipitation amount of calcite is large, which is mainly controlled by AOM. When the methane leakage rate is 20 times higher than the base case, aragonite starts to precipitate. During the simulation, oligoclase, k-feldspar, smectite-Na, smectite-Ca, chlorite dissolved. Our study specific to this area as a starting point may provide a quantitative approach for investigating carbonate and pyrite formation in hydrate-bearing sediments accounting for methane oxidation and sulfate reduction. The method presented here and the model built in this study can be used for other sites with similar conditions. In addition, this study may serve as an indication for the potential natural gas hydrate reservoir in depth, and is also significant for marine carbon and sulfur cycle.  相似文献   

Does compaction-induced subsidence control accommodation space at the top of prograding carbonate platforms? Constraints from the numerical modelling of the Triassic Esino Limestone (Southern Alps,Italy)     

《Marine and Petroleum Geology》2016

The demise of the high-relief, steep-slope, prograding Ladinian-Early Carnian carbonate platforms of the Esino Limestone (Central Southern Alps of Italy) is marked by subaerial exposure of the platform top associated with different erosional (mainly karst-related), depositional and diagenetic processes (Calcare Rosso). The exposure-related deposits consist of three major facies associations: 1) residual soils with thin lenses of conglomerates with black pebbles, and, locally, weathered vulcanites; 2) chaotic breccia lenses irregularly distributed in the uppermost part of the Esino Limestone carbonate platform, interpreted as collapse breccias in karstic setting: 3) inter-supratidal carbonate cycles with dissolution and development of paleosols and tepee structures.Facies distribution follows the sub-environments of the underlying Esino Limestone. Facies 1 and 2 typically characterize the core of the platform, covering the underlying inner platform facies. Facies 3 instead develops toward the edge of the platform, above reef-upper slope facies of the prograding facies of the Esino Limestone. The thickness of facies 3 decreases toward the core of the platform. Facies distribution reflects differences in the accommodation space and sedimentary processes from the rim (highest accommodation, favouring the deposition of peritidal-supratidal carbonates) to the core (reduced accommodation, causing pedogenesis and karstification) of the carbonate system.The observed thickness changes may be controlled by different factors: 1) syndepositional tectonics, 2) subsidence induced by magmatic activity or 3) differential subsidence controlled by the stratigraphic architecture of the Esino Limestone platform and adjoining basins. As evidence of tectonics was not observed and the presence of volcanic bodies is only documented tens of km away from the study area, the scenario involving the creation of accommodation space by compaction of the basinal sediments (resedimented, fine-grained calciturbidites) during the progradation of the carbonate platform is here investigated. Numerical modelling was performed to verify the compatibility of compaction-induced subsidence with the observed depositional architecture. The models were built to simulate the architectural evolution of the platform by progressively adding layers from deepest to shallowest, while compacting the underlying sediments, in order to evaluate compaction-induced subsidence (and accommodation space for the Calcare Rosso) after the deposition of the youngest platform strata. Modelling results allow us to conclude that the wedge geometry of the Calcare Rosso, deposited on top of the extinct Esino carbonate platform, can be explained by subsidence controlled by compaction of the basinal sediments present below the early-cemented, fast prograding platform slope deposits.  相似文献   

2D and 3D coexisting modes of thermal convection in fractured hydrothermal systems - Implications for transboundary flow in the Lower Yarmouk Gorge     

《Marine and Petroleum Geology》2016

Numerical investigations of 2D and 3D modes of large-scale convection in faulted aquifers are presented with the aim to infer possible transport mechanisms supporting the formation of thermal springs through fractured aquicludes. The transient finite elements models are based on idealized structural features that can characterize many hydrothermal systems. The sensitivity analysis of the fault permeability showed that faults cross-cutting the main regional flow direction allow groundwater to be driven laterally by convective forces within the fault planes. Therein thermal waters can either discharge along the fault traces or exit the fault through adjacent permeable aquifers. In the latter case, the resulting flow is helicoidally and transient. The location and the spacing between discharge areas can migrate with time, is not strictly constrained to the damage zones and reflects the wavelength of the multicellular regime in the fault zone.An illustrative example based on simplified structural data of the Lower Yarmouk Gorge (LYG) is presented. The numerical calculations indicate that crossing flow paths result from the coexistence of fault convection, developing for example along NE-SW oriented faults within the Gorge, and additional flow fields. The latter are induced either by topography NS gradients, e.g. perpendicular to the major axe of the Gorge, or by local thermal convection in permeable aquifers below the Eocene aquiclude. Sensitivity analysis of fault hydraulic conductivity (K) and the analytical solutions based on viscous-dependent Rayleigh theory show that K values between 2.3e⁻⁷ m/s and 9.3e^{– 7} m/s (i.e. 7 m/yr and 30 m/yr, respectively) favor coexisting transport processes. The uprising thermal plumes spread over several hundred meters forming clusters of springs, in agreement with observation, and which temperature fall within the measured ranges, i.e. 20 °C−60 °C. To some extent the models also reproduced the transient behavior of the spring temperature. Owing to the idealized nature of the presented models, the numerical results and the associated analytical solution can be applied to study the onset of thermal convection and resulting flow patterns of any fractured hydrothermal basin.  相似文献   

FLAC3D在边坡数值分析及位移监测中应用     

张颖  孙兆涛  裴仲文 《探矿工程》2018,45(3):57-61

以辽宁省岫岩县某山头边坡为例,通过槽探、钻探、大型剪切试验等多种手段,找出合理潜在滑动面。利用FLAC3D软件对该边坡进行三维数值模拟,分析应力应变关系,得出变形较大及应力集中区域,同时利用强度折减法得出边坡的安全系数,评价其稳定性。将从边坡后缘张拉裂隙上布设的位移监测点采集到的数据,与数值模拟结果进行对比分析,验证数值模拟的可靠性,为边坡抗滑设计提供依据。  相似文献   

Seismicity and fault aseismic deformation caused by fluid injection in decametric in-situ experiments     

《Comptes Rendus Geoscience》2018,350(8):464-475

Seismicity induced by fluid perturbations became an important societal concern since felt earthquakes (M_w up to 6) occurred after anthropogenic activities. In order to mitigate the risks associated with undesired seismicity, as well as to be able to use the micro-seismicity as a probe for in-depth investigation of fluid-driven processes, it is of crucial importance to understand the links between seismicity, fluid pressure and flow. We have developed a series of in-situ, decameter-scale experiments of fault zone reactivation by controlled fluid injection, in order to improve the near-source geophysical and hydromechanical observations. The deployed geophysical monitoring close to the injection allows one to cover the full frequency range of the fault responses from the static deformation to the very high-frequency seismic emissions (up to 4 kHz). Here, we focus on the microseismicity (M_w ∼ –4 to –3) recorded during two fluid injection experiments in low-permeable shale and highly-fractured limestone formations. In both experiments, the spatio-temporal distribution of the seismic events, the energy balance, and the seismic velocity changes of the fractured medium show that most of the deformation does not actually emit seismic signals. The induced deformation is mainly aseismic. Based on these high-resolution multiparametric observations in the near-field, we therefore proposed a new model for injection-induced seismicity: the seismicity is not directly induced by the increasing fluid pressure, but it is rather triggered by the stress perturbations transferred from the aseismic motion caused by the injection.  相似文献   

基于多源数据约束的成矿动力学模拟——以宁芜盆地钟姑矿田典型矿床为例   总被引：1，自引：0，他引：1       下载免费PDF全文

贾 蔡 袁 峰李晓晖  胡训宇廖宝胜孙维安  张明明沈 乐 《地质科学》2018,(4):1327-1346

数据驱动是大数据分析的特点之一。成矿过程的数值模拟即采用数据驱动方式,在基本的物理化学方程及地质模型基础上,通过定量的手段,模拟不同阶段、不同条件下成矿系统的演化过程,从而能够更深入、直观的讨论各种地质要素对成矿系统的影响,验证已有的成矿理论,获得新的认识。姑山铁矿床和杨庄铁矿床是宁芜火山岩盆地钟姑矿田内的典型矿床。前人研究显示,杨庄铁矿床为典型的充填-接触交代型热液矿床,而姑山铁矿床目前在成因上仍存在争议,存在岩浆成因和热液成因等多种观点。为了更好的对比和了解姑山铁矿床和杨庄铁矿床的成因及成矿过程,本文基于多源数据的综合解译剖面,采用力-热耦合数值模拟方法对与热液矿床密切相关的扩容空间的形成机制和特点开展研究。模拟结果显示,姑山铁矿床和杨庄铁矿床可能分属不同的成因类型;随着岩体的逐渐冷却,杨庄铁矿床内岩体和地层接触带会出现明显的扩容空间,可为成矿流体的汇聚及矿质的沉淀提供良好的条件,有利于热液矿床的形成;而姑山铁矿床内岩体和地层接触带只有局部区域出现扩容空间,同时体应变数值较小且形态变化大,并不利于大规模的流体汇聚及矿质的沉淀,姑山铁矿床可能主要为岩浆成因,并受少量热液的叠加改造。  相似文献   

基于数值预报模式的乡镇温度预报方法研究     

吴春英  刘多文  钟博  蒋婷婷  于蕙箐  陈佳美  高燕 《气象与环境学报》2019,(1):108-112

利用2015年8月至2017年7月长兴岛站和交流岛站日最高气温、日最低气温实况资料,对ECMWF细网格2 m温度预报值和日本FSFE02(24 h地面形势场预报)、FSFE03(36 h地面形势场预报)进行了检验。结果表明:根据历史回归统计检验,ECMWF细网格模式24 h的2 m最高气温、最低气温预报效果显著,通过了0. 05信度显著性检验。对各月做相关分析,相关性均较好。利用前一日ECMWF细网格2 m温度预报值与长兴岛站实况差值,根据统计的ECMWF细网格2 m温度预报订正值,做出长兴岛站未来24 h的气温预报。交流岛站温度预报是在长兴岛站温度预报的基础上订正做出,经统计分析,交流岛站和长兴岛站的气温差值与地面形势场和风场有较好的对应关系,根据不同类型的地面形势场和风场订正值,做出交流岛站的温度预报。应用Matlab计算机语言的开发功能,提取ECMWF细网格2 m温度预报的最高、最低气温值,并录入当日长兴岛站和交流岛站最高、最低气温实况值,自动预报各站未来24 h最高气温、最低气温。创建可视化预报工作界面,实现乡镇温度预报自动化。  相似文献   

Numerical and experimental study of a flapping foil generator     

《Applied Ocean Research》2017

The energy extraction performance of a flapping foil generator is studied through experiment and numerical simulation. A practical flapping foil generator has been proposed. The heave motion and the pitch motion of the foil are adjusted through a crankshaft-like structure. The heave and pitch motions of the foil are transferred to the rotational motion of the main shaft. A pair of gears is adopted to increase the pitch angle. A prototype with pitch amplitude θ₀ = 60^∘ has been built and the experiment is carried out in a tunnel. The overall performance of the mechanism has been analysed. Good agreement of numerical results and experiment data has been found. Further simulations with larger pitch amplitudes are carried out. It is found that higher efficiency can be achieved with larger pitch amplitude at medium frequency.  相似文献   

Stress‐dependent hardening and failure surfaces of dry sand     

W. Ehlers  O. Avci 《国际地质力学数值与分析法杂志》2013,37(8):787-809

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

A survey of numerical methods for stochastic differential equations     

P. E. Kloeden  E. Platen 《Stochastic Hydrology and Hydraulics》1989,3(3):155-178

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.  相似文献