Computation of surface energies in an electrostatic point charge model: II. Application to zircon (ZrSiO4)

The attachment energies, the slice energies and the specific surface energies can be calculated in an electrostatic point charge model using the formula derived by Madelung for the potential introduced by an infinite row of equally spaced point charges. Power series are given for the Hankel function iH ₍₀₎ ⁽¹⁾ (iy) and (x)=d ln x!/dx. The logarithmic expression in the Madelung formula converges rapidly when applying a power series, which combines equally charged cations and anions. Besides the specific surface energy ( ^hkl), the slice energy (E _s ^hkl ) and the attachment energy (E _a ^hkl ) can be considered as special categories of surface energies as they depend on surface configurations as well. The specific surface energy is the energy per unit area of surface needed to split the crystal parallel to a face (hkl). The attachment energy (E _a) is the energy released per mole, when a new slice of thickness d _hkl crystallizes on an already existing crystal face (hkl). The growth rate of the crystal face (hkl) is a function of its attachment energy. The slice energy (E _s) is the energy released per mole, when a new slice d _hkl is formed from the vapour neglecting the influence of edge energies. The lattice energy (E _c) which is the energy released per mole of a crystal crystallizing from the vapour, is given by the following relation: E _c=E _a+E _s.     

渤海湾盆地东营凹陷现河矿区古近系断裂体系及形成演化

受双重动力背景控制,渤海湾盆地东营凹陷现河矿区古近系伸展与走滑作用叠加明显,明确断裂体系叠加关系对油气勘探意义重大。文章运用高精度地震资料及钻井资料,通过井震结合解释、构造样式识别、断裂活动性定量表征及平衡剖面恢复等技术手段,从断裂类型与组合样式、断裂叠加关系、断裂活动及演化、区域动力学背景等方面开展了深入研究。结果表明:现河矿区古近系发育伸展与伸展-走滑两套断裂体系,伸展-走滑断裂基于先期伸展断裂叠加走滑作用形成。古近系主干断裂存在沙三中期和东营早期两次活动高峰,受控于伸展-走滑应力叠合作用,现河矿区古近系共经历初始伸展阶段（Ek—Es₄）、强伸展阶段（Es₃下—Es₂下）、走滑-弱伸展阶段（Es₂上—Es₁）、再活化伸展阶段（Ed）和定型保留阶段（N—Q）5个主要的断裂演化阶段。      

Computation of surface energies in an electrostatic point charge model: I. Theory

The attachment energies, the slice energies and the specific surface energies can be calculated in an electrostatic point charge model using the formula derived by Madelung for the potential introduced by an infinite row of equally spaced point charges. Power series are given for the Hankel function iH ₍₀₎ ⁽¹⁾ (iy) and Ψ(x)=d ln x!/dx. The logarithmic expression in the Madelung formula converges rapidly when applying a power series, which combines equally charged cations and anions. Besides the specific surface energy (γ ^hkl), the slice energy (E _s ^hkl ) and the attachment energy (E _a ^hkl ) can be considered as special categories of surface energies as they depend on surface configurations as well. The specific surface energy γ is the energy per unit area of surface needed to split the crystal parallel to a face (hkl). The attachment energy (E _a) is the energy released per mole, when a new slice of thickness d _hkl crystallizes on an already existing crystal face (hkl). The growth rate of the crystal face (hkl) is a function of its attachment energy. The slice energy (E _s) is the energy released per mole, when a new slice d _hkl is formed from the vapour neglecting the influence of edge energies. The lattice energy (E _c) which is the energy released per mole of a crystal crystallizing from the vapour, is given by the following relation: E _c=E _a+E _s.     

Effect of the form of data on the quality of mine tremors hazard forecasting using neural networks

During hard coal mining operations conducted under conditions of rockburst hazard, one of the most important preventive measures can be the prediction of occurrence time and location of the strong seismic mine tremors of energy E _s ⩾ 10⁴ J. This is a very difficult task and the way it is being currently performed appears to be unsatisfactory. Therefore, attempts have been made to use neural networks, specifically trained for this application. The paper presents an approach for determining an influence of the type and shape of the input data on the efficiency of such a prediction. The considerations are based on a selected example of the seismic activity recorded during longwall mining operations conducted in one of the Polish mines.     

