In this paper, we present the finite cube elements method (FCEM); a novel numerical tool for calculating the gravity anomaly g and structural index SI of solid models with defined boundaries and variable density distributions, tilted or in normal position (e.g. blocks, faulted blocks, cylinders, spheres, hemispheres, triaxial ellipsoids). Extending the calculation to fractal objects, such as Menger sponges of different orders and bodies defined by polyhedrons, demonstrates the robustness of FCEM. In addition, approximating the cube element by a sphere of equal volume makes the calculation of gravitation and related derivatives much simpler. In gravity modelling of a sphere, cubes with edges of 100 m and 200 m achieve a good compromise between running time and overall error. Displaying the distribution of SI of the studied models on contour maps and profiles will have a strong impact on the forward and inverse modelling of potential field data, especially for Euler deconvolution. For Menger sponges, plots of gravity elements g and its derivatives show similar patterns independent of fractal order. Moreover, both the pattern and magnitude of SI are independent of fractal order, allowing the use of SI as a new invariant measure for fractal objects. However, SI pattern and magnitude strongly depend on the depth to the buried bodies as do other elements In this study, we also present a new type of plot; the structural index against distance variation diagrams from which we extract the three critical SI ( CSI ) values, one per axis. The inversion of gravity anomaly data at CSI values gives the optimal mean location of the buried body. 相似文献
The cumulative probability distributions for stream order, stream length, contributing area, and energy dissipation per unit length of channel are derived, for an ordered drainage system, from Horton's laws of network composition. It is shown how these distributions can be related to the fractal nature of single rivers and river networks. Finally, it is shown that the structure proposed here for these probability distributions is able to fit the observed frequency distributions, and their deviations from straight lines in a log-log plot. 相似文献
Examples of situations are presented where the grading of a soil changes during its lifetime either by crushing of particles
leading to an increase of fine material or by slow transport of fine particles with seepage leading to a decrease of fine
material. Such grading changes influence the basic constitutive properties of the soil, in particular properties such as critical
states which are dependent on the available range of densities of packing. Discrete element modelling is used to show the
dependence of critical state conditions on grading and the way in which the particle assembly seeks out new critical state
conditions as the grading changes. 相似文献
Upper Pollara eruption products (13 ka, Salina Island, Italy) include both homogeneous and heterogeneous pumices resulting from mixing/mingling processes between an HK andesite and a high-SiO2 rhyolite. Representative samples of heterogeneous pumices are collected and analyzed in order to check the correspondence between glass composition and morphological features of the mingling/mixing structures. Image analysis techniques are applied and eight grey color ranges (classes) are extracted from high-resolution scans of pumice. Class 1 (lighter colors) and class 8 (darker colors) show end-member glass compositions, i.e. HK andesite and high-SiO2 rhyolite, respectively. These two classes show spot- to cluster-like morphological structures. Intermediate classes show an HK dacitic to rhyolitic composition and a banding- to fold-like morphology. Fractal analysis by box-counting of the boundary pattern of eight grey classified images is performed over a length scale of 0.028–1.8 cm. Fractal dimension D is between 1.01 and 1.84. Coupled fractal analysis and geochemical data reveal that D increases as the degree of magma interaction (homogenization) increases. This feature well fits the results from numerical models on the convective mixing of fluids driven by thermal convection. We conclude that the increase of D observed in the Upper Pollara samples reflects the transition from fractal mixing to homogenization. End-member magmas (HK andesite and high-SiO2 rhyolite) represent isolated mixing regions, while homogenized magmas represent active mixing regions. In the analyzed pumices, isolated and active mixing regions coexist at scales between 10−4 and 10−2 m. Morphological and compositional features of the Upper Pollara pumices result from turbulence. 相似文献
Magma mixing structures from the lava flow of Lesbos (Greece) are analyzed in three dimensions using a technique that, starting from the serial sections of rock cubes, allows the reconstruction of the spatial distribution of magmas inside rocks. Two main kinds of coexisting structures are observed: (i) “active regions” (AR) in which magmas mix intimately generating wide contact surfaces and (ii) “coherent regions” (CR) of more mafic magma that have a globular shape and do not show large deformations. The intensity of mingling is quantified by calculating both the interfacial area (IA) between interacting magmas and the fractal dimension of the reconstructed structures. Results show that the fractal dimension is linearly correlated with the logarithm of interfacial area allowing discrimination among different intensities of mingling.
The process of mingling of magmas is simulated using a three-dimensional chaotic dynamical system consisting of stretching and folding processes. The intensity of mingling is measured by calculating the interfacial area between interacting magmas and the fractal dimension, as for natural magma mixing structures. Results suggest that, as in the natural case, the fractal dimension is linearly correlated with the logarithm of the interfacial area allowing to conclude that magma mixing can be regarded as a chaotic process.
Since chemical exchange and physical dispersion of one magma inside another by stretching and folding are closely related, we performed coupled numerical simulations of chaotic advection and chemical diffusion in three dimensions. Our analysis reveals the occurrence in the same system of “active mixing regions” and “coherent regions” analogous to those observed in nature. We will show that the dynamic processes are able to generate magmas with wide spatial heterogeneity related to the occurrence of magmatic enclaves inside host rocks in both plutonic and volcanic environments. 相似文献
In order to model flow and transport in fractured rocks it is important to know the geometry of the fracture network. A stochastic approach is commonly used to generate a synthetic fracture network from the statistics measured at a natural fracture network. The approach presented herein is able to incorporate the structures found in a natural fracture network into the synthetic fracture network. These synthetic fracture networks are the images generated by Iterated Function Systems (IFS) as introduced by Barnsley (1988). The conditions these IFS have to fulfil to determine images resembling fracture networks and the effects of their parameters on the images are discussed. It is possible to define the parameters of the IFS in order to generate some properties of a fracture network. The image of an IFS consists of many single points and has to be suitably processed for further use. 相似文献
The distributions of contact areas in single, natural fractures in quartz monzonite (Stripa granite) are found to have fractal dimensions which decrease fromD=2.00 to values nearD=1.96 as stress normal to the fractures is increased from 3 MPa up to 85 MPa. The effect of stress on fluid flow is studied in the same samples. Fluid transport through a fracture depends on two properties of the fracture void space geometry. the void aperture; and the tortuosity of the flow paths, determined through the distribution of contact area. Each of these quantities change under stress and contribute to changes observed in the flow rate. A general flow law is presented which separates these different effects. The effects of tortuosity on flow are largely governed by the proximity of the flow path distribution to a percolation threshold. A fractal model of correlated continuum percolation is presented which quantitatively reproduces the flow path geometries. The fractal dimension in this model is fit to the measured fractal dimensions of the flow systems to determine how far the flow systems are above the percolation threshold. 相似文献