层状方位各向异性介质的视电阻率计算

从电性各向异性的欧姆定律出发，推导了直流电法层状方位各向异性介质中的电位分布、边界条件及视电阻率计算公式．以四极对称装置系统为例，对具有相同各向异性系数的4层模型采用核函数递推法作了理论数值模拟，得到了不同方向的电阻率测深曲线及其等值线形态．结果表明理论公式是正确的，测深曲线既反映了分层介质的电阻率差异，又反映了各层中电阻率的各向异性特征．      

层状介质任意形状回线源瞬变电磁全区视电阻率的研究

利用时域垂直磁场求取全区视电阻率.任意形状回线的垂直磁场可以看成是多个水平电偶极子的垂直磁场累加而成;因此,不论回线是何形状,只要将其分解为水平电偶极子,就可利用水平电偶极子的垂直磁场求取任意形状回线的全区视电阻率.     

电阻率各向异性的研究

概述了国内外有关地电阻率各向异性实验方面的研究现状，总结了迄今所提出的地电阻率各向异性产生的机理；提出存在的问题和以后的发展方向等。     

Crustal anisotropy of Taihang Mountain Range using azimuthal variation of receiver functions

We discussed the possibility of studying crust anisotropy by analyzing azimuthal variation of the receiver functions and presented a technique for computing the transmission response of a flat-layered medium with arbitrarily ori-ented hexagonally symmetric anisotropy using the reflectivity algorithm. Using this method we investigated the crust anisotropy of Taihang Mountain Range (TMR). Our result shows that there is significant anisotropy with a slow symmetry axis in the upper crust and a fast symmetry axis in the lower crust. The anisotropy magnitude of about 8%~15% is found in the upper crust and a smaller magnitude of about 3%~5% in the lower crust. Orienta-tion of the symmetry axes and the depth of anisotropy appearance as deduced from the seismic records of four individual seismic stations are different from each other. The crust anisotropy beneath the four stations may be associated with the local crustal fabrics in a small area.     

地震地电阻率各向异性变化数值模拟

将APE模式中差应力与裂隙的演化理论与前期研究中建立的张量形式的各向异性电性模型相结合,推导出了裂隙演化参数(密度、孔压、纵横比、闭合角)与电导率关系的动态方程,用数值模拟方法研究了饱水岩石在水平向应力作用下电阻率的各向异性变化．并结合实验和观钡5结果,讨论了电阻率各向异性变化的可能机理。     

转换波方位各向异性裂缝检测技术研究及应用

HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果.     

利用接收函数方位变化研究太行山区地壳各向异性

探讨了利用接收函数方位变化研究各向异性复杂地区速度结构的可能性. 详细介绍了具有任意对称轴方向的横向各向同性系统中介质响应的计算过程,并将利用接收函数方位变化研究地壳各向异性的方法应用于华北太行山地区. 结果表明,研究台站下方地壳具有明显的各向异性特征,其中地表盖层表现为8％——15％的慢轴型各向异性;深部地壳各向异性略小,表现为3％——5％的快轴型各向异性特征. 不同台站获得的各向异性对称轴的方向及深度分布不具有明显的一致性. 区域构造环境的变化可能是导致地壳各向异性存在的原因.     

Low-frequency anisotropy in fractured and layered media

Computing effective medium properties is very important when upscaling data measured at small scale. In the presence of stratigraphic layering, seismic velocities and anisotropy parameters are scale and frequency dependent. For a porous layer permeated by aligned fractures, wave-induced fluid flow between pores and fractures can also cause significant dispersion in velocities and anisotropy parameters. In this study, we compare the dispersion of anisotropy parameters due to fracturing and layering at low frequencies. We consider a two-layer model consisting of an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce anisotropy parameters dispersion due to fractures (meso-scale) in the sand layer. This intrinsic dispersion is added to anisotropy parameters dispersion induced by layering (macro-scale) at low frequencies. We derive the series coefficients that control the behaviour of anisotropy parameters at low frequencies. We investigate the influences of fracture length and fracture density on fracturing effect, layering effect and combined effect versus frequency and volume fraction of sand layer. Numerical modelling results indicate that the frequency dependence due to layering is not always the dominant effect of the effective properties of the medium. The intrinsic dispersion is not negligible compared with the layering effect while evaluating the frequency-dependent properties of the layered medium.     