Delamination and ultra‐deep subduction of continental crust: constraints from elastic wave velocity and density measurement in ultrahigh‐pressure metamorphic rocks

Thirty‐three samples, including 22 eclogites, collected from the Dabie ultrahigh‐pressure (UHP) metamorphic belt in eastern China, have been studied for seismic properties. Compressional (V_p) and shear wave (V_s) velocities in three mutually perpendicular directions under hydrostatic pressures up to 1.0 GPa were measured for each sample. At 1.0 GPa, V_p (7.5–8.4 km s^?1), V_s (4.2–4.8 km s^?1), and densities (3.2–3.6 g cm^?3) in the UHP eclogites are higher than those of UHP orthopyroxenite (7.3–7.5 km s^?1, 4.1–4.3 km s^?1, 3.2–3.3 g cm^?3, respectively) and HP eclogites (7.1–7.9 km s^?1, 4.0–4.5 km s^?1, 3.1–3.5 g cm^?3, respectively). Kyanitites (with 99.5% kyanite) show extremely high velocities and density (9.37 km s^?1, 5.437 km s^?1, 3.581 g cm^?3, respectively). The eclogites show variation of V_p‐ and V_s‐anisotropy up to 9.70% and 9.17%, respectively. Poisson’s ratio (σ) ranges from 0.218 to 0.278 (with a mean of 0.255) for eclogites, 0.281–0.298 for granulites and 0.248 to 0.255 for amphibolites. The σ values for serpentinite (0.341) and marble (0.321) are higher than for other lithologies. The elastic moduli K, G, E of kyanitite were obtained as 163, 102 and 253 GPa, respectively. The V_p and density of representative UHP metamorphic rocks (eclogite & kyanitite) were extrapolated to mantle depth (15 GPa) following a reasonable geotherm, and compared to the one dimension mantle velocity and density model. The comparison shows that V_p and density in eclogite and kyanitite are greater than those of the ambient mantle, with differences of up to ΔV_p > 0.3 km s^?1 and Δρ > 0.3–0.4 g cm^?3, respectively. This result favours the density‐induced delamination model and also provides evidence in support of distinguishing subducted high velocity materials in the upper mantle by means of seismic tomography. Such ultra‐deep subduction and delamination processes have been recognized by seismic tomography and geochemical tracing in the postcollisional magmatism in the Dabie region.     

原子折射度与电子构型的关系

原子折射度与其电子构型密切相关。元素在周期表中按最外层电子数可分为s、p、sd、df4个区,原子折射度的变化在各个区内是不同的。原子折射度可表示为主量子数及最外层电子数的函数,根据这个函数计算出来的值与文献给出的值十分吻合。     

Upper crustal seismic structure of the Mazury complex and Mazowsze massif within East European Craton in NE Poland

The large-scale seismic experiment POLONAISE '97 (POlish Lithospheric ONsets—An International Seismic Experiment) was carried out in May 1997 in Poland, Lithuania, and Germany. Its main purpose was to investigate the structure of the crust and the uppermost mantle in the region of the Trans European Suture Zone (TESZ) that lies between the East European Craton (EEC) and the Palaeozoic Platform. This paper covers the interpretation of seismic data along the NW–SE-trending, 180-km-long profile P5 located on the EEC. The recordings were of a high quality with seismic energy clearly visible along the whole profile. We have not found waves refracted below the upper crust in first arrivals. In the NW part of the profile, we have delineated a high-velocity body with the P-wave velocity in the range of 6.5–6.75 km/s in the upper crust. It corresponds to the K trzyn anorthosite massif within the Mazury complex. The Mazowsze massif is rather uniformly characterized by P-wave velocities 5.9–6.05 and 6.2–6.35 km/s in two layers, respectively. Sufficient S-wave data were available to estimate the V_p/V_s ratio (as well as the Poisson ratio), being 1.80 (0.277) in the high-velocity body and 1.67 (0.220) in the upper crust.Apart from the 2-D model along the profile, results of 3-D modelling in the area of the P5 profile are presented. Using off-line recordings, we got P-wave velocity field up to 8 km/s below the P5 profile at the depth of about 40 km as well as horizontal extent of the high-velocity body.     

