Multifractal modeling and spatial statistics

In general, the multifractal model provides more information about measurements on spatial objects than a fractal model. It also results in mathematical equations for the covariance function and semivariogram in spatial statistics which are determined primarily by the second-order mass exponent. However, these equations can be approximated by power-law relations which are comparable directly to equations based on fractal modeling. The multifractal approach is used to describe the underlying spatial structure of De Wijs 's example of zinc values from a sphalerite-bearing quartz vein near Pulacayo, Bolivia. It is shown that these data are multifractal instead of fractal, and that the second-order mass exponent (=0.979±0.011 for the example) can be used in spatial statistical analysis.     

Discrete multifractals

The concept of multifractal modeling has been used intensively in various fields of science for characterizing measures with self- similarity. It has been shown that multifractal modeling provides powerful tools for characterizing patterns in the spatial distribution of geological quantities and objects. Existing multifractal models were proposed for the purpose of handling spatially intertwined fractals with continuous fractal spectrum f(α) (or continuous codimension function C(γ)). In this paper, these conventional multifractals are termed “continuous multifractals” whereas multifractals with discrete fractal dimensions are termed “discrete multifractals.” The properties of discrete multifractals are investigated. It is shown by various artificial examples and a case study of stratigraphy of Ocean Drilling Program (ODP) site 645 that spatially intertwined fractals/multifractals indeed can have discrete fractal dimensions. Histogram-and moment-based techniques are proposed for discrete multifractal modeling and illustrated using the artificial examples. The new concept of discrete multifractals and associated multifractal modeling yields not only techniques for characterizing multifractals with discrete fractal dimensions but it also provides insight into the relationships between fractals, bifractals, and multifractals.     

Characterization of Geochemical Distributions Using Multifractal Models

The use of multifractals in the applied sciences has proven useful in the characterization and modeling of complex phenomena. Multifractal theory has also been recently applied to the study and characterization of geochemical distributions, and its relation to spatial statistics clearly stated. The present paper proposes a two-dimensional multifractal model based on a trinomial multiplicative cascade as a proxy to some geochemical distribution. The equations for the generalized dimensions, mass exponent, coarse Lipschitz–Hölder exponent, and multifractal spectrum are derived. This model was tested with an example data set used for geochemical exploration of gold deposits in Northwest Portugal. The element used was arsenic because a large number of sample assays were below detection limit for gold. Arsenic, however, has a positive correlation with gold, and the two generations of arsenopyrite identified in the gold quartz veins are consistent with different mineralizing events, which gave rise to different gold grades. Performing the multifractal analysis has shown problems arising in the subdivision of the area with boxes of constant side length and in the uncertainty the edge effects produce in the experimental estimation of the mass exponent. However, it was possible to closely fit a multifractal spectrum to the data with enrichment factors in the range 2.4–2.6 and constant K₁ = 1.3. Such parameters may give some information on the magnitude of the concentration efficiency and heterogeneity of the distribution of arsenic in the mineralized structures. In a simple test with estimated points using ordinary lognormal kriging, the fitted multifractal model showed the magnitude of smoothing in estimated data. Therefore, it is concluded that multifractal models may be useful in the stochastic simulation of geochemical distributions.     

川西北地区控矿断裂的分形特征及其预测意义

断裂构造是控制川西北地区金矿形成与分布的主导因素。运用分形方法定量计算了该区金矿化异常和控矿断裂体系的计盒维数和信息维数。研究结果表明 ,利用地球化学数据二维空间序列的R/S分形方法可以较为精确有效地厘定金矿化的空间结构特征 ,分数维D值的大小能表征控矿断裂体系的复杂性。金矿化带整体上受NW向断裂控制 ,但金矿体则产出于NW向断裂与NE向断裂的复合部位。根据对三个矿化区断裂体系分维特征与金矿发育特征关系的分析 ,发现断裂体系的分维高值区与金矿分布密集区对应。阿坝地块西南缘的找矿前景不如东北缘     

中国金矿床分布的分形研究