The variations of the degree of ground resistivity anisotropy during the Tangshan earthquake

ＴｈｅｖａｒｉａｔｉｏｎｓｏｆｔｈｅｄｅｇｒｅｅｏｆｇｒｏｕｎｄｒｅｓｉｓｔｉｖｉｔｙａｎｉｓｏｔｒｏｐｙｄｕｒｉｎｇｔｈｅＴａｎｇｓｈａｎｅａｒｔｈｑｕａｋｅＴｏｎｇ－ＥｎＭＡＯ（毛桐恩）；Ｔｉｅ－ＣｈｅｎｇＷＡＮＧ（王铁成）；Ｊｉａ－ＬｉｕＹＡＯ...     

横波分裂和纵波方位各向异性在油气勘探中的应用

本文首先回顾了横波分裂和纵波方位各向异性应用于油气勘探的发展历程; 在此基础上重点讨论了二者在上地壳中的分布情况、相互关系及应用条件.横波分裂是地震波传播经过裂缝时所独有的横波双折射现象, 40年的科学研究及生产应用表明这一物理现象在上地壳中广泛存在, 并主要存在于近地表1200 m以浅.纵波方位各向异性从浅到深在上地壳中分布差异较大; 应用于检测裂缝时仅对含气裂缝发育区敏感, 对含其他流体的裂缝区不敏感.针对裂缝发育区的预测, 本文利用理论和实验数据分析论证了横波分裂与纵波方位各向异性的关系; 阐明在油气勘探中横波分裂分析与纵波方位各向异性预测裂缝储层须首先消除近地表的影响, 这在一定程度上限制了二者的应用.同时阐明二者的应用还需要满足近地表相对简单、裂缝储层相对较厚等条件; 且需采集非零偏或环形垂直地震剖面(VSP)资料进行标定.即多分量地震数据品质、近地表地震地质条件及储层方位各向异性发育程度是影响横波分裂及纵波方位各向异性应用效果的3个关键因素.

     

视电阻率及其差分曲线在岩土分层中的应用

将视电阻率曲线与视电阻率一阶差分曲线应用于岩土工程施工中．对两种曲线比较后发现,两种曲线均为有效的视电阻率信息．视电阻率一阶差分曲线具有较视电阻率曲线更高的分辨率．两种方法结合使用能较好地对岩土进行分层．当视电阻率存在明显差异且岩层厚度较大时,可用视电阻率曲线对岩土进行分层;当视电阻率差异较小时,则用视电阻率一阶差分曲线对岩土进行分层．     

Research on dependence of resistivity changing anisotropy on microcracks extending in rock with experiment

Intr0ducti0nTheproblemsaboutcrustalrocksfracturingandfricti0nalslippingalongthefracturesurfaCearecloselyrelatedtotheresearchesofearthquakemechanismandgeol0gichazards,andalsoareoneofthosebeinggivenmoreconcernbyscientistsinearthscience-Discussingab0utthegenera-tionanddevelopmentprocess0fintraPlateearthquakes,Wang,etal(l999)proPOsedthatearth-quake'spregnant,mechanicspropertyandgeometricstructureofmatterinearthquake'sgestatingareaaredifferinthousandswaysf0rdifferentplaceinsideslab-Thus,thepregna…     

Dynamic evolution patterns of the degree of ground resistivity anisotropy and the seismogenic process

Following a new train of thinking, this paper has explored first the potential information in the ground resistivity data observed by the existing geoelectric observation system, investigated and proposed a new dimensionless geoelectric precursor factor, the degree of ground resistivity anisotropy, S, and studied the characteristics of dynamic evolution pattern of S during the seismogenic process. The results show that, during the seismogenic process, the degree of ground resistivity anisotropy (S) displays a process of 'normal' →'abnormal strengthening (amplitude, range)' →abnormal weakening' → 'earthquake occurrence'→ 'normal'. The earthquake would occur at the time when the S value has entered the late stage of strengthening and turns to weaken and in the gradient belt on the margin of S anomaly region. The dynamic evolution pattern of S reflects the changes of the tectonic stress field during the seismogenic process. Therefore, it would be possible to trace the process of earthquake generation and occurrence from the dynamic evolution pattern of S so as to service earthquake prediction.     