Dynamic Systems of Petroleum Accumulation in the Nanpu Depression,Bohai Bay Basin,China

<正>It is significant to distinguish the dynamic systems of petroleum accumulation(DSPA) for the understanding of petroleum accumulation and distribution.According to the formation pressure framework,genetic types of petroleum and characteristics of conduit systems,three dynamic systems of petroleum accumulation were identified in the vertical profile in the Nanpu depression,Bohai Bay basin.The deeper DSPA(including formations Es_3 to Es_2) is a sealed system with high-overpressure and high-mature self-sourced oil.Most of the crude oil in the system accumulated in the periods of late Oligocene(23.5 Ma) and late Pliocene(2.4 Ma).The middle DSPA(including formations Es_1 to Ed_1) is an overpressured half-sealed system with mature or lower-mature self-sourced oil.The accumulation of oil in the system also occurred in the late Oligocene(23.5 Ma) and late Pliocene(2.4 Ma).The shallower DSPA(including formations Ed_2 to Q) is a hydrostatic system with lower-mature aliensourced oil from the middle system.Oil within this system accumulated only in the late Pliocene period.The oil in the shallower system migrated vertically along the faults from the formerly accumulated oil in the middle system by lateral migration along the sandbodies,whereas petroleum accumulation in the deeper system was mainly derived from the system itself by lateral migration along the sandbodies and rarely migrated out of the system.In this case,it seems that the deeper system is a more potential exploration prospect in addition to the other two proved favorable systems.     

Process Analysis of In-situ Strain during the Ms8.0 Wenchuan Earthquake—Data from the Stress Monitoring Station at Shandan

Abstract: There were huge life and property losses during the Ms8.0 Wenchuan earthquake on May 12, 2008. Strain fluctuation curves were completely recorded at stress observatory stations in the Qinghai-Tibet plateau and its surroundings in the process of the earthquake. This paper introduces the geological background of the Wenchuan earthquake and the profile of in-situ stress monitoring stations. In particular, data of 174 earthquakes (Ms4.0-Ms8.5) were processed and analyzed with various methods, which were recorded at the Shandan station from August 2007 to December 2008. The results were compared with other seismic data, and further analyses were done for the recoded strain seismic waves, co-seismic strain stepovers, pre-earthquake strain valleys, Earth’s free oscillations before and after the earthquake and their physical implications. During the Wenchuan earthquake, the strainmeter recorded a huge extensional strain of 70 seconds, which shows that the Wenchuan earthquake is a rupture process predominated by thrusting. Significant precursory strain anomalies were detected 48 hours, 30 hours, 8 hours and 37 minutes before the earthquake. The anomalies are very high and their forms are very similar to that of the main shock. Similar anomalies can also be found in strain curves of other shocks greater than Ms7.0, indicating that such anomalies are prevalent before a great earthquake. In this paper, it is shown that medium aftershocks (Ms5.5-6.0) can also cause Earth’s free oscillations. Study of free oscillations is of great significance to understand the internal structure of the Earth and focal mechanisms of earthquakes and to recognize slow shocks, thus providing a scientific basis for the prevention and treatment of geological disasters and the prediction of future earthquakes.     

Late-Quaternary Slip Rate and Seismic Activity of the Xianshuihe Fault Zone in Southwest China

The Xianshuihe fault zone is a seismo-genetic fault zone of left-lateral slip in Southwest China. Since 1725, a total of 59 Ms ≥ 5.0 earthquakes have occurred along this fault zone, including 18 Ms 6.0–6.9 and eight Ms ≥ 7.0 earthquakes. The seismic risk of the Xianshuihe fault zone is a large and realistic threat to the western Sichuan economic corridor. Based on previous studies, we carried out field geological survey and remote sensing interpretation in the fault zone. In addition, geophysical surveys, trenching and age-dating were conducted in the key parts to better understand the geometry, spatial distribution and activity of the fault zone. We infer to divide the fault zone into two parts: the northwest part and the southeast part, with total eight segments. Their Late Quaternary slip rates vary in a range of 11.5 mm/a –(3±1) mm/a. The seismic activities of the Xianshuihe fault zone are frequent and strong, periodical, and reoccurred. Combining the spatial and temporal distribution of the historical earthquakes, the seismic hazard of the Xianshuihe fault zone has been predicted by using the relationship between magnitude and frequency of earthquakes caused by different fault segments. The prediction results show that the segment between Daofu and Qianning has a possibility of Ms ≥ 7.0 earthquakes, while the segment between Shimian and Luding is likely to have earthquakes of about Ms 7.0. It is suggested to establish a GPS or In SAR-based real-time monitoring network of surface displacement to cover the Xianshuihe fault zone, and an early warning system of earthquakes and post seismic geohazards to cover the major residential areas.     

