A new multidimensional scaling (MS) technique, referred to as the Pijk model, is formulated on the basis of associations among triple objects (samples or variables), instead of pairs of objects as used in the usual MS methods, such as factor analysis. The computational scheme provided for this method is the reduction of an original problem to a standard eigenvalue-eigenvector problem. The major goal of the technique is simplification and reduction of data structures and the rescaling of original objects into a new and reduced space, so that patterns and relations of the original objects can be conventiently examined in two-dimensional factor plots. The Pïjk method is illustrated and tested by using a set of geochemical data related to the epithermal gold and silver vein deposits in the Walker Lake quadrangle of Nevada and California. The characteristics of element associations suggested in the Pijk analysis are consistent with field observations. A preliminary comparison between the new method and the ordinary factor analysis also is made on the basis of the same data set. Results are encouraging in that analysis by the Pijk model captures triple-object associations that might be missed by the ordinary factor analysis which considers only pair-variable correlations 相似文献
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. 相似文献
Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments. 相似文献
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. 相似文献
Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure
of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile
passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust
(MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out
of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric
strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented,
which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate,
(iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the
region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that
the major features of the geoelectrical model are relevant and desired by the MT data. 相似文献
We monitored near-surface atmospheric fallout (15-cm above ground) and soil solution (at 15, 35 and 55 cm below ground) derived nanoparticles over an 8-month period by collecting the particles directly onto TEM grids in anthropogenically-influenced (vineyard) and pristine (native forest) sites in France. Particle clusters trapped on the grid were selected randomly and individual particles were binned into eight different groups (euhedral clays, weathered clays, salts, oxi-hydroxides, bacteria, non-living organic matter, aggregates and undetermined). Bacteria represent 9–23% of the collected nanoparticle area (ave. 9.4% and 18% for two atmospheric collection sites and ave. 23% for soil infiltration samples). Bacteria were most often associated with non-living organic matter and comprised a variety of morpho-types. Interestingly, 45% of all the bacteria analyzed by transmission electron microscopy and electron dispersive spectroscopy (TEM-EDX) showed the presence of intracellular grains significantly enriched in lead and phosphorus. Intracellular sequestration of Pb into polyphosphate bodies has been observed in the laboratory, but this is the first observation of this phenomenon in a natural environment. Furthermore, this suggests that microbial-bound Pb may be an important transport mechanism in subsurface environments. 相似文献
The backward particle tracking method, an effective and powerful tool that can be used to delineate groundwater protection
zones, is presented. The theoretical background and insights on the applicability of this method are provided. Moreover, the
present work enriches the backward particle tracking method with an uncertainty analysis concerning the porosity values, applying
a Monte Carlo (MC) approach, coupled with the use of geographical information systems (GIS). As an application example, a
wellfield in the Komotini area, Greece, is investigated. The present study may serve as a potential guideline for wellfield
delineation, particularly in areas like Greece where lack of data related to the hydrogeological system is often a problem. 相似文献