Separation of magnetic anomalies into induced and remanent magnetization contributions

Inversion of magnetic data is complicated by the presence of remanent magnetization, and it provides limited information about the magnetic source because of the insufficiency of data and constraint information. We propose a Fourier domain transformation allowing the separation of magnetic anomalies into the components caused by induced and remanent magnetizations. The approach is based on the hypothesis that each isolated source is homogeneous with a uniform and specific Koenigsberger ratio. The distributions of susceptibility and remanent magnetization are subsequently recovered from the separated anomalies. Anomaly components, susceptibility distribution and distribution of the remanent and total magnetization vectors (direction and intensity) can be achieved through the processing of the anomaly components. The proposed method therefore provides a procedure to test the hypotheses about target source and magnetic field, by verifying these models based on available information or a priori information from geology. We test our methods using synthetic and real data acquired over the Zhangfushan iron-ore deposit and the Yeshan polymetallic deposit in eastern China. All the tests yield favourable results and the obtained models are helpful for the geological interpretation.     

三维地震勘探技术在新疆沙吉海井田中的应用

针对沙吉海井田以往未开展过地震勘探,煤层多、局部间距小特点,介绍了研究区主要地震反射波与煤层对应关系.通过对T13、T12、T10、T6反射波研究,采用地震波阻抗反演技术对该区煤层厚度及煤层结构进行解释,查明井田构造形态和主要煤层赋存变化特征,取得了较好地震勘探效果.     

The Pindos Fold-and-thrust belt (Greece): inversion kinematics of a passive continental margin

Continuous exposure in the Pindos mountain chain (Greece) and the detailed stratigraphic measurements in the area enable us to construct eight balanced cross sections across the Pindos Fold-and-thrust belt (PFTB) and to approach quantitatively some parameters which controlled foreland evolution. The 160-km-wide passive continental margin of the Apulian continent in Greece was progressively shortened from east to west at rates of 6 mm/year between the Early Oligocene and Late Eocene. From the rear to the frontal part of the wedge, fault-bend folds, duplexes and imbricates were formed, while strain was partitioned into faulting (~34%), layer parallel shortening (~23%) and buckling (~9%). Foreland subsidence and internal deformation of the orogenic wedge are strongly affected by two parameters of equal importance: the thrust load of the overthrusted microcontinents and the rigidity of the underthrusted Apulian passive margin. Changes in the thickness of the pre-orogenic sediments and reactivated transform faults induced salients. During the Lower Miocene, the orogenic wedge in the Peloponnese suffered additional uplift and westward gravitational gliding induced by the intracontinental subduction of the Palaeozoic rift zone of the Phyllite-Quartzite Series, which was reactivated and returned to the earths surface during the Hellenic orogeny.     

New basin modelling results from the Polish part of the Central European Basin system: implications for the Late Cretaceous–Early Paleogene structural inversion

Combined subsidence and thermal 1D modelling was performed on six well-sections located in the north-western Mid-Polish Trough/Swell in the eastern part of the Central European Basin system. The modelling allowed constraining quantitatively both the Mesozoic subsidence and the magnitude of the Late Cretaceous–Paleocene inversion and erosion. The latter most probably reached 2,400 m in the Mid-Polish Swell area. The modelled Upper Cretaceous thickness did not exceed 500 m, and probably corresponded to 200–300 m in the swell area as compared with more than 2,000 m in the adjacent non-inverted part of the basin. Such Upper Cretaceous thickness pattern implies early onset of inversion processes, probably in the Late Turonian or Coniacian. Our modelling, coupled with previous results of stratigraphic and seismic studies, demonstrates that the relatively low sedimentation rates in the inverted part of the basin during the Late Cretaceous were the net result of several discrete pulses of non-deposition and/or erosion that were progressively more pronounced towards the trough axis. The last phase of inversion started in the Late Maastrichtian and was responsible for the total amount of erosion, which removed also the reduced Upper Cretaceous deposits. According to our modelling results, a Late Cretaceous heat-flow regime which is similar to the present-day conditions (about 50 mW/m²) was responsible for the observed organic maturity of the Permian-Mesozoic rocks. This conclusion does not affect the possibility of Late Carboniferous–Permian and Late Permian–Early Triassic thermal events.     

Faulted backfold versus reactivated backthrust: the role of inherited structures during late extension in the frontal Piémont nappes east of Pelvoux (Western Alps)

In the central part of the internal Western Alps, widespread multidirectional normal faulting resulted in an orogen-scale radial extension during the Neogene. We revisit the frontal Piémont units, between Doire and Ubaye, where contrasting lithologies allow analysing the interference with the N–S trending Oligocene compressive structures. A major extensional structure is the orogen-perpendicular Chenaillet graben, whose development was guided by an E–W trending transfer fault zone between the Chaberton backfold to the north and the Rochebrune backthrust to the south. The Chaberton hinge zone was passively crosscut by planar normal faults, resulting in a E–W trending step-type structure. Within the Rochebrune nappe, E–W trending listric normal faults bound tilted blocks that slipped northward along the basal backthrust surface reactivated as an extensional detachment. Gravity-driven gliding is suggested by the general northward tilting of the structure in relation with the collapse of the Chenaillet graben. The stress tensors computed from brittle deformation analysis confirm the predominance of orogen-parallel extension in the entire frontal Piémont zone. This can be compared with the nearby Briançonnnais nappe stack where the extensional reactivation of thrust surfaces locally resulted in prominent orogen-perpendicular extension. Such a contrasting situation illustrates how the main direction of the late-Alpine extension may be regionally governed by the nature and orientation of the pre-existing structures inherited from the main collision stage.     

