斜向逆冲断层相关褶皱的正演模型与实例分析

斜向逆冲作用在自然界普遍存在,研究斜向逆冲断层相关褶皱的构造几何学特征,识别断层相关褶皱是否存在斜向逆冲有重要意义。文章采用Trishear 4.5、Gocad以及Trishear3D软件构建一系列不同滑移量的断层转折褶皱和断层传播褶皱的二维正演剖面,通过连接一系列不同排列方式的二维剖面建立了三种不同逆冲滑移方向的断层转折褶皱和断层传播褶皱的假三维模型,通过不同假三维模型的比较分析来探讨斜向逆冲断层相关褶皱的构造几何学特征。研究发现,斜向逆冲断层相关褶皱区别于正向逆冲断层相关褶皱的特征主要有两点：① 正向逆冲断层相关褶皱层面等高线图上的最高点与后翼等高线中点的连线以及水平切面上的核心点与后翼中点的连线方向均与断层走向垂直,而斜向逆冲断层相关褶皱的最高点以及核心点与后翼中点的连线方向均与断层走向斜交,并且最高点与后翼等高线中点的连线方向或者核心点与后翼中点的连线方向均与逆冲滑移方向一致;② 在褶皱平行断层走向纵剖面上,正向逆冲断层相关褶皱各个层面最高点的连线是直立的,而斜向逆冲断层相关褶皱各个层面最高点的连线发生倾斜。通过这两个特征可以判别褶皱是否存在斜向逆冲以及逆冲的方向。将模型分析结果运用到四川盆地西南部三维地震勘探资料所覆盖的邛西背斜和大兴西背斜的实例中。研究结果表明,两个背斜均存在右旋斜向逆冲,逆冲方向与各自断层走向的夹角均为70°左右,邛西背斜和大兴西背斜的逆冲方向分别是NE79°和NE77°左右,这与龙门山南段晚上新世以来的主应力方向以及反演的汶川地震最大主应力方向一致。     

2008年汶川地震破裂的滑移向量

估计同震滑移向量对于认识和理解破裂方式和破裂过程具有重要意义。2008年汶川大地震在青藏高原东缘龙门山推覆构造带的中央断裂和前山断裂上各形成了一条长250 km和72 km的地表破裂带。地震发生后至今,已经发表了大量有关同震位错沿破裂带分布的论文和报告,但绝大部分都仅仅是破裂的走向位错和垂直位错,极少有同震滑移向量的报道。这不仅是因为野外难以直接测量到水平缩短量(或拉张量),而且还因为这些走滑位错实际上是视走滑位错,部分或全部来自水平缩短或拉张。因此,仅仅根据视走滑同震位错和垂直同震位错估计的同震总滑移量肯定包含了相当大的误差。尝试利用据不同走向参考线测量到的一组(两个以上)视走滑位错来计算水平滑移向量的这一新方法,获得了中央破裂带上的7个水平同震滑移向量,并结合垂直位错量进一步计算了走滑、倾滑和水平缩短三个同震滑移分量以及断层倾角和破裂面上的同震滑移向量,综合出露破裂面的擦痕所指示的滑移向量,并对比根据矩张量解获得的震源深度的滑移向量,得出以下认识:(1)破裂南段的地表滑移向量的方位角明显小于震源深度滑移向量的方位角,表明在破裂从震源向地表传播过程中破裂面上的滑移向量发生了逆时针旋转;(2)滑移方位角向北东方向逐渐增大,表明地平面上水平滑移向量表现出顺时针旋转的趋势,而且在破裂向北东方向传播过程中近地表的走滑分量逐渐减小而倾滑分量逐渐增大;(3)几乎在每一个观测点倾滑分量都大于走滑分量,表明汶川地震的破裂方式在任何地点都是以逆冲运动为主;(4)破裂面倾角在10.4°～64.7°,平均值为41°,与天然破裂露头和探槽揭示的结果基本一致;(5)滑移向量沿破裂带的分布显示,走滑分量中段大而两端小,倾滑分量则相反,中段小两端大。     

