离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.     

Crustal structure of the Ruby Mountains and southern Carlin Trend region, Nevada, from magnetotelluric data

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings reside in three east–west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass, and Secret Pass latitudes). One short cross-line was also taken to assess an east–west structure to the north of the northern profile. Model resistivity cross-sections were derived from the MT data using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity–thickness product) from east to west. These increases are attributed to graphitization caused by Elko–Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, steeply dipping fault zones. The zones provide pathways to the lower crust for heterogeneous, upper crustal induced, electric current flow. Resistive core complex crust appears steeply bounded under the middle of the neighboring grabens and not to deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness (T_e) estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east–west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to nearby graphitic metasediments, with possible alteration by middle Tertiary magmatism. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (<1 vol.%) containing hypersaline brines and possible water-undersaturated crustal melts, residual to the mostly Miocene regional extension. The resistivity expression of the southern Carlin Trend (CT) in the Pinon Range is not a simple lineament but rather a family of structures attributed to Eocene intrusion, stratal deformation, and alteration/graphitization. Substantial reactivation or overprinting by core complex uplift or Basin–Range extensional events seems likely. We concur with others that the Carlin Trend may result in part from overlap of the large Eocene Northeast Nevada Volcanic Field with Precambrian–Paleozoic deep-water clastic source rocks thickening abruptly to the west of the Pinon Range, and projecting to the north–northwest.     

2000年4月29日内乡ML4.7级地震前潢川地电异常变化

2000年4月29日11时54分，河南省南阳市内乡县马山口镇与镇平县寺山乡之间发生一次ML4.7级地震。地震前3天，潢山台视电阻率ρs异常较为明显。分析了此次异常与内乡ML4.7级地震的关系，认为此次异常很有可能是内乡ML4.7级地震的短临前兆异常。当观测台站与震源位于同一构造带时，即便离震中较远，地电阻率对地震异常反映仍较灵敏，且与构造带同一走向的测向尤为明显。     

Mechanism of the June 4, 2000 Southern Sumatra, Indonesia, earthquake

Introduction The June 4, 2000 southern Sumatra, Indonesia, earthquake (16h28min26.2s UTC, 4.72, 102.09, 33 km, MS=8.0 [IRIS]) occurred under the Indian Ocean, near the Mentawai fault, along the well-known Sumatran subduction zone and the great Sumatran fault, all of which trend northwest-southeast. 1 800 houses were totally destroyed, 10 196 were heavily damaged, and 18 378 were slightly damaged by the earthquake. At least 97 people were killed, 1 900 were injured, and 122 000 were left…     

低电阻率油气层物理参数变化机理研究

研究了淡水环境下富含粘土矿物、细颗粒砂岩低电阻率油气层. 分析了粘土的附加导电性,认为粘土质量百分数高且富含伊利石和蒙脱石的地层中,若地层水质量浓度小于15g/L等效NaCl盐溶液时,粘土的附加导电性是造成油气层电阻率低的首要原因. 同时分析了细颗粒、粘土质量百分数不高的储层,高的毛细管压力使束缚水饱和度很高,同样能形成低电阻率油气层. 储层中粘土质量百分数的多少对储层岩石电性特征(I-Sw曲线形态)有重要影响,可用于指导测井资料解释模型的选择.     

Moment-tensor statistics of the 1997 Long Valley microearthquake swarm

Abstract. A simple, fast, moment-tensor inversion method using bandpass-filtered P-amplitudes was used to study the moment-tensor statistics of Long Valley caldera microearthquakes. The events were recorded in the summer of 1997, during a swarm in the caldera. The swarm was associated with geodetic extension, uplift, and subsequent moderate earthquake activity. The moment tensor solutions for 1,993 events were calculated using the new method. The majority of the resulting focal mechanisms appear to be explained in terms of double couple mechanisms. Since some events did exhibit considerable deviation from double-couples, the moment data were studied for their statistical significance. The moments of the actual data were compared to the moments of synthetic data with varying degrees of random noise in their spectra. The results of this study suggested that unless data from more than 20 stations are used and the earthquake epicenter is located inside or very close to the network area, moment-tensor inversion does not correctly resolve the non-double-couple components of microearthquakes. Analysis of the inversion residuals shows that the average noise in the P-wave spectra was close to 20%. The fluctuations of the volumetric components of the moment-tensor are in good agreement with those of the synthetic pure double-couples with 20% of added noise. Thus the moment-tensor statistics suggests that little if any volume change is required to explain the observed seismic energy release in the swarm. However, the statistics do show that a significant compensated-linear-vector-dipole component maybe present in the bulk of the seismicity. Given the network used in this study, such a component could not be precisely resolved for individual earthquakes. This possibility deserves further investigation because of its bearing on the nature of fluid-fault-earthquake processes in swarms.     

山西地电中短期前兆异常识别方法浅析

用傅氏滑动法、动态距平法、速率法、均值法等对山西大同台、代县台、临汾台的地电异常进行了分析 ,对地电阻率异常的判定、预报方法效能进行了评价 ,对不同的台站 ,从正常的变化形态中找出了突出异常的有效方法 ,从中提取出因地壳应力变化引起的地震前兆异常 ,结果表明 ,5 .0级以上地震在Δ=15 0 km范围内有明显的中短期前兆异常     

Relation of the Wave-Propagation Metric Tensor to the Curvatures of the Slowness and Ray-Velocity Surfaces

The contravariant components of the wave-propagation metric tensor equal half the second-order partial derivatives of the selected eigenvalue of the Christoffel matrix with respect to the slowness-vector components. The relations of the wave-propagation metric tensor to the curvature matrix and Gaussian curvature of the slowness surface and to the curvature matrix and Gaussian curvature of the ray-velocity surface are demonstrated with the help of ray-centred coordinates.