根据1:40O万中国岩金成矿图,在矿床地质研究基础上,用分形理论探讨了中国岩金矿的空间分布。研究表明,在两个尺度范围内,金矿床在空间上呈分形分布,其中数盒子法的分维数D分别为0.3333(5～80km)和1.3259(80～2400km);密度分布的分维数D分别为1.2033(5～80km)和1.5459(80～2400km)。金矿床的分形分布给超大陆旋回会聚构造边缘控制金矿的分布及成矿作用的观点以有力的佐证。     

Correlation Fractal and Multifractal Characteristics of Seismic Activity in the Taiwan Area, China

Based on the analysis of newly collected data of plate tectonics, distribution of active faults and crustal deformation, the Taiwan area is divided into two seismic regions and six seismic belts. Then, correlation fractal dimensions of all the regions and belts are calculated, and the fractal characteristics of hypocenteral distribution can be quantitatively analyzed. Finally, multifractal dimensions Dq and f(α) are calculated by using the earthquake catalog of the past 11 years in the Taiwan area. This study indicates that (1) there exists a favorable corresponding relationship between spatial images of seismic activity described with correlation fractal dimension analysis and tectonic settings; (2) the temporal structure of earthquakes is not single but multifractal fractal, and the pattern of Dq variation with time is a good indicator for predicting strong earthquake events.     

Multifractal Modeling and Lacunarity Analysis

The so-called gliding box method of lacunarity analysis has been investigated for implementing multifractal modeling in comparison with the ordinary box-counting method. Newly derived results show that the lacunarity index is associated with the dimension (codimension) of fractal, multifractal and some types of nonfractals in power-law relations involving box size; the exponent of the lacunarity function corresponds to the fractal codimension (E – D) for fractals and nonfractals, and to the correlation codimension (E – lpar;2)) for multifractals. These results are illustrated with two case studies: De Wijs's zinc concentration values from the Pulacayo sphalerite-quartz vein in Bolivia and Cochran's tree seedlings example. Both yield low lacunarities and slightly depart from translational invariance.     

Markov Processes and Discrete Multifractals

Fractals and multifractals are a natural consequence of self-similarity resulting from scale-independent processes. Multifractals are spatially intertwined fractals which can be further grouped into two classes according to the characteristics of their fractal dimension spectra: continuous and discrete multifractals. The concept of multifractals emphasizes spatial associations between fractals and fractal spectra. Distinguishing discrete multifractals from continuous multifractals makes it possible to describe discrete physical processes from a multifractal point of view. It is shown that multiplicative cascade processes can generate continuous multifractals and that Markov processes result in discrete multifractals. The latter result provides not only theoretical evidence for existence of discrete multifractals but also a fundamental model illustrating the general properties of discrete multifractals. Classical prefractal examples are used to show how asymmetrical Markov process can be applied to generate prefractal sets and discrete multifractals. The discrete multifractal model based on Markov processes was applied to a dataset of gold deposits in the Great Basin, Nevada, USA. The gold deposits were regarded as discrete multifractals consisting of three spatially interrelated point sets (small, medium, and large deposits) yielding fractal dimensions of 0.541 for the small deposits (<25 tons Au), 0.296 for the medium deposits (25--500 tons Au), and 0.09 for the large deposits (>500 tons Au), respectively.     

Discrete multifractals

The concept of multifractal modeling has been used intensively in various fields of science for characterizing measures with self- similarity. It has been shown that multifractal modeling provides powerful tools for characterizing patterns in the spatial distribution of geological quantities and objects. Existing multifractal models were proposed for the purpose of handling spatially intertwined fractals with continuous fractal spectrum f(α) (or continuous codimension function C(γ)). In this paper, these conventional multifractals are termed “continuous multifractals” whereas multifractals with discrete fractal dimensions are termed “discrete multifractals.” The properties of discrete multifractals are investigated. It is shown by various artificial examples and a case study of stratigraphy of Ocean Drilling Program (ODP) site 645 that spatially intertwined fractals/multifractals indeed can have discrete fractal dimensions. Histogram-and moment-based techniques are proposed for discrete multifractal modeling and illustrated using the artificial examples. The new concept of discrete multifractals and associated multifractal modeling yields not only techniques for characterizing multifractals with discrete fractal dimensions but it also provides insight into the relationships between fractals, bifractals, and multifractals.     