3D form line construction by structural field interpolation (SFI) of geologic strike and dip observations

Interpreting and modelling geometries of complex geologic structures from strike/dip measurements using manually-drafted structural form lines is labour intensive, irreproducible and inherently limited to two dimensions. Herein, the structural field interpolation (SFI) algorithm is presented that overcomes these limitations by constructing 3D structural form lines from the vector components of strike/dip measurements. The SFI interpolation algorithm employs an anisotropic inverse distance weighting scheme derived from eigen analysis of the poles to strike/dip measurements within a neighbourhood of user defined dimension and shape (ellipsoidal to spherical) and honours younging directions, when available. The eigen analysis also provides local estimates of the plunge vector and associated Woodcock distribution properties to assure plunge-normal structural form line reconstruction with unidirectional propagation of form lines across fold and fan structures. The method is advantageous for modelling geometries of geologic structures from a wide range of structurally anisotropic data. Modelled vector fields from three case studies are presented that reproduce the expected bedding-foliation geometry and provide reasonable representation of complex folds from local to regional scales. Results illustrate the potential for using vector fields to support geologic interpretation through the direct visualization of geometric trends of structural features in 3D.     

A simple graphical method for estimating the components of the fault-slip vector

A simple calculation of the components of the total slip vector (D) on a fault plane allows the relationships between the magnitudes of the three slip components of D, the lateral horizontal displacement (L), the transversal horizontal displacement (T) and the vertical offset (V), to be determined. The contribution of each slip component to the total slip can be plotted jointly in a ternary diagram, assuming a unit length of the vector modulus and a suitable normalization for D. Each component equals the magnitude of D at the vertices of the diagram, hence it is possible to estimate the percentage contribution of each slip component to the total movement of a fault or a set of faults. The dip of the fault surface and the rake of the slickenlines are the input data required for displaying L, T and V relationships in these diagrams. This information may be useful in the analysis of movement geometry for different fault populations and in the determination of D by measuring just one of the fault-slip components, such as vertical slip associated with a fault scarp.     

Role of shear along horizontal plane in the formation of helicoidal structures

An unusual structural paragenesis, complicated by brachyanticlines, is revealed for the first time in the sedimentary cover of the West Siberian Plate by 3D seismic surveying. These are linear (in plan view) systems of en-echelon arranged low-amplitude normal faults related to wrench faults in the basement. On different sides off a wrench fault, the planes of normal faults dip in opposite directions, forming a helicoidal structure that resembles the blades of a propeller. In the section parallel to the wrench fault, the boundaries of the beds and normal fault planes dip in opposite directions as well. In the section across the strike of the normal faults converging toward the basement, the beds take the shape of an antiform with a crest sagged along the normal faults (flower structure). This structural assembly was formed as a result of interference of stress fields of horizontal shear in the vertical plane (induced by faulting in the basement) and in the horizontal plane (caused by gravity resistance of the cover). In this case, the displacements along the normal faults develop in both the vertical and, to a greater extent, horizontal directions, so that the faults in cover are actually characterized by normal-strike-slip kinematics. The regional N-S-trending compression of the West Siberian Plate is the main cause of shearing along the NW- and NE-trending faults in the basement, which make up a rhomb-shaped system in plan view. Petroliferous brachyanticlines, whose axes, notwithstanding tectonophysical laws, are oriented in the direction close to the maximum compression axis, are known in the large wrench fault zones of Western Siberia. Our experiments with equivalent materials showed that a local stress field arising at the ends of echeloned Riedel shears within a wrench fault zone may be a cause of the formation of such brachyanticlines. The progressive elongation of Riedel shears leads to the corresponding elongation of the brachyanticlines located between their ends. The performed study has shown that the known types of interference of elementary geodynamic settings such as horizontal shear along the vertical plane + horizontal compression (transpression) and horizontal shear along the vertical plane + horizontal extension (transtension) may be supplemented by combination of horizontal shears along the vertical and horizontal planes, resulting in tectonic lamination. By analogy, we propose to name this type of interference of elementary shear settings translamination. Petroliferous helicoidal structures arise in the given geodynamic setting of translamination.     