Polarized X-ray absorption spectroscopy of metal ions in minerals

We have made use of the nearly complete linear polarization of synchrotron radiation to study the polarization dependence of X-ray absorption near-edge structure (XANES) and extended fine structure (EXAFS) in oriented single crystals of gillespite (BaFe²⁺ Si₄O₁₀; Fe^{2 +} in square-planar coordination, point symmetry C ₄), anatase (TiO₂; Ti⁴⁺ in octahedral coordination, point symmetry D _2d), and epidote (Ca₂(Al, Fe³⁺)₃SiO₄)₃(OH); Fe³⁺ in distorted octahedral coordination, point symmetry (C _s). For gillespite, the Fe K-XANES spectrum varies strongly with E-vector orientation of the incident X-ray beam. When the E-vector lies in the plane of the FeO₄ group (i.e., perpendicular to the c-axis), multiple-scattering features at 7127 and 7131 eV intensify, whereas when the E-vector is perpendicular to the plane of the FeO₄ group (i.e., parallel to the c-axis), a strongly-polarized 1s to 4p bound state transition occurs at 7116 eV and a localized continuum resonance occurs at 7122 eV. The Fe-K-EXAFS spectrum of gillespite is also highly polarization dependent. When the E-vector is perpendicular to c, all four nearest-neighbor oxygens around Fe²⁺ contribute to the EXAFS signal; when E is parallel to c, the EXAFS signal from nearest-neighbors is reduced by at least 86%. The unpolarized Ti K-XANES spectrum of anatase has three relatively strong pre-edge features at 4967.1, 4969.9, and 4972.7 eV which have resisted definitive interpretation in past studies. The lowest energy feature has a strong xy polarization dependence, suggesting a large amount of 4p _x,y character, and it is also very sharp, indicating a well-defined transition energy. Both of these observations are consistent with an excitonic state with a binding energy of 2.8 eV. The two higher energy features, which are characteristic of octahedrally-coordinated Ti⁴⁺, show little polarization dependence and are probably due to 1s to 3d bound-state transitions, with a small degree of np character in the final state wavefunction. Interpretation of the polarization dependence of Fe K-XANES spectra for epidote is not as straightforward due to the lower space group symmetry (P2₁/m) relative to gillespite (P4/ncc) and anatase (I4₁/amd) and the lower point group symmetry (C _s) of the M(3) site which contains most of the Fe³⁺ in the epidote structure. However, the presence of a shoulder at 7121 eV in the E parallel to b spectrum and its absence in the E normal to bc spectrum are consistent with it being a 1s to 4p _z bound-state transition. Strong, weakly x, y polarized features near 7126 eV in both spectra are most likely due to localized continuum transitions. Also, the 1s to 3d pre-edge intensity varies in intensity with E-vector orientation which is consistent with displacement of Fe³⁺ from the center of the M(3) octahedral site. Analysis of EXAFS spectra of epidote in these two polarizations yields bond distances which are within 0.04 Å of previous single-crystal X-ray diffraction analysis. This study demonstrates the utility of polarized X-ray absorption spectroscopy in quantifying the energies and orbital compositions of final state wavefunctions associated with various X-ray induced transitions in transition-metal containing minerals. It also shows that reasonably accurate M-O distances can be obtained for individual bonds oriented in crystallographically non-equivalent directions.     

Estimating tsunami potential of earthquakes in the Sumatra–Andaman region based on broadband seismograms in India