空间模式的广义自相似性分析与矿产资源评价

尺度不变性(scale invariance)包括自相似性(各向同性)、自仿射性(成层结构)、广义自相似性(各向异性标度不变性),是由各种地质过程和地质事件所产生的地质特征和模式的本质属性．尺度不变性可用分形和多重分形模型来表征．这些尺度特征的定量化可为刻画地质空问模式和模式识别提供有力的工具．例如。热液矿床的群聚现象可以用局部分形特征(局部奇异性)来刻画．通过在特征空问中(如频率空问)识别空问模式的广义自相似性．可以将空间混合模式进行分解或异常的识别．介绍了几种相关的分形模型和方法。包括度量空问模式广义尺度独立性(GSI)的线性模型;基于广义尺度独立性的异常分解S—A方法;度量空问模式的局部奇异性方法;以及如何利用分形特征预测未发现矿床的2种方法．有些方法已应用于许多矿产资源评价实例中．给出了对加拿大Nova Scotia省西南部湖泊沉积物样品中的4种元素As、Pb、Zn和Cu的地球化学数据处理分析结果。证明了局部奇异性分析和S—A异常分解方法对地球化学异常的增强和分离的有效性．研究表明：由S—A方法分解的异常往往具有多重分形的特点,而且普遍具有局部奇异性．研究区内具有明显奇异性的地区(元素含量富集区)是金矿异常区域。它们与金矿成矿作用和已知矿床的赋存密切相关．     

地层含金性评价的分数维研究———以琼西抱板群为例

摘 要　 利用分形理论研究了琼西抱板群金含量分布的分数维结构特征‚并对其含金性作出 评价。结果表明‚较小的分维及具有多标度分形‚有利于金矿化的产出。     

四川地区断层空间分布的多重分形特征

断层系统多重分形避免了容量维Ｄ０的不足，考虑了每一条断层对整体分维的贡献，反映了断层的空间分布特征。根据四川地区断层分布图，运用多重分形理论检验了断层空间分布的自相似性，分区计算了四川地区各带的多重分维值。结果表明，断层破裂系统为分形结构，具有很好的自相似性，相关系数均达０．９９。不同带分维值的较大差异，反映不同构造区域的断层分布特点。最后，对分维的含义进行了深入探讨，对油气勘探的有利区域进行了分析。     

Multifractality as a Measure of Spatial Distribution of Geochemical Patterns

To quantify the spatial distribution of geochemical elements, the multifractality indices for Zn, Cu, Pt, Pd, Cr, Ni, Co, Pb, and As in lake-sediment samples in the Shining Tree area in the Abitibi area of Ontario are determined. The characterization of multifractal distribution patterns is based on the box-counting moment method and involves three functions: a mass exponent function (q); Coarse Hölder Exponent (q); and fractal dimension spectrum f( (q)). Properties of these functions at different values of q, characterize the spatial distribution of the variable under study. It is shown that the degree of multifractality defined by (1) can be used as a measure of irregularity of geochemical spatial dispersion patterns. The variations of Zn and Cu in the study area are characterized by relatively low degree of multifractality, whereas those for Pt, Pd, Cr, Ni, and Co; and particularly for As and Pb are characterized by higher multifractality indices.In the case of Zn and Cu, singularity spectra are close to a monofractal compared to the ones for As an Pb. The determination of multifractality indices allows us, in a quantitative way, to study the pattern of metal dispersions and link them to different physical processes, such as metal adsorption by organic material or glaciogenic processes.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

基于分形理论的成矿空间分析 ——以新疆东天山康古尔塔格金矿带为例

利用分形理论中的盒式分维探究新疆东天山康古尔塔格金矿带构造及金矿床的空间分布,得出两者在空间上具有双重分形的特点。结合前人的研究成果,认为研究区在区域上是韧性的剪切带,而在局部体现为脆性地质体,金矿的分布在一定范围是随机的,但在区域上受韧性剪切带的控制,呈带状或线状分布,更为重要的是线性构造与金矿可能是同一时期同一构造运动的结果。在区域上,脆韧性剪切转换带是金矿的有利成矿部位。     

滑坡位移多重分形特征与滑坡演化预测

