Optimized arrays for 2D cross‐borehole electrical tomography surveys

The use of optimized arrays generated using the ‘Compare R’ method for cross‐borehole resistivity measurements is examined in this paper. We compare the performances of two array optimization algorithms, one that maximizes the model resolution and another that minimizes the point spread value. Although both algorithms give similar results, the model resolution maximization algorithm is several times faster. A study of the point spread function plots for a cross‐borehole survey shows that the model resolution within the central zone surrounded by the borehole electrodes is much higher than near the bottom end of the boreholes. Tests with synthetic and experimental data show that the optimized arrays generated by the ‘Compare R’ method have significantly better resolution than a ‘standard’ measurement sequence used in previous surveys. The resolution of the optimized arrays is less if arrays with both current (or both potential) electrodes in the same borehole are excluded. However, they are still better than the ‘standard’ arrays.     

唐山黄庄岩溶塌陷形成条件和机理分析

为查明唐山黄庄岩溶塌陷形成地质条件与机理,文章从该区岩溶塌陷的地质、地层岩性、水文地质特征对岩溶塌陷的影响等方面入手,研究岩溶地面塌陷的成因诱发条件,利用综合物探方法、钻探及地下水长期观测等技术手段,采用二维高密度测网数据,进行三维高密度数据处理和三维结果多角度分析,研究岩溶塌陷区地下三维空间分布形态、水流运移通道和基岩分布等。研究区岩溶地面塌陷应属于渗透潜蚀作用,塌陷坑分布于地下水的降落漏斗范围内和流水通道区域,特别集中于地下水主要径流方向和水交替强烈的地下水联系密切地段,雨季水气压力变化极为频繁、剧烈,加剧了土洞形成。      

Non‐uniqueness with refraction inversion – the Mt Bulga shear zone

The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable.     

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

Dynamic interaction of twin vertically overlapping lined tunnels in an elastic half space subjected to incident plane waves

The scattering of plane harmonic P and SV waves by a pair of vertically overlapping lined tunnels buried in an elastic half space is solved using a semi-analytic indirect boundary integration equation method. Then the effect of the distance between the two tunnels, the stiffness and density of the lining material, and the incident frequency on the seismic response of the tunnels is investigated. Numerical results demonstrate that the dynamic interaction between the twin tunnels cannot be ignored and the lower tunnel has a significant shielding effect on the upper tunnel for high-frequency incident waves, resulting in great decrease of the dynamic hoop stress in the upper tunnel; for the low-frequency incident waves, in contrast, the lower tunnel can lead to amplification effect on the upper tunnel. It also reveals that the frequency-spectrum characteristics of dynamic stress of the lower tunnel are significantly different from those of the upper tunnel. In addition, for incident P waves in low-frequency region, the soft lining tunnels have significant amplification effect on the surface displacement amplitude, which is slightly larger than that of the corresponding single tunnel.     

青藏高原东南缘及邻区地壳介质散射强度变化特征

基于远震扰动场方法,利用47个数字地震观测台记录的10次远震记录对青藏高原东南缘及邻区地壳介质散射强度进行了研究.结果表明研究区内地壳介质横向非均匀性强烈,扬子块体地壳介质散射强度与青藏块体存在较为明显的差异;地壳介质散射强度梯度带与主要断裂带分布趋势一致,龙门山断裂带、安宁河断裂带与地壳介质散射强度梯度带具有较强的相关性,散射强度高值区偏向青藏块体一侧;强震震中沿散射强度梯度带分布,其震中偏向散射强度高值区一侧;散射强度呈现纵向非均匀性特征,下地壳介质散射强度高于上地壳.研究区强、弱散射强度的空间分布可能与区内强烈的地质构造运动、频繁的地震及火山活动所引起深部物质破碎、熔融及其沿断裂带的上涌与运移有关.     

Laboratory measurements of guided‐wave propagation within a fluid‐saturated fracture

A fluid‐saturated flat channel between solids, such as a fracture, is known to support guided waves—sometimes called Krauklis waves. At low frequencies, Krauklis waves can have very low velocity and large attenuation and are very dispersive. Because they propagate primarily within the fluid channel formed by a fracture, Krauklis waves can potentially be used for geological fracture characterization in the field. Using an analogue fracture consisting of a pair of flat slender plates with a mediating fluid layer—a trilayer model—we conducted laboratory measurements of the velocity and attenuation of Krauklis waves. Unlike previous experiments using ultrasonic waves, these experiments used frequencies well below 1 kHz, resulting in extremely low velocity and large attenuation of the waves. The mechanical compliance of the fracture was varied by modifying the stiffness of the fluid seal of the physical fracture model, and proppant (fracture‐filling high‐permeability sand) was also introduced into the fracture to examine its impact on wave propagation. A theoretical frequency equation for the trilayer model was derived using the poroelastic linear‐slip interface model, and its solutions were compared to the experimental results.     

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

圆弧状凹陷地形对平面SV波散射问题的宽频带解析解

利用波函数的Fourier-Bessel级数展开法,推导了具有不同深宽比的圆弧状凹陷地形对入射平面SV波二维散射问题的解析解。区别于现有其他解析解,利用柱函数的渐近性质,使得散射波的待定系数得以直接确定,避免了线性方程组的求解以及相应的高频波入射下的数值计算问题,从而拓展了解析解适用的频带范围。通过与已有解析解的比较论证了该解析解的正确性,并在一个较宽的频带范围内研究了圆弧状凹陷地形对入射平面SV波的散射效应。