基于GIS城市洪水淹没模拟分析

城市洪水淹没范围动态模拟分析是城市防洪规划与防治决策的一个重要基础工作。本文主要研究洪水从翻堤口出发在地理空间蔓延、扩散、动态行进及确定淹没范围的数字模拟方法。为此,根据数学形态学及测地圆概念,研究设计了洪水扩散范围的“膨胀”模拟算法和淹没范围搜索算法,用于查询淹没通块中从翻堤处到任一点之间的淹没路径和ｔ时刻洪水扩散范围。     

The effects of carbon dioxide leakage on fractures, and water quality of potable aquifers during geological sequestration of CO2

Since industrial revolution, the ＂greenhouse effect＂ is one of the most important global environmental issues. Of all the greenhouse gases, CO2 is responsible for about 64% of the enhanced ＂greenhouse effect＂, making it the target for mitigation, so reducing anthropogenic discharge of carbon dioxide attracts more and more attention. Geological sequestration of CO2 in deep saline aquifers is one of the most promising options. But because unknown fractures and faults may exist in the caprock layers which can prevent the leakage of CO2, CO2 will leak upward into upper potable aquifers, and lead to adverse impacts on the shallow potable aquifers. In order to assess the potential effect of CO2 leakage from underground storage reservoirs on fractures and water quality of potable aquifers, this study used the non-isothermal reactive geochemical transport code TOUGHREACT developed by Xu et al to establish a simplified 2-D model of CO2 underground sequestration system, which includes deep saline aquifers, caprock layers, and shallow potable aquifers, and study and analyze the changes of mineral and aqueous components. The simulation results indicated that the minerals of deep saline aquifers and fractures should be mainly composed of aluminosilicate and silicate minerals, which not only enhance the mass of CO2 sequestrated by mineral trapping, but also decrease the porosity and permeability of caprock layers and fractures to prevent and reduce CO2 leakage. The results from deep saline aquifers showed that the mass of carbon dioxide trapped by minerals and solution phases is limited, the rest remained as a supercritical phase, and so once the caprock aquifers have some unknown fractures, the free carbon dioxide phase may leak from CO2 geologic sequestration reservoirs by buoyancy.     

The bimodal rift-related Juscelândia volcanosedimentary sequence in central Brazil: Mesoproterozoic extension and Neoproterozoic metamorphism

The Barro Alto Complex and Juscelândia volcanosedimentary sequence are exposed in the central part of the Neoproterozoic Brasília belt of central Brazil. The former is a large (approximately 150 km long), boomerang-shaped, mafic-ultramafic, layered complex formed by two different intrusions metamorphosed under granulite facies. These rocks are tectonically overlain by rocks of the Juscelândia volcanosedimentary sequence, represented mainly by biotite-gneiss and amphibolite, or amphibolite facies metamorphic equivalents of rhyolite and basalt, respectively. New SIMS U–Pb zircon data and Sm–Nd isochron data presented herein help clarify the igneous and metamorphic evolution of the Juscelândia volcanosedimentary sequence, as well as its relationship with the Barro Alto Complex. Zircon grains from two biotite gneisses were analyzed by SIMS (SHRIMP) and indicate Mesoproterozoic dates, approximately 1.28 Ga, interpreted as the time of bimodal volcanism in a tectonic setting transitional between a continental rift and an ocean basin. Metamorphism is constrained by Sm–Nd garnet-whole-rock isochrons for garnet amphibolite and pelitic schists of the Juscelândia sequence, as well as for clinopyroxene-garnet amphibolite and garnet granulite of the Barro Alto Complex, which give ages between 0.74 and 0.76 Ga, in agreement with SIMS dates for metamorphic zircon rims. These new data are significant, because they establish that a single metamorphic event affected both the Barro Alto Complex and the Juscelândia sequence. Based on these new data, we present a modified tectonic model for the Brasília belt.     

Sequence stratigraphy of a Late Devonian ramp-situated reef system in the Western Canada Sedimentary Basin: dynamic responses to sea-level change and regressive reef development

The Alexandra Formation, located in the Northwest Territories of Canada, is formed of a Late Devonian (Frasnian) reef system that developed on a gently sloping, epicontinental ramp in the Western Canada Sedimentary Basin. High‐resolution sequence stratigraphic analysis of its deposits delineates two reef complexes that are separated by a Type I sequence boundary. The second reef complex developed on the outer ramp, basinward of the first, after sea‐level fell ≈17 m. Stratigraphic complexity of the second reef complex was a result of its initiation during forced regression, and its development through an entire cycle of sea‐level rise followed by sea‐level fall. Its highstand systems tract was not characterized by high rates of carbonate production or sediment shedding. Rather, these features took place as sea‐level fell, after its highstand systems tract. The sequence stratigraphic framework of this regressive reef system highlights a number of depositional parameters that differ from high‐relief, shelf‐situated reef systems with steep, narrow margins. These have implications for understanding the controls on the development of ramp‐situated reef systems, and the nature of reef systems with gently sloping profiles. This study demonstrates that the development of stromatoporoid reef systems may be far more complex than generally realized, and that high‐resolution sequence stratigraphy may provide the tools for better understanding of complex, often enigmatic, aspects of these systems.     