In this paper, we report that the ratio of broadband energy (0.01?C2?Hz) to high-frequency energy (0.3?C2?Hz), E _r, estimated from regional seismograms of India, might be a useful parameter in estimating tsunami potential of earthquakes in the Sumatra?CAndaman region. E _r is expected to be sensitive to the depth as well as to the source characteristics of an earthquake. Since a shallow and slow earthquake has a greater tsunamigenic potential, E _r may be a useful diagnostic parameter. We base our analysis on broadband seismograms of the great earthquakes of Sumatra?CAndaman (2004, M _w?~?9.2) and Nias (2005, M _w 8.6), 41 of their aftershocks, and the earthquakes of north Sumatra (2010, M _w 7.8) and Nicobar (2010, M _w 7.4) recorded at VISK, a station located on the east coast of India. In the analysis, we also included the two recent, great strike-slip earthquakes of north Sumatra (2012, M _w 8.6, 8.2) recorded at VISK and three south Sumatra earthquakes (2007, M _w 8.5; 2007, M _w 7.9; 2010, M _w 7.8) recorded at PALK, a station in Sri Lanka. We find that E _r is a function of depth; shallower earthquakes have higher E _r values than the deeper ones. Thus, E _r may be indicative of tsunamigenic potential of an earthquake. As M _w and E _r increase so does the tsunami potential. In addition to the parameter E _r, the radiated seismic energy, E _s, may be estimated from the regional seismograms in India using empirical Green??s function technique. The technique yields reliable E _s for the great Sumatra and Nias earthquakes. E _r and E _s computed from VISK data, along with M _w and focal mechanism, may be useful in estimating tsunami potential along the east coast of India from earthquakes in the Sumatra?CAndaman region in less than ~20?min.     

Extension of global creepex definition (Ms-mb) to local studies (Md-ML): the case of the Italian region

Creepex, defined as the systematic deviation of the magnitude of a single earthquake from the linear orthogonal regression between local magnitude M_L and coda duration magnitude M_d, calculated for the whole region, is used as a measure of the frequency content of the seismic sources in the Italian region. Predominantly high-frequency events are found in the two areas of Quaternary tectonic shortening in North-Central Italy and in the Calabrian Arc. This result, confirmed by two independent statistical tests, is in agreement with the global pattern obtained from the study of the regression between body-wave magnitude, m_b, and surface-wave magnitude M_s: systematic shift to high frequencies in the energy release of seismic sources located in subduction zones and to low frequencies in zones of spreading. The analysis of the correlation between the patterns of heat flow and of seismic source spectral properties indicates that these source properties, in general, do not reflect thermal conditions in the lithosphere, but rather represent the result of tectonic processes.     

Source parameters of the Burma-India border earthquake of July 29, 1970, from body waves

The source parameters are determined for the Burma-India border earthquake of July 29, 1970, from body-wave spectra. We obtain seismic moment [ , ] · 10²⁶ dyne cm, source dimension [ ] km, radiated energy [ , −E_R (S) = 1.35] · 10²⁰ ergs and the stress drop = 11 bars.     

国内外地震滑坡研究：现状、问题与展望

地震滑坡一直是地学界关注的研究热点之一,特别是2005年10月8日克什米尔M_s7.8地震、2008年5月12日汶川M_s8.0地震、2010年玉树M_s7.1地震、2014年云南鲁甸M_s6.5地震、2015年4月25日尼泊尔M_s8.1地震以及2017年8月8日九寨沟M_s7.0地震等一系列强震活动,都不同程度地触发了大量地震滑坡,并伴生了许多次生灾害,从而促进地震滑坡研究进入了新的阶段。本文通过总结国内外地震滑坡研究现状,归纳了地震滑坡的特征和分布规律、地震动参数与地震滑坡关系、边坡地质条件与地震滑坡关系、地震滑坡的动力响应特征、地震滑坡预测以及地震滑坡危险性与风险区划等6方面的研究现状,在此基础上进一步提出了未来地震滑坡应关注的主要研究方向和重点。

     

围压下裂隙对地震波速度影响的物理实验研究

