This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological. 相似文献
Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper. 相似文献
We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.
We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.
The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.
ZS9 is made out of 36 zones where earthquakes with Mw > = 5 are expected. It also assumes that earthquakes with Mw up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.
Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates. 相似文献
MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because
of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for
water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing
with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and
domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development
of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology.
Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing
of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship
between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater
flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater
resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water
budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December
2003. During this period, the total recharge of the groundwater system was 49,374 × 106 m3 and the total discharge was 56,530 × 106 m3 the budget deficit was −7,156 × 106 m3. In this integrated system, the original data including graphs and attribution data could be stored in the database. When
the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core
program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online. 相似文献
We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock
along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for
component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition,
THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used
to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled
for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and
growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G
produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium
phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a
pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the
combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion.
Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.
Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors
governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts
were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using
high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological
methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural
analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast
to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast
size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast
size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization
annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust
pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these
calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism
bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures
for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens,
respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet
crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations
of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal
size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily
to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important
in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites. 相似文献
Flood mapping requires the combination and integration of geomorphological and hydrological-hydraulic methods; however, despite
this, there is very little scientific literature that compares and validates both methods. Two types of analysis are addressed
in the present article. On the one hand, maps of flood plains have been elaborated using geomorphological evidence and historical
flood data in the mountainous area of northwestern Spain, covering an area of more then 232 km2 of floodplains. On the other hand, a hydrometeorological model has been developed (Clark semidistributed unit hydrograph)
in the Sarria River basin (155 km2, NW Spain). This basin is not gauged, hence the model was subjected to a goodness-of-fit test of its parameter (curve number)
by means of Monte Carlo simulation. The peak flows obtained by means of the hydrological model were used for hydraulic modeling
(one-phase, one-dimensional and steady flow) in a 4 km2 urban stretch of the river bed. The delineation of surface areas affected by floods since 1918, as well as those analyzed
subsequent to the geomorphological study, reveals a high degree of reliability in the delineation of the flooded areas with
frequent recurrence intervals (<50 years). If we compare these flooded surface areas with the estimate obtained by the hydrological-hydraulic
method we can see that the latter method overestimates the extent of the surface water by 144% for very frequent recurrence
intervals (>10 years) and underestimates it as the recurrence interval increases, by up to 80% less floodplain for exceptional
events (>500 years). Finally, a management map is put forth combining the most reliable results available by integrating both
methods.
Originally presented at the Sixth International Conference on Geomorphology. 相似文献