异源点集下三维断层建模方法

用于断层建模的异源点集,因其具有不同特征而影响到建模途径和方法。以三维断层模型的3D GIS应用为目标,分析了原始地质观测点集、钻孔点集、剖面点集和平面点集的信息构成特征,提出了异源点集下三维断层数据模型的构建方法和流程。原始地质观测点集和钻孔点集间接表达了断层的空间特征,通过空间平面、曲面和曲线构建了断层面、地层面和断层交线的数学模型,并借助于连续模型的离散化,实现了独立断层建模点集的重建;剖面点集和平面点集可直接用于断层建模。以双TIN面模型为例阐述了断层三维模型的构建方法。      

塔里木盆地走滑带碳酸盐岩断裂相特征及其与油气关系

通过露头与井下资料的综合分析,塔里木盆地奥陶系碳酸盐岩走滑断裂带断裂相具有多样性,根据内部构造发育程度可以分为断层核发育、断层核欠发育两类。露头走滑带断层核部以裂缝带、透镜体、滑动面等断裂相发育为特征,断裂边缘的破碎带发育裂缝带、变形带。裂缝带主要分布在断层核附近50m的破碎带内,裂缝多开启,渗流性好。断裂核部透镜体发育,在破碎带也有分布,破碎角砾组合的透镜体多致密。滑动面具有平直截切型、渐变条带型等两种类型,多为开启的半充填活动面。变形带多为方解石与碎裂岩充填,破碎带局部部位裂缝与溶蚀作用较发育。利用地震剖面、构造图、相干图等资料可以判识塔里木盆地内部奥陶系碳酸盐岩走滑断裂相的特征及其发育程度,沿走滑断裂带走向上断裂相具有分段性与差异性,根据渗流性可以定性区分高渗透相、致密相区。沿断裂带高渗透相区是碳酸盐岩缝洞体储层发育的有利部位。断裂相的横向变化造成油气分布的区段性,形成高渗透相输导模式、致密相遮挡模式等两类成藏模式。走滑断裂带碳酸盐岩断裂相的特征及其控藏作用对油气勘探开发储层建模具有重要意义。     

三峡九湾溪断裂带断层泥的特征及分段研究

在分析研究前人资料的基础上,结合野外实际工作和室内研究,从三峡九湾溪断裂带内各段断层泥的特征研究分析出发,探讨了九湾溪断裂带的分段性,结果发现:九湾溪断裂带存在明显的分段性,中段和南段地震活动频度和强度大,应变强度高,断层活动以粘滑为主,其分维值在2.2～2.6之间,石英碎砾表面发育各种撞击揳入现象、线状擦痕等显微构造;北段地震活动频度和强度较小,应变强度弱,断层活动以稳滑为主,兼具粘滑,其分维值在2.4～2.8之间,石英颗粒多为磨圆球砾,其刻蚀形貌以裂而不破现象为代表;中段是九湾溪断裂带最有可能发生诱发地震的地段和我们的重点监测地段。     

塔里木地块北部横向构造及断条模式

提要:塔里木盆地北部与西南天山毗邻区域发育的浅部构造系统是喀什坳陷—柯坪塔格逆冲推覆构造带—库车前陆坳陷冲断带,它们在山前以1～3排方式平行造山带展布,是一套由底板滑脱断坪和断坡构成的逆冲推覆系统,具有明显的横向分段性;该区深部构造系统却是一组以北西向为主展布的断条构造,并可以划分出4个一级断条构造。笔者从构造层变形、重力、航磁异常场及天然地震平面分布密度等几个角度,对区内深部横向构造系统进行了研究,并利用天然地震深度/频次统计结果,识别出圈层拆离解耦面的深度分布,进而探讨了深部横向构造运动的岩石圈地壳圈层归属,对区域北西向横向构造系统做了详细的论述。本文选择喀拉玉尔滚断裂带和库尔勒断裂带进行重点解剖,研究了横向构造与浅部构造的转换关系,最终提出了塔里木地块北部断条构造几何学运动学模式,即塔里木地块深部向北西方向运动受到古西南天山的阻挡性约束而“被迫”俯冲,岩石圈地壳发生拆离解耦,原有的横向构造——北西向构造带(断裂带)被激活,使得俯冲系统以断条为单位进行俯冲;在俯冲过程中,岩石圈地壳的拆离及横向构造被激活的方向是从约束体(南天山)附近开始,向塔里木盆地(北部)的腹地方向推进,因此,越是接近造山带区域,断条俯冲状态的差异越明显,横向分段性越突出;而正因为塔里木地块以断条形式向南天山的俯冲行为,使得山前坳陷冲断带(及天山)发生分段。     

基于水平集的复杂三维地层模型建模及在地下工程中的应用研究