构建了等直径不同裂隙密度和等裂隙密度不同裂隙直径两组物理模型,进行不同围压条件下多方向的超声波速度测试,并运用Hudson理论进行了理论模型计算。结果显示,计算与实测结果吻合较好。随围压的增大,纵、横波速度均近线性增加,纵、横波各向异性基本保持不变;裂隙密度从2%增大到6%,纵波速度不同程度降低,其中慢纵波降低幅度相对较大,快横波变化不明显,而慢横波则大幅降低。随着裂隙密度的增大,纵、横波各向异性均增大,且横波各向异性增加速率大于纵波;裂隙直径从2 mm增大到3 mm,快纵波速度增加很小,慢纵波增加明显,横波速度均不发生改变。随着裂隙直径的增大,纵波各向异性逐渐降低,横波各向异性保持不变。最后结合试验结果分析了Hudson理论在不同深度进行参数预测的必要条件。研究结果有助于油气生产、地下水的开采与控制、污染处理等。      

Static stress drop inferred from near-fault parameters for the 1999 Chi–Chi earthquake in Taiwan

In 2001, a special issue of the Bulletin of the Seismological Society of America (BSSA) featured seismological research for the 1999 Chi–Chi Taiwan earthquake. This study uses source parameters suggested by the first author in this special issue to estimate static stress drop associated with the Chi–Chi earthquake. The waveform simulation method was used to carefully examine these source parameters. The simulation results indicate that source parameters, inferred from near-fault observations, are well determined. According to the rupture area and slip, the static stress drops (Δσ_s) obtained were distributed between a small value of 47 bars near the epicentral region and a much larger value (>200 bars) to the north. Similar trends in dynamic stress drop (Δσ_d) were also recognized by the first author in his paper published in 2001 BSSA special issue. Comparing the Δσ_s with Δσ_d, satisfies the relation Δσ_s/Δσ_d ≈ 1. This relation suggests that fault motion is mostly spent releasing seismic wave energy during the rupture process of the Chi–Chi earthquake. The consistency between static and dynamic stress drops thus provides a measure of energy-moment (E_s/M₀) ratios, which range from 9.0 × 10⁻⁵ to 6.5 × 10⁻⁴. The average E_s/M₀ ratio estimated for the northern portions of the fault is 3.4 × 10⁻⁴, which is about 3 times that of the south. Such a high E_s/M₀ ratio can be interpreted as having low strength in the rupture for the northern portions of the fault, where the fault would release less energy per unit rupture surface to create the new rupture.     

东濮凹陷北部沙三段古压力恢复及油气运聚动力构成

综合应用流体包裹体法和盆地模拟法, 恢复了东濮凹陷北部沙三段古压力, 并分析了成藏期油气运聚动力构成。研究结果表明, 沙三段超压分布受构造格局、沉降中心、生烃中心控制明显, 超压幅度表现为洼陷区大、中央隆起带次之、西部斜坡带最小, 受盐岩层发育影响, 濮卫-文留地区盐岩下部层系表现为压力系数高值区。成藏期超压和浮力是研究区沙三段油气运聚的主要动力, 压力过渡带和正常压力带是油气的主要聚集场所。研究区主要存在超压驱动、超压-浮力联合驱动和浮力驱动等3种类型的驱动机制, 其中斜坡带和洼陷带等超压带主要为超压驱动, 部分中央隆起带上的压力过渡带为超压-浮力联合驱动, 西部斜坡带和部分中央隆起带等正常压力带主要为浮力驱动。      

Empirical modeling of plate load test moduli of soil via gene expression programming

New empirical models were developed to predict the soil deformation moduli using gene expression programming (GEP). The principal soil deformation parameters formulated were secant (E_s) and reloading (E_r) moduli. The proposed models relate E_s and E_r obtained from plate load-settlement curves to the basic soil physical properties. The best GEP models were selected after developing and controlling several models with different combinations of the influencing parameters. The experimental database used for developing the models was established upon a series of plate load tests conducted on different soil types at depths of 1–24 m. To verify the applicability of the derived models, they were employed to estimate the soil moduli of a part of test results that were not included in the analysis. The external validation of the models was further verified using several statistical criteria recommended by researchers. A sensitivity analysis was carried out to determine the contributions of the parameters affecting E_s and E_r. The proposed models give precise estimates of the soil deformation moduli. The E_s prediction model provides considerably better results in comparison with the model developed for E_r. The simplified formulation for E_s significantly outperforms the empirical equations found in the literature. The derived models can reliably be employed for pre-design purposes.