在系统分析滑坡位移监测资料和位移演化特征的基础上,根据多重分形理论基本原理,对滑坡位移演化所具有的复杂性、突变性和非线性特征进行了分析和研究。单一分形维数对滑坡位移的演化趋势预测存在不足,文中分别以新滩滑坡、丹巴滑坡和黄蜡石滑坡为例,计算了滑坡位移时序演化的多重分维数演化特征。分析和评价位移演化规律与多重分维数演化特征的关系发现,多重分维数D1 > D2 > … > D∞时,滑坡趋于稳定;D1 < D2 < … < D∞时,滑坡向失稳破坏演化。当滑坡位移时序多重分维数演化特征出现拐点时,即分维数由D1 > D2 > … > D∞,经D1 > … > Dn < Dn+1 < … < D∞到 D1 < D2 < … < D∞的演化过程时,滑坡向不稳定的状态演化;当分维数由D1 < D2 < … < D∞,经D1 > … > Dn D2 > … > D∞的演化过程时,滑坡向趋于稳定的状态演化。研究表明,可以运用多重分维数演化特征对滑坡位移演化趋势与规律进行评价与预测。     

Elasticity of single-crystal aragonite by Brillouin spectroscopy

The elastic constants of natural single-crystal aragonite (CaCO₃) have been measured by Brillouin spectroscopy at ambient conditions. The elastic constants C₁₁, C₂₂, C₃₃, C₄₄, C₅₅, C₆₆, C₁₂, C₁₃ and C₂₃ are 171.1±1.0, 110.1±0.9, 98.4±1.2, 39.3±0.6, 24.2±0.4, 40.2±0.6, 60.3±1.0, 27.8±1.6 and 41.9±2.0 GPa, respectively, for aragonite. The linear compressibilities of the a-, b- and c-axis for aragonite at ambient conditions were derived from our measured data to be 3.0±0.2, 4.2±0.2 and 7.3±0.6×10^–3 GPa^–1, respectively. The aggregate bulk and shear moduli for aragonite using the Voigt-Reuss-Hill (VRH) scheme are thus calculated to be 68.9±1.4 and 35.8±0.2 GPa, respectively. The value of bulk modulus is in remarkable contrast to the literature value of 46.9 GPa measured almost a century ago. Our new datum, however, is closer to that derived from recent atomistic simulation and static compression studies.     

新疆东准噶尔地区断裂的多重分形研究

结合东准噶尔地区强应变构造带的空间展布情况和断裂的空间分布特征 ,用多重分形模型测算了各带的多重分维值。研究结果表明 :强应变构造带在不同尺度上具有自相似性 ;断裂的空间分布为多重分形结构 ,断裂的演化有从复杂几何结构的次级断裂组合向单一连续的大断裂过渡而分维值逐渐减小的趋势 ;不同的大地构造单元有不同的分维值 ,而且有靠近板块缝合带 (或海沟 )分维值减小、远离板块缝合带 (或海沟 )分维值增大的趋势     

Multifractality in seismic sequences of NW Himalaya

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.

     

The fractal properties of geochemical landscapes as an indicator of weathering and transport processes within the Eastern Alps

Variations in surface morphology and lithology provide an opportunity to study lithologic and morphologic influences on the spatial pattern of stream-sediment geochemistry within two contrasting environments of the Eastern Alps (Hohe Tauern Range and Gurktaler Alpen Range). The fractal dimension, a measure of surface roughness over a variety of scales, is used to model the dissipation of erosive products due to climatic controlled denudation and fluvial mass transport. Based on a spatial correlation analysis, specific elemental concentrations are used as indicators for a dominant lithotype. Fractal geometry of these elements has been estimated by sequential Gaussian simulation of the area/perimeter relationship (Dal) and by the estimation of multifractal spectra. It is shown that within a 510–780 km² survey area the spatial variations of Al, Ga, Ni and Ca can be approximated by single fractals but for those of Ag and Sn multifractal models must be used. Fractal properties derived from simulated surfaces are explainable by the process controlling the spatial structure of the data. Climatic and tectonic parameters apparently influences Dal at large scales. At smaller scales rock-type variation exert an additional influence on Dal.