由于地质体存在断层、尖灭、出露等复杂地质现象,在三维地层建模时,为了表达这些现象,无论是面模型还是体模型或混合模型,都存在空间分割或曲面求交线的问题。由于地质体拓扑关系的复杂性、数据误差以及计算机精度问题,使得这些模型在实际建模过程中常常失效。运用水平集理论可以有效解决这一问题,水平集用隐函数表达曲面（或超曲面）,可以实现复杂地质体的表达及并、交、差等拓扑运算。在三维地层建模中,插值生成各种地质界面后,用水平集表达这些地质界面,利用水平理论完成各种复杂的拓扑操作,建立以水平集表达的三维地层模型。在此基础上,插入水平集表达的各种工程活动界面,利用水平集理论进行拓扑操作,可构建各种工程活动后的地层模型。再利用Marching cube算法抽取各种地质界面或工程活动界面,构建可用于实时可视化或用于工程评估（如有限元计算）等的NMTINF-BR地层模型或工程活动后的NMTINF-BR地层模型。     

The corrugation and bifurcation of fault surfaces by cross-slip

The formation of corrugated and stepped fault surfaces by conservative shear is discussed. It is suggested that fault boundaries which lie parallel to the direction of shear are able to deviate from a plane as a fault spreads, like dislocations that cross-slip in a crystal. It is explained how deviations from planar shear can bifurcate fault surfaces and create en-échelon faults. The residual stresses produced by such deviations are tensions and compressions in the direction of shear.     

Relation Between Level Set and Truncated Pluri-Gaussian Methodologies for Facies Representation

Truncated Gaussian/pluri-Gaussian representations and level set representations are two methodologies for implicit representation that can be applied to large-scale subsurface structures. Identification of facies in a petroleum reservoir from dynamic well data is a common field of application for these representation methodologies. No comparison of the methodologies has appeared in the literature to date. The paper seeks to improve on this situation by comparing selected level set and truncated Gaussian/pluri-Gaussian representations in detail. Strong similarities are found, in the sense that every truncated Gaussian/pluri-Gaussian representation considered has a level set counterpart. Furthermore, the transition from a truncated Gaussian/pluri-Gaussian representation to the corresponding level set representation is easily accessible. In addition to the comparison aspect, this paper also introduces a novel level set representation—the hierarchical level set representation—that removes a difficulty present in existing level set representations in a shape estimation setting. It is shown that the hierarchical level set representation corresponds to a well-known truncated pluri-Gaussian representation.     

Regional dome evolution and its control on ore-grade distribution: Insights from 3D implicit modelling of the Navachab gold deposit,Namibia

We introduce a novel approach to analyse and assess the structural framework of ore deposits that fully integrates 3D implicit modelling in data-rich environments with field observations. We apply this approach to the early Palaeozoic Navachab gold deposit which is located in the Damara orogenic belt, Namibia. Compared to traditional modelling methods, 3D implicit modelling reduces user-based modelling bias by generating open or closed surfaces from geochemical, lithological or structural data without manual digitisation and linkage of sections or level plans. Instead, a mathematically defined spatial interpolation is used to generate 3D models that show trends and patterns that are embedded in large drillhole datasets. In our 3D implicit model of the Navachab gold deposit, distinctive high-grade mineralisation trends were identified and directly related to structures observed in the field. The 3D implicit model and field data suggest that auriferous semi-massive sulphide ore shoots formed near the inflection line of the steep limb of a regional scale dome, where shear strain reached peak values during fold amplification. This setting generated efficient conduits and traps for hydrothermal fluids and associated mineralisation that led to the formation of the main ore shoots in the deposit. Both bedding-parallel and highly discordant sets of auriferous quartz-sulphide veins are interpreted to have formed during the later lock-up stage of the regional scale dome. Additionally, pegmatite dykes crosscut and remobilise gold mineralisation at the deposit scale and appear to be related to a younger joint set. We propose that kilometre-scale active folding is an important deformation mechanism that influences the spatial distribution and orientation of mineralisation in ore deposits by forming structures (traps and pathways for fluids) at different preferred sites and orientations. We also propose that areas that experience high shear strain, located along the inflection lines of folds can act as preferred sites for syn-deformational hydrothermal mineralisation and should be targeted for regional scale exploration in fold and thrust belts. Our research also suggests that examination of existing drillhole datasets using 3D implicit modelling is a powerful tool for spatial analysis of mineralisation patterns. When combined with fieldwork, this approach has the potential to improve structural understanding of a variety of ore deposits.     

三肇凹陷断层垂向分段生长与扶杨油层油源断层的厘定