Detecting rupture precursors and determining the main fracture spread direction of rock with dynamic rock resistivity change anisotropy

In the past, domestic and abroad scientific workers have done a large number of experiment researches on a great number of researches of experiment for the relationship between electrical resistivity change and the load stress on rock (Geoelectric Testing Group, Department of Geophysics, Peking University and Research Division, Seismological Brigade of Lanzhou, 1978; Zhang, 1981; Chen, et al, 1983; Zhao, et al, 1983; Zhang, Lu, 1983; Zhang, et al, 1985; Brace, Orange, 1968; Kurite, 1986; Teisseyre, 1989), such as imitating dynamic fracture of rock, frictional sliding course, and observing its resistivity change shape with a lot of experimental workssuch as simulating the dynamic course of crustal rocks fracturing and frictional slipping along the fracture surface, and observing its resistivity change shape with experiment. They found that, in the fracture process of rock under loading, not only its resistivity shape change notably, and its direction changes also very clearly, and therefore many scientific workers have shown strong interest in anisotropy of resistivity change, and have done some researches onfor it.     

Traveltime approximation in vertical transversely isotropic layered media

The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging.     

欧亚大陆及西太平洋边缘海地区瑞雷波群速度结构与各向异性分布

根据欧亚大陆及西太平洋地区102个数字化台站记录的近万次地震事件的长周期垂直向瑞雷波资料,利用时频分析方法测量并筛选后共得到11213条质量较高的基阶瑞雷波群速度频散资料.纯路径频散的反演中同时计算方位各向异性,反演获得了欧亚大陆及西太平洋边缘海地区(10°E—150°E,10°S—80°N)8—200s共28个周期的瑞雷波群速度及各向异性空间分布图象.瑞雷波高分辨率层析成像表明,30—60s周期,以青藏高原为中心呈极低速分布;100—120s周期,速度差异幅度较大,在东亚东部及西太平洋边缘海,自北向南显示出一条宽2500—4000km,长约8000km的巨型低速异常带.相对海洋来说,欧亚大陆各向异性强度较弱且快波方向较复杂.由于受到印度板块与欧亚板块的碰撞,中国大陆西部的各向异性强度明显大于东部.     

Dynamic evolution patterns of the degree of ground resistivity anisotropy and the sets mogenic process

IntroductionTothefinalanalysis,earthquakepredictionisstillamarginalscienceofextremedifficultyintheworld.Tomakesomebreakthroughorpromoteearthquakeprediction,itisnecessarytostrengthentheresearchonnewtrainofthinking,newtheoryandnewmethods.Theresearchontheanisotropyofcrustalmediaduringtheseismogenicprocesshasbeentakenmoreandmoreseriouslybyseismologists.Inthepast,manyscholarsoftenmadestudiesonthefeaturesofmediainresponsetoseismicwavesandachievedsomeprogress(Teng,etal,1992).Thispaperstudiestheelect…     

three-dimensional sensitivity coefficients of apparent resistivity and preliminary application

We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.     

Anisotropy of the apparent resistivity variation rate near the epicentral region for strong earthquakes

IntroductionThe geo-electrical anisotropy before earthquake is of special significance for the prediction of earthquake location and for the research on crustal stress status near the epicentral region. The anisotropy has been studied from various aspects (CHEN, et al, 1983; QIAN, et al, 1985, 1996; LU, et al, 1990a; MAO, et al, 1995; RUAN, et al, 1999). But the anisotropic phenomena presented by geo-electrical data are not practically satisfactory and not clear as compared with the res…     

A method for inversion of layered shear wavespeed azimuthal anisotropy from Rayleigh wave dispersion using the Neighborhood Algorithm

Seismic anisotropy provides important constraints on deformation patterns of Earth's material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm(NA) for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depthdependent VSV(or the elastic parameter L) as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polarization azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.