层状岩质边坡破坏模式及稳定性的数值分析

运用FLAC3D模拟层状岩质边坡的破坏模式,并采用强度折减法分析结构面倾角与稳定性之间的关系。结果表明:①水平层状边坡坡顶变形破坏早于坡面和坡脚,形成拉破坏区。当结构面倾角较小时,顺倾向层状边坡主要发生滑移破坏;当结构面倾角较大时,发生弯折-溃曲破坏。直立层状边坡主要发生弯曲-板间拉裂-塌落破坏。逆倾向层状边坡的破坏形式为对于小倾角为滑移破坏,对于陡倾角为倾倒破坏。②对于顺倾向边坡,整体安全系数随结构面倾角先减小后增大,呈现两头高、中间低的形态,在倾角为30°时安全系数最小;对于逆倾向边坡,曲线呈现增大-减小-增大的态势,并且大部分高于顺倾向边坡曲线,符合实际情况。③所采用的低强度弹塑性单元能够比较真实地模拟软弱结构面的变形。     

利用成像测井资料分析汶川地震断裂带科学钻探3号孔(WFSD-3)裂缝特征

WFSD-3孔是汶川地震断裂带科学钻探主要钻孔之一,全井段(终孔深度1 502.30 m)实施了连续取心(累计取心进尺1 548.44 m)和测井作业。采集到的成像测井资料包含丰富的原位地质特征信息,对裂缝、破碎带识别和构造应力场分析具有重要作用。利用该钻孔的电阻率成像测井和超声成像测井资料,结合岩心资料进行了裂缝特征分析。结果表明,WFSD-3孔岩层高角度斜交缝最为发育,其次是低角度斜交缝,垂直缝和水平缝极少;25～200 m和900～1 000 m深度范围内裂缝尤为发育;不同深度的裂缝倾向存在明显差异:410 m之上主要分布于260°～290°,410～730m集中于330°～360°,730～960 m主要分布于210°～240°,960～1 185 m与410 m之上基本一致,1 410～1 450 m与前述各深度段不同,集中于180°～200°;裂缝与破碎带、层理密切相关,宏观分布受构造控制。     

Recent stress field and neotectonics in the Eastern Jura Mountains,Switzerland

The regional stress field and its local variation were determined for the northern part of central Switzerland (Fig. 1) by using overcoring techniques (doorstopper, triaxial strain cell) and observations of breakouts in deep boreholes. The results are compared with fault plane solutions of earthquakes and with the orientation of horizontal stylolites.In the northern part of central Switzerland the NW-SE-orientation of the maximum horizontal stress (S_H) which is characteristic for Central Europe was observed only in the crystalline basement. In the Folded Jura and south of it in one well the greatest principal horizontal stress above the Triassic decollement horizon is oriented approximately in a N-S to NNE-SSW direction.This direction persists into the western Tabular Jura and the southernmost Rhine Graben. Only in the eastern part of the Tabular Jura the greatest principal horizontal stress shows a NNW-SSE to NW-SE orientation. Comparison of the near surface stress field as determined by in situ stress measurements and borehole breakouts with the directions of horizontal stylolites generated during the evolution of the Folded Jura, indicates that the orientation of the recent stress field near the earth's surface is the same as that which prevailed during the Upper Miocene to Lower Pliocene.The central part of northern Switzerland is therefore the first area in Central Europe where it is possible to demonstrate that the near surface stress field is decoupled from that in the crystalline basement. The difference in the orientation of the greatest principal horizontal stress is about 50°.     

Determination of a spatial fracture network