A sequence stratigraphic model for an intensely bioturbated shallow-marine sandstone: the Bridport Sand Formation, Wessex Basin, UK

The Bridport Sand Formation is an intensely bioturbated sandstone that represents part of a mixed siliciclastic‐carbonate shallow‐marine depositional system. At outcrop and in subsurface cores, conventional facies analysis was combined with ichnofabric analysis to identify facies successions bounded by a hierarchy of key stratigraphic surfaces. The geometry of these surfaces and the lateral relationships between the facies successions that they bound have been constrained locally using 3D seismic data. Facies analysis suggests that the Bridport Sand Formation represents progradation of a low‐energy, siliciclastic shoreface dominated by storm‐event beds reworked by bioturbation. The shoreface sandstones form the upper part of a thick (up to 200 m), steep (2–3°), mud‐dominated slope that extends into the underlying Down Cliff Clay. Clinoform surfaces representing the shoreface‐slope system are grouped into progradational sets. Each set contains clinoform surfaces arranged in a downstepping, offlapping manner that indicates forced‐regressive progradation, which was punctuated by flooding surfaces that are expressed in core and well‐log data. In proximal locations, progradational shoreface sandstones (corresponding to a clinoform set) are truncated by conglomerate lags containing clasts of bored, reworked shoreface sandstones, which are interpreted as marking sequence boundaries. In medial locations, progradational clinoform sets are overlain across an erosion surface by thin (<5 m) bioclastic limestones that record siliciclastic‐sediment starvation during transgression. Near the basin margins, these limestones are locally thick (>10 m) and overlie conglomerate lags at sequence boundaries. Sequence boundaries are thus interpreted as being amalgamated with overlying transgressive surfaces, to form composite erosion surfaces. In distal locations, oolitic ironstones that formed under conditions of extended physical reworking overlie composite sequence boundaries and transgressive surfaces. Over most of the Wessex Basin, clinoform sets (corresponding to high‐frequency sequences) are laterally offset, thus defining a low‐frequency sequence architecture characterized by high net siliciclastic sediment input and low net accommodation. Aggradational stacking of high‐frequency sequences occurs in fault‐bounded depocentres which had higher rates of localized tectonic subsidence.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Modeling of the classification behavior of the diaphragms used in multi-chamber cement mills

A new approach has been developed for describing the classification that arises at the diaphragm in multi-chamber ball mills. Fourteen industrial scale cement mills have been sampled in this study and after mass balancing, model fits of the equipment were done. Successive mill product predictions for the sampled conditions were achieved, and requirements to build up an optimization capable model are discussed.     

虚拟样机技术在地质工程机械领域的应用研究

基于对国内地质工程机械的研发现状的考虑，根据虚拟样机技术的特点，讨论了在地质工程机械领域实施虚拟样机技术的必要件，并提出了建立地质工程机械虚拟样机的策略。联合PRO／E、ANSYS及ADAMS建立了气动潜孔锤的虚拟样机模型，对气动潜孔锤虚拟样机进行了仿真，阐明了虚拟样机技术在地质工程机械领域的应用前景。     

Paleoenvironmental records of water level and climatic changes from the middle to late Holocene at a Lake Erie coastal wetland, Ontario, Canada

Pollen and diatom assemblages, and peat stratigraphies, from a coastal wetland on the northern shore of Lake Erie were used to analyze water level and climatic changes since the middle Holocene and their effects on wetland plant communities. Peat deposition began 4700 cal yr B.P. during the Nipissing II transgression, which was driven by isostatic rebound. At that time, a diatom-rich wild rice marsh existed at the site. Water level dropped at the end of the Nipissing rise at least 2 m within 200 yr, leading to the development of shallower-water plant communities and an environment too dry for most diatoms to persist. The sharp decline in water level was probably driven primarily by outlet incision, but climate likely played some role. The paleoecological records provide evidence for post-Nipissing century-scale transgressions occurring around 2300, 1160, 700 and 450 cal yr B.P. The chronology for these transgressions correlates with other studies from the region and implies climatic forcing. Peat inception in shallow sloughs across part of the study area around 700 cal yr B.P. coincides with the Little Ice Age. These records, considered alongside others from the region, suggest that the Little Ice Age may have resulted in a wetter climate across the eastern Great Lakes region.     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.