野外观察、典型地震解剖和物理模拟证实, 只要岩石存在能干性差, 断层具有典型垂向分段生长特征.岩石能干性差异决定断层演化历经3个阶段: 下部断层形成、断裂上下分段和贯通性断裂形成.以三肇凹陷为例, 应用岩石力学特征和断层相关褶皱理论, 证实扶杨油层砂岩层段普遍形成断裂, 青一段泥岩阻止下部断层向上传播, 从而形成断层垂向分段生长现象.结合油藏精细解剖表明, 三肇凹陷扶杨油层"倒灌"运移的主要油源断层是成藏期活动、沟通源储且垂向分段生长的断裂.      

徐家围子断陷深层结构形成与演化的探讨

徐家围子断陷是松辽盆地北部深层最重要的含气断陷。本文根据最新地震和钻井资料综合研究成果,结合前人研究认识,探讨了徐家围子断陷的形成与演化,认为徐家围子断陷在区域构造位置上处于东北地区松辽盆地的北部,处于松辽盆地区域莫霍面隆起的西斜坡上,特定的位置决定了其形成演化的各个阶段受东北地区区域板块构造运动的控制,与松辽盆地的形成与演化具有一定的关系。徐家围子断陷形成与演化可以划分为5个阶段,是太平洋构造域板块间的相互作用和深部热力作用发育演化的结果,这种过程造就了现今南北分块、东西分带、凹隆相间、构造复杂和沉积岩与火山岩并存的地质结构。徐家围子断陷深层现今构造格局受控于近NW向的徐中断裂、近SN向的徐西断裂、近NW向的徐东断裂带和四组近NE向的断裂。本次研究,对以往无法解释的很多地质现象给予了很好的解释,得出了一些崭新的认识,也从宏观上为寻找天然气资源给予了有力指导。     

Sequence Stratigraphic Features of Mesozoic－Cenozoic Nonmarine Sediments in Eastern China

ＳｅｑｕｅｎｃｅＳｔｒａｔｉｇｒａｐｈｉｃＦｅａｔｕｒｅｓｏｆＭｅｓｏｚｏｉｃ－ＣｅｎｏｚｏｉｃＮｏｎｍａｒｉｎｅＳｅｄｉｍｅｎｔｓｉｎＥａｓｔｅｒｎＣｈｉｎａ￥ＸｕＨｕａｉｄａ；ＷｅｉＫｕｉｓｈｅｎｇ（ＤｅｐａｒｔｒｍｅｎｔｏｆＥｎｅｒｇｙＲｅｓｏ...     

Structural Traces of Secondary Faults (Fractures) along the Main Faults and Their Reliability as Kinematic Indicators

Secondary/minor structures occurring along the main fault surfaces are important indicators for judging the kinematic characteristics of faults. However, many factors can lead to the formation of these structures, which results in the difficulty for rapid judgment and application in the fields. A series of secondary faults/fractures developed due to the movement of main faults are the most important and widespread phenomena in the scope of brittle deformation. The morphology of the main fault surfaces is various, and former researchers mainly discussed the structures on the main even fault surfaces. However, the fluctuation of fault surfaces is the intrinsic character of the faults, and the intersection between the main fault and secondary faults/fractures can produce a series of kinematic indicators on the main fault surfaces. Based on previous studies and our observations, i.e. the structural traces of the P, R, R’, T and X shears/faults along the main faults, some indicators which are rarely reported previously, are described in the paper. Furthermore, their reliabilities are also discussed, and more practical and reliable criteria are brought forward. We suggest that the simple application of congruous and incongruous steps without knowing their exact origins should be abandoned in the fields, and several types of indicators along one fault surface should be checked with each other as much as possible. Meanwhile, the origins of some other arcuate indicators on the fault surfaces are also discussed, and new models are brought forward.     

阿尔金活断层的古地震与分段

活断层的破裂过程往往是通过多个独立破裂段落的组合而进行的。活断层的分段就是对断层上各个稳定的独立破裂单元的识别。这是深入认识断层的活动习性特征和进行潜在地震危险性评估的重要途径。古地震资料是识别破裂单元最直接的依据。本文在根据障碍构造对阿尔金活断层进行几何结构分段的基础上，结合古地震资料的分析进行了破裂分段的识别，划分了１１个破裂段。在约２００００年来各个段落活动时段、强度、古地震复发间隔都有所不同。总趋势是由西向东逐渐减弱。在约４５００—２０００ａＢ．Ｐ．存在一个各个段落都比较活跃的阶段。这与沿青藏高原边缘的一些活动带上揭露出的情况很相似，看来这是与整个高原的构造活动在这个时段曾有所加强有关。     

Continuous Facies Updating Using the Ensemble Kalman Filter and the Level Set Method

We present a methodology based on the ensemble Kalman filter (EnKF) and the level set method for the continuous model updating of geological facies with respect to production data. Geological facies are modeled using an implicit surface representation and conditioned to production data using the ensemble Kalman filter. The methodology is based on Gaussian random fields used to deform the facies boundaries. The Gaussian random fields are used as the model parameter vector to be updated sequentially within the EnKF when new measurements are available. We show the successful application of the methodology to two synthetic reservoir models.     

Crustal structure, fault segmentation, and activity of the Median Tectonic Line in Shikoku, Japan

We conducted a seismic tomographic analysis to estimate the crustal structure beneath the Shikoku and Chugoku regions in Japan. The Philippine Sea slab (PHS slab) subducts continuously in a SE–NW direction beneath this region, and the crustal structure is complex. Furthermore, the Median Tectonic Line (MTL), one of the longest and most active arc-parallel fault systems in Japan (hereafter, the MTL active fault system), is located in this area, and the right-lateral strike–slip movement of this fault system is related to the oblique subduction of the PHS slab. The MTL active fault system has ruptured repeatedly during the last 10 000 years, and has high seismic potential. Our tomographic analysis clarified the heterogeneous crustal structure along the MTL active fault system. This fault system in Shikoku can be divided into two segments, an east segment and a west segment, on the basis of the velocity structure. This segmentation model is consistent with other such models that have been determined from geological and geomorphological data such as fault geometry, slip rate, and faulting history. This consistency suggests that the surface characteristics of the MTL active fault system are related to structural properties of the crust. In particular, a prominent low-velocity (low-V) zone is present in the lower crust beneath the east segment. Our tomographic images show that the lower crust structure beneath the east segment is obviously different from that of the other segment. Furthermore, this low-V zone may indicate the presence of fluid, possibly related to dehydration of the PHS slab. As the presence of fluid in the lower crust affects the activity of the fault, stress accumulation and the fault failure mechanism may differ between the two segments of the MTL active fault system.     

准南前陆冲断带构造分段及其与油气关系

通过大量地面地质、地震和非地震资料的综合解释,分析了准南前陆冲断带的构造分段特征.认为乌鲁木齐-米泉断裂和红车断裂等盆内大型基底断裂是控制构造分段的主要因素,其形成演化决定了构造分段特征的出现或消失;盆缘调节构造是造成次级构造分段的主要原因.准南前陆冲断带构造分段对该区主要烃源岩与生储盖组合的形成及分布、对油气成藏期等方面具有重要的控制作用.     

Geometry and segmentation of an evaporite-detached normal fault array: 3D seismic analysis of the southern Bremstein Fault Complex,offshore mid-Norway

The Halten Terrace, offshore mid-Norway, is underlain by a Triassic evaporitic package that is rheologically weak, and led to decoupling of fault systems during Middle Jurassic to Early Cretaceous rifting. We use 2D and 3D reflection seismic data, constrained by wells, from the southern Bremstein Fault Complex of the Halten Terrace to map faults and key stratigraphic horizons, and analyse throw variations along faults, allowing us to constrain patterns of fault segmentation and linkage within the complex. The Bremstein Fault Complex has an overall tilted monoclinal geometry with localised fault systems at base salt level associated with overlying, highly distributed systems of normal faults. Vertical strain partitioning across the evaporite package means that sub-evaporite and supra-evaporite fault populations acted as semi-independent fault systems. Supra-evaporite faults are partly gravity-driven, and controlled by sub-evaporite faulting and consequent tilting of the evaporitic package. This behaviour leads to a wide variety of possible vertical linkage patterns of faults across the evaporite package. A greater variety of lateral segment linkage patterns occurs in evaporite-detached normal fault systems than in normal fault systems developed in the absence of evaporite units. Segment boundary styles can also be modified by migration of evaporite. Some segment boundaries are associated with a footwall anticline and hanging-wall syncline, in contrast to the footwall synclines and hanging-wall anticlines widely described in studies of normal fault systems.     

Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.     

Geophysical evidence on segmentation of the Tancheng-Lujiang fault and its implications on the lithosphere evolution in East China

The north–south trending Tancheng-Lujiang (Tanlu) fault belt extends from northeast China to the Dabie–Sulu orogenic belt, for a length of more than 3000 km. This fault belt probably has close links with the lithosphere evolution, seismic activity and mineral resource concentration in East China. Surface geological mapping and studies on sedimentation and basin formation have indicated segmentation at the southern, middle and northern domains of the fault. Here we employ geophysical constraints to evaluate these fault segments. Unlike previous geophysical studies focused on laterally varying crust/mantle seismic velocity structure across the fault, in this study we have integrated a variety of geophysical data sets, such as crustal P-wave velocity, earthquake occurrence and released seismic energy, seismogenic layer thickness, surface heat flow and geothermal field, to understand the deep structure and strength of the lithosphere along the Tanlu segmented fault belt. The results demonstrate remarkable crustal-scale north-to-south segmentation this major fault. The geophysical evidence and some geochemical constraints suggest that the Tanlu fault belt probably served as a channel for melt and fluid percolation, and exerted a significant control on the lithosphere evolution in East China.