Hydraulic and mechanical properties of rock masses are largely dependent upon the distribution and variation of fracture areas. Consequently, it is essential to develop good models for analysing the spatial geometric characteristics of fracture fields. The latter, however, may be ascertained only by measuring open cuts or boreholes. The most important factor from the analysis of areas is to be able to determine spatial networks. Without an exhaustive and expensive inventory, the classical methods involving measurements of dip and strike do not provide either a correct estimation of fracture orientations in space, or a good estimation of their distribution. At present, automatic analysis of linear traces on large plane areas is used by several authors^1,2 and the field strike and dip measurements are applied only for verification purposes. This method of computing spatial networks is especially suitable because of the nature of data so obtained, since it permits various simulations on computed networks. It is based upon classical methods of vector geometry that involve a measurements data matrix and provide a means of solving a double system of equations. Results are the directional spatial fracture densities. The object of this paper is to define the geometric parameters of this problem and the computation method, to present an example and, finally, to conclude with critical comments supported by certain proposals for a variety of applications.     

基于自适应加权广义逆矢量方向滤波估计地震同相轴倾角

在断层等不连续位置处估计局部地震同相轴的方向较为困难。在断层、地质边界等处,同相轴不连续,结构不稳定,导致倾角等地震属性变化大。为了克服这个问题,本文提出一种自适应加权广义逆矢量方向滤波方法,通过光滑经过逆转的梯度矢量场来估计局部地震倾角剖面。计算过程是：首先由改进的梯度算子计算梯度矢量场;然后设置垂直向下为参考方向,对梯度矢量场进行逆转;最后通过自适应加权广义逆矢量方向滤波方法估计倾角剖面。自适应加权广义逆矢量方向滤波方法在不连续位置处具有传统倾角估计方法不具有的高分辨率,能够稳健高效地估计地震数据的倾角剖面。理论模型和实际地震资料处理得到的局部地震倾角剖面表明,所提出的方法能够高分辨率地保持构造细节信息。     

A preliminary study of the dynamics of a modern draa, Algodones, southeastern California, USA

A transverse crescentic draa in the Algodones dune field, California, was monitored for a year using surface process mapping, aerial photography and supplemental wind measurement. The draa is oriented by the long-term resultant wind, whereas its superimposed features are in equilibrium with the bedform-modified secondary airflow. Surface airflow and the movement of superimposed bedforms is typically oblique or parallel to the draa brinkline, particularly on the lee slope. Comparison of measurements of draa movement and sand deposition on the lee slope, with expected rates calculated from wind data and draa size, confirm that there is a significant component of sand flow parallel to the draa brinkline. The internal structure being generated at the base of the draa lee slope is inferred from the surface processes active there. Within the space of a kilometer two types of compound cross-strata, separated by an area of simple cross-strata, are being produced. This has significant implications for interpretations of ancient aeolian strata. Variations in internal structure types found in lateral sequence may be generated by one complex bedform, and these cross-strata may be simple or compound. Second-order bounding surface orientations indicate resultant primary palaeowind directions; compound cross-strata dip directions indicate secondary flow conditions. The existence of cross-strata dip directions oblique or perpendicular to the second-order surface indicates longitudinal secondary flow on the lee face, but not necessarily a longitudinal or oblique draa. Without further detailed knowledge about various draa configurations and behaviour, stratification attributed to draas can be used only to interpret activity on the lower draa lee face.     

小角度成层倾斜场地动力响应分析的大型振动台试验研究

基于大型振动台试验,研究小角度成层倾斜场地在倾向、走向、垂向、坡面方向、坡面垂直方向上的地震动响应特征,以及地层倾角对场地反应谱的影响规律。试验结果表明：在倾向方向,随着地层倾角的增大,倾斜场地的放大效应增大;在走向方向,当地层倾角小于12.5°时,放大效应强于水平成层场地;在倾向、走向方向,当地层倾角达到12.5°时,场地对反应谱短周期T≤0.1 s部分具有放大作用,对T＞0.1 s的部分具有削弱作用;在垂向方向,随着地层倾角的增加,倾斜场地的放大效应强于水平成层场地;在坡面方向,倾斜场地的放大效应弱于水平成层场地;在坡面垂向方向,倾斜场地的放大效应强于水平成层场地。该研究成果对小角度成层倾斜场地上的建筑物抗震设计具有一定的参考意义。     