港口工程中非均质层状地基承载力计算方法

在深水码头、防浪堤等大型港口工程建设中常常会遇到构筑物下覆地基由非均质的层状土构成的情况。对于这类非均质层状地基,我国港口工程地基规范推荐采用Hansen 61公式、Hansen 69公式和Sokolovskii数值解法3种方法计算地基的极限承载能力。然而,工程实践中由这3种方法计算得出的地基极限承载力数值往往有较大的差别。考虑到规范的继承性,基于牛顿迭代法,针对3个理论计算公式编制了计算程序,根据长江口深水治理二期工程地质实测资料,运用该计算程序对该地区地基的极限承载能力进行了计算比较。计算结果表明,对于非均质层状地基,在相同地质条件、相同荷载条件下Sokolovskii数值解法得到的承载力数值最大,Hansen 69公式得到的承载力数值最小。     

Application of the normalized surface magnetic charge model to UXO discrimination in cases with overlapping signals

This paper presents an application of the normalized surface magnetic charge (NSMC) model to discriminate objects of interest, such as unexploded ordnance (UXO), from innocuous items in cases when UXO electromagnetic induction (EMI) responses are contaminated by signals from other objects. Over the entire EMI spectrum considered here (tens of Hertz up to several hundreds of kHz), the scattered magnetic field outside the object can be produced mathematically by equivalent magnetic charges. The amplitudes of these charges are determined from measurement data and normalized by the excitation field. The model takes into account the scatterer's heterogeneity and near- and far-field effects. For classification algorithms, the frequency spectrum of the total NSMC is proposed and investigated as a discriminant. The NSMC is combined with the differential evolution (DE) algorithm in a two-step inversion procedure. To illustrate the applicability of the DE–NSMC algorithm, blind test data are processed and analyzed for cases in which signals from nearby objects frequently overlap. The method was highly successful in distinguishing UXO from accompanying clutter.     

Scattering pattern analysis and true-amplitude generalized Radon transform migration for acoustic transversely isotropic media with a vertical axis of symmetry

In multi-parameter ray-based anisotropic migration/inversion, it is essential that we have an understanding of the scattering mechanism corresponding to parameter perturbations. Because the complex nonlinearity in the anisotropic inversion problem is intractable, the construction of true-amplitude linearized migration/inversion procedures is needed and important. By using the acoustic medium assumption for transversely isotropic media with a vertical axis of symmetry and representing the anisotropy with P-wave normal moveout velocity, Thomsen parameter δ and anelliptic parameter η, we formalize the linearized inverse scattering problem for three-dimensional pseudo-acoustic equations. Deploying the single-scattering approximation and an elliptically anisotropic background introduces a new linear integral operator that connects the discontinuous perturbation parameters with the multi-shot/multi-offset P-wave scattered data. We further apply the high-frequency asymptotic Green's function and its derivatives to the integral operator, and then the scattering pattern of each perturbation parameter can be explicitly presented. By naturally establishing a connection to generalized Radon transform, the pseudo-inverse of the integral operator can be solved by the generalized Radon transform inversion. In consideration of the structure of this pseudo-inverse operator, the computational implementation is done pointwise by shooting a fan of rays from the target imaging area towards the acquisition system. Results from two-dimensional numerical tests show amplitude-preserving images with high quality.     

3D simultaneous joint PP‐PS prestack seismic inversion at Schiehallion field,United Kingdom Continental Shelf

A workflow for simultaneous joint PP‐PS prestack inversion of data from the Schiehallion field on the United Kingdom Continental Shelf is presented and discussed. The main challenge, describing reasonable PS to PP data registration before any prestack or joint PP‐PS inversion, was overcome thanks to a two‐stage process addressing the signal envelope, then working directly on the seismic data to estimate appropriate time‐variant time‐shift volumes. We evaluated the benefits of including PS along with PP prestack seismic data in a joint inversion process to improve the estimated elastic property quality and also to enable estimation of density compared with other prestack and post‐stack inversion approaches. While the estimated acoustic impedance exhibited a similar quality independent of the inversion used (PP post‐stack, PP prestack or joint PP‐PS prestack inversion) the shear impedance estimation was noticeably improved by the joint PP‐PS prestack inversion when compared to the PP prestack inversion. Finally, the density estimated from joint PP and PS prestack data demonstrated an overall good quality, even where not well‐controlled. The main outcome of this study was that despite several data‐related limitations, inverting jointly correctly processed PP and PS data sets brought extra value for reservoir delineation as opposed to PP‐only or post‐stack inversion.     

Stability of source mechanisms inverted from P‐wave amplitude microseismic monitoring data acquired at the surface

We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots.