Increased correlation range of seismicity before large events manifested by earthquake chains

“Earthquake chains” are clusters of moderate-size earthquakes which extend over large distances and are formed by statistically rare pairs of events that are close in space and time (“neighbors”). Earthquake chains are supposed to be precursors of large earthquakes with lead times of a few months. Here we substantiate this hypothesis by mass testing it using a random earthquake catalog. Also, we study stability under variation of parameters and some properties of the chains. We found two invariant parameters: they characterize the spatial and energy scales of earthquake correlation. Both parameters of the chains show good correlation with the magnitudes of the earthquakes they precede. Earthquake chains are known as the first stage of the earthquake prediction algorithm reverse tracing of precursors (RTP) now tested in forward prediction. A discussion of the complete RTP algorithm is outside the scope of this paper, but the results presented here are important to substantiate the RTP approach.     

Alternatives to tap water: a case study of the Gaza Strip, Palestine

Gaza Strip is a highly populated, small area in which the groundwater is the main water source. During the last few decades, groundwater quality has deteriorated to a limit that the municipal tap water became brackish and unsuitable for human consumption in most parts of the strip. To overcome this serious situation, several attempts were made to replace the tap water or to improve its quality by water desalination applying the reverse osmosis (RO) technology, bottled water, importing water, and storm water harvesting. Water desalination, which is widely applied in the strip, has environmental and health risks. Brine water resulting from the RO systems has adverse environmental effects, whereas the produced fresh water with very low chemical concentrations may have harmful effects on human health. Therefore, introducing international water quality standards to delineate the harmless minimum limits is necessary. Storm water harvesting can provide a partial solution for the water problems in the strip.     

Salt tectonics and tear faulting in the central part of the Zagros Fold-Thrust Belt,Iran

The central part of the Zagros Fold-Thrust Belt is characterized by a series of right-lateral and left-lateral transverse tear fault systems, some of them being ornamented by salt diapirs of the Late Precambrian–Early Cambrian Hormuz evaporitic series. Many deep-seated extensional faults, mainly along N–S and few along NW–SE and NE–SW, were formed or reactivated during the Late Precambrian–Early Cambrian and generated horsts and grabens. The extensional faults controlled deposition, distribution and thickness of the Hormuz series. Salt walls and diapirs initiated by the Early Paleozoic especially along the extensional faults. Long-term halokinesis gave rise to thin sedimentary cover above the salt diapirs and aggregated considerable volume of salt into the salt stocks. They created weak zones in the sedimentary cover, located approximately above the former and inactive deep-seated extensional faults. The N–S to NNE–SSW direction of tectonic shortening during the Neogene Zagros folding was sub-parallel with the strikes of the salt walls and rows of diapirs. Variations in thickness of the Hormuz series prepared differences in the basal friction on both sides of the Precambrian–Cambrian extensional faults, which facilitated the Zagros deformation front to advance faster wherever the salt layer was thicker. Consequently, a series of tear fault systems developed along the rows of salt diapirs approximately above the Precambrian–Cambrian extensional faults. Therefore, the present surface expressions of the tear fault systems developed within the sedimentary cover during the Zagros orogeny. Although the direction of the Zagros shortening could also potentially reactivate the basement faults as strike-slip structures, subsurface data and majority of the moderate-large earthquakes do not support basement involvement. This suggests that the tear fault systems are detached on top of the Hormuz series from the deep-seated Precambrian–Cambrian extensional faults in the basement.     

Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.     

基于起伏地形平化策略的弹性波逆时偏移成像方法

随着能源和资源勘查开采工作的深入,地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一,能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓,由于传统的地震波场延拓技术往往基于水平地表条件,相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法：基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法,采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造,获得了相应起伏地形条件下的复杂几何模型,开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性,可适应不同类型起伏地表采集的地震资料,显示出该方法在地震勘探领域的良好应用前景.     

人工源宽角地震实验的谱元法模拟和逆时偏移成像方法研究

人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。     

Multipathing in three‐parameter common‐image gathers from reverse‐time migration

Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region.     

矿液致裂与川口三角潭大脉型钨矿床的形成

川口三角潭钨矿床与华南众多同类型矿床有别 ,矿床不具“五层楼”特征 ,而具逆向分带现象 ,是一个受岩浆岩和断裂构造双重控制 ,矿液致裂作用形成的高—中温 ,以黑钨矿为主的钨—多金属矿床。     

红砂和珍珠对极端环境的生理响应与调节机制

红砂和珍珠的生理生态特性随其生长的环境条件而变化;长期生长在一起的红砂和珍珠其生理生态特性具明显的差异,两者适应逆境的机制亦不相同。红砂和珍珠作为超旱生植物,其生理生态学特性明显区别于其他沙生植物;生长在山前地带的红砂和珍珠,可以在极低的水势条件下保证其生命活动,低的水势可能与渗透调节物质（如大量的脯氨酸）的存在有关,这类物质增加了植物的吸水能力;红砂和珍珠的光抑制程度较其他沙生植物要高得多,且光抑制的产生已对两种灌木光合色素产生了破坏,然而在条件稍好的早晚光抑制可得到恢复。     

基于保辛算法的声波叠前逆时偏移

叠前逆时偏移是目前成像精度最高的地震偏移方法之一,其实现过程中的一个重要步骤是数值求解全波方程,所以快速有效求解全波方程的数值算法对逆时偏移至关重要. 四阶近似解析辛可分Runge-Kutta (NSPRK) 方法是近年发展的一种具有高效率、高精度的数值求解波动方程的保辛差分方法, 能在粗网格条件下有效压制数值频散, 从而提高计算效率, 节省计算机内存需求量. 本文利用四阶NSPRK方法构造的基本思想,发展了具有六阶空间精度的NSPRK方法,并对新的六阶NSPRK方法进行了详细的稳定性和数值频散分析,以及计算效率比较和波场模拟. 同时将该方法用于声波叠前逆时偏移中, 得到一种时间上保辛、空间具有六阶精度、低数值频散、可应用大步长进行波场延拓并能长时计算的叠前逆时偏移方法,对Sigsbee2B模型进行了偏移成像, 并和四阶NSPRK方法、传统的六阶差分方法、四阶Lax-Wendroff correction (LWC) 方法进行了对比. 数值结果表明, 基于六阶NSPRK方法的叠前逆时偏移能得到更好的成像结果, 是一种优于四阶NSPRK方法、传统的六阶差分方法、四阶LWC叠前逆时偏移的方法, 尤其是在粗网格情况下具有更明显的优越性.