陡倾坡外弱面控制的斜坡滑移-剪损变形破坏模式

斜坡变形破坏和稳定性分析是各类工程建设中高度关注的问题.采用实例调查、理论分析、数值计算等技术方法,以雅砻江某水电工程坝址区右岸顺层岩质斜坡为例,研究总结了斜坡发育滑移-剪损变形破坏的成因机理、发育特征以及与弱面倾角和发育深度的关系.研究表明,滑移-剪损变形破坏通常发育在力学性能相对较差的薄层、互层状结构的顺层岩质斜坡或斜坡强-弱风化带内.斜坡发育滑移-剪损变形破坏与陡倾坡外弱面的倾角和发育位置密切相关.倾角在45°~65°之间或距斜坡表部水平距离小于80 m的弱面对斜坡发育滑移-剪损变形破坏的影响控制作用明显,并且弱面距斜坡表部水平距离比弱面倾角对斜坡发育滑移-剪损变形破坏的影响控制作用更强.研究成果可补充完善岩质斜坡变形破坏类型,具有重要的工程意义和实践价值.      

Late‐Holocene wind‐field evolution at the southern Baltic coast as revealed by GPR data from the Mrzeżyno dunefield,NW Poland

All dunes within a dunefield usually show a similar response to the wind regime. However, in the late‐Holocene coastal dunefield in NW Poland the surface topography suggests that slipface orientation may vary significantly with distance from the coast, which is rarely reported in the literature. The dunefield was stabilized with forest in the mid‐19th century, preserving a unique record of atmospheric circulation in the South Baltic region at the end of the Little Ice Age. To elucidate the Holocene processes occurring along the study site a pseudo‐3D GPR data set was collected. Six grids of parallel GPR lines combined with 400 m of GPR profiles (2D) were collected across the dunefield and displayed using GOCAD for interpretation and geostatistical analysis. The geophysical data revealed that the larger dunes almost entirely consist of steeply laminated facies. Most importantly the pseudo‐3D data, supported by geostatistical estimates of strata dip directions, revealed the existence of three zones parallel to the coastline with the mean dip direction almost perpendicular to the coast in the northern coastal strip and almost parallel to it in the southern part. Spreads of the dip directions in pseudo‐3D GPR data sets recorded on the stoss slopes of dunes and crests suggested initial deposition on transverse dunes, which later were transformed into barchanoid dunes. This can probably be linked to changes in the wind regime, i.e. reduction in velocity of the northern and northwestern winds leading to reduced sediment supply from the coast. While the data provide a new interpretation of Holocene dunefield dynamics at this site, they also suggest that the minimum number of pseudo‐3D GPR grids required to establish general trends using geostatistical analysis should be at least 10, with even more data needed at larger, or more complex dunefields.     

Structural data extraction from stereographic or equal area projections using image processing

Structural orientation data from the field are commonly represented on stereographic or equal area projections, and the raw data are usually not accessible to other workers for integration into larger datasets or comparative study. We have developed a technique for retrieving raw planar structural data from the above projection types using image analysis techniques. The principle involves two steps: (i) conversion of the foliation pole plots in the image into Cartesian coordinates, and (ii) inversion of the data soobtained into a geographic reference frame. The process involves either scanning of the target projection, or cropping of the projection figure from a pdf file using image processing software, followed by tracing of the data points and generation of an image file that can be subsequently input into a MATLAB supported program, “StrucExtract”. The results from the program are displayed in the form of dips and dip directions. The program has been tested with synthetic and natural examples, and may be invaluable for retrieval of structural data from older structural literature.     

重力场数据处理关键技术应用

重力场数据处理涉及的内容及相关处理软件很多。本文利用我省"十一五"期间自主开发的"SPDCASPS"软件对已知矿区重力场数据进行了曲化平、向下延拓、局部场和背景场的分离等关键技术的应用。应用结果表明,经三种重力场数据处理关键技术的应用,取得了良好的地质效果。     

Composite detrital and thermal remanent magnetization in tuffs from Aeolian Islands (southern Tyrrhenian Sea) revealed by magnetic anisotropy

This paper reports on the complex relation between rock emplacement and remanence acquisition in tuffs deposited by pyroclastic density currents, disclosed by systematic measurements of the anisotropy of magnetic susceptibility and natural remanent magnetization (NRM). Thermal demagnetization shows that the NRM consists of two components with different blocking-temperature spectra. The direction of the low-temperature component is consistent with the geocentric axial dipole value, whereas the high-temperature component has dispersed directions. The magnetic fabric is oblate, the magnetic foliation is close to the bedding and the lineations are generally dispersed along a girdle within the foliation plane. The directions of the magnetic lineation and the high-temperature remanence component of individual specimens are close to each other. This correspondence suggests that the high blocking-temperature grains acquired a remanence aligned to their long dimension before deposition, while cooling within the explosive cloud and the moving pyroclastic current. Thereafter, during deposition, the traction processes at the base of the current oriented the grains along the flow direction and affected both fabric and high-temperature remanence. This NRM component results from mechanical orientation of previously magnetized grains and is thus detrital in origin. A second, thermal component was then acquired during the cooling of the low blocking-temperature grains after deposition. These results show that NRM in fine-grained pyroclastic rocks is affected by the Earth’s magnetic field as well as the emplacement processes and that magnetic fabric data are essential to unravel its complex nature.     

Rock bursting phenomena in a superficial rock mass in southern central Sweden

SummaryRock Bursting Phenomena in a Superficial Rock Mass in Southern Central Sweden Rock bursting phenomena have been observed in connection with excavation of two shallow sited rock tunnels in the Forsmark area in southern, central Sweden. The tunnels are excavated in gneiss granite and the rock cover varies between 5–15 m. The stress release occurred locally in the tunnels which indicates significant variations in rock stresses.In situ measurements of the triaxial stress situation carried out in the area show values of the principal stresses in the horizontal plane of a magnitude in excess of 20 MPa in the superficial rock mass. The rock stress measurements indicate a clear correspondence of the directions of the highest compressive stress and the rock foliation. There is also a clear correspondence between the directions of a vertical and a horizontal joint set and the highest compressive stress. The three-dimensional stress tensor indicates that the horizontal fractures are exposed to the smallest closure pressure, i. e. the vertical stress, which set should therefore be most open. The horizontal joint set in the superficial rock mass has a significantly higher value of aperture than the vertical sets which may be a consequence of the stress situation in the rock mass. The rock bursting is not attributed to a breakage of the rock matrix, but merely to a propagation of already existing small fractures and potential fractures. The stress situation in the Forsmark area may be regarded as an uniaxial stress field in the horizontal plane.With 12 Figures     

地磁场数据处理关键技术应用

地磁场数据处理涉及的内容颇多。实际工作中,应该针对地磁场异常特点和需要解决的地质问题,必须抓住数据处理的关键技术。本文通过对已知矿区地磁场数据的曲化平、向下延拓、局部场和背景场的分离等关键技术的应用,结果表明,经三种地磁场数据处理关键技术的应用,取得了良好的地质效果。     

井下超声波成像提取岩层节理产状信息

井下超声波电视能够提取工程勘察中岩层产状信息并完成井眼岩层3维（3D）信息图。在某地下石油储备库的工程勘察中,完成7口井的超声波电视成像测井,以其中的钻孔ZK10为例,研究并分析了井下超声波电视成像所获取的结构面倾角和倾向、岩体完整性、结构面线密度、3D岩心图和井壁图以及节理极坐标图、节理方位频率图。统计分析结果表明：井下超声波电视技术可以对钻孔所揭露的岩体质量进行初级定量化分析,为区域稳定分析、工程岩体质量分级及工程设计提供有价值的基础数据,对提高工程勘察质量具有重要意义。