Linear relationships in geomagnetic variation studies

The linear properties of Maxwell's equations enable the primary and secondary fields of geomagnetic variations to be related together by an induction tensor. Variations of total field which are confined on or near a plane, as described first by Parkinson, can occur formally only under certain conditions on either the induction tensor (that it have a simple degeneracy), or on the primary variation field components (that they consistently obey the same linear relationship over all time). The second condition is implausible for typical transient variations observed on earth, and examination of the first condition demonstrates that for variations with no consistent correlation of the horizontal components, the existence of a well-determined Parkinson vector means that observed vertical variations are secondary fields, induced by observed horizontal variations. This conclusion has important implications for the theoretical explanation of the geomagnetic “coast effect”.That primary fields be spatially uniform is a tacit condition of much theoretical work on geomagnetic variations, though the condition can in fact be relaxed to require spatial uniformity of the primary fields just within some “induction region”. If however the primary fields are not uniform over the induction region a poor Parkinson-vector determination may result. This may be most noticeable at sites which have their induction regions greatly extended by a process such as current channelling.     

HORIZONTAL STACKING AND MULTICHANNEL FILTERING APPLIED TO COMMON DEPTH POINT SEISMIC DATA*

The common depth point method of shooting in oil exploration provides a series of seismic traces which yield information about the substrata layers at one location. After normal moveout and static corrections have been applied, the traces are combined by horizontal stacking, or linear multichannel filtering, into a single record in which the primary reflections have been enhanced relative to the multiple reflections and random noise. The criterion used in optimum horizontal stacking is to maximize the signal to noise power ratio, where signal refers to the primary reflection sequence and noise includes the multiple reflections. It is shown when this criterion is equivalent to minimizing the mean square difference between the desired signal (primary reflection sequence) and the weighted horizontally stacked traces. If the seismic traces are combined by multichannel linear filtering, the primary reflection sequence will have undergone some phase and frequency distortion on the resulting record. The signal to noise power ratio then becomes less meaningful a criterion for designing the optimum linear multichannel filter, and the mean square criterion is adopted. In general, however, since more a priori information about the seismic traces is required to design the optimum linear multichannel filter than required for the optimum set of weights of the horizontal stacking process, the former will be an improvement over the latter. It becomes evident that optimum horizontal stacking is a restricted form of linear multichannel filtering.     

HORIZONTAL STACKING AND MULTICHANNEL FILTERING APPLIED TO COMMON DEPTH POINT SEISMIC DATA*

The common depth point method of shooting in oil exploration provides a series of seismic traces which yield information about the substrata layers at one location. After normal moveout and static corrections have been applied, the traces are combined by horizontal stacking, or linear multichannel filtering, into a single record in which the primary reflections have been enhanced relative to the multiple reflections and random noise. The criterion used in optimum horizontal stacking is to maximize the signal to noise power ratio, where signal refers to the primary reflection sequence and noise includes the multiple reflections. It is shown when this criterion is equivalent to minimizing the mean square difference between the desired signal (primary reflection sequence) and the weighted horizontally stacked traces. If the seismic traces are combined by multichannel linear filtering, the primary reflection sequence will have undergone some phase and frequency distortion on the resulting record. The signal to noise power ratio then becomes less meaningful a criterion for designing the optimum linear multichannel filter, and the mean square criterion is adopted. In general, however, since more a priori information about the seismic traces is required to design the optimum linear multichannel filter than required for the optimum set of weights of the horizontal stacking process, the former will be an improvement over the latter. It becomes evident that optimum horizontal stacking is a restricted form of linear multichannel filtering.     

Decomposing InSAR LOS displacement into co-seismic dislocation with a linear interpolation model： A case study of the Kunlun Mountain Ms=8.1 earthquake

It has always been a difficult problem to extract horizontal and vertical displacement components from the InSAR LOS （Line of Sight） displacement since the advent of monitoring ground surface deformation with InSAR technique. Having tried to fit the firsthand field investigation data with a least squares model and obtained a preliminary result, this paper, based on the previous field data and the InSAR data, presents a linear cubic interpolation model which well fits the feature of earthquake fracture zone. This model inherits the precision of investigation data; moreover make use of some advantages of the InSAR technique, such as quasi-real time observation, continuous recording and all-weather measurement. Accordingly, by means of the model this paper presents a method to decompose the InSAR slant range co-seismic displacement （i.e. LOS change） into horizontal and vertical displacement components. Approaching the real motion step by step, finally a serial of curves representing the co-seismic horizontal and vertical displacement component along the main earthquake fracture zone are approximately obtained.     

Interpretation of VLF Resistivity Data for Ground Water Contamination Surveys

Very low frequency (VLF) military communications systems provide a primary field that can be used for shallow geophysical surveys to locate ground water contamination and vertical geologic contacts. Useful properties that can be easily obtained from the interaction of the earth and the primary field are the magnitude of the vertical secondary magnetic field, the surface impedence, and the phase angle between the electrical and magnetic horizontal components. The variations in the secondary magnetic field can be related to vertical geologic contacts, such as the edges of landfill trenches. The surface impedence yields an apparent terrain conductivity, which can be used to locate low-resistivity anomalies often associated with contaminated ground water. The phase angle gives information on vertical variations in resistivity, phase angles less than 45° indicating increasing resistivity with depth. The depth of penetration of the VLF field is about one skin depth. For a frequency of 20 kHz, the skin depth in meters is approximately equal to 3.67 where p is terrain resistivity in ohmmeters.     

层状介质中水平磁偶极子电磁场空间分布特征

针对岩层中介电常数和磁导率随深度变化的情况,给出水平磁偶极子源电磁场的水平或垂直空间分布特征及相应的变化规律。采用Kong给出的汉克尔（Hankel）J0变换线性滤波器（241点）和汉克尔J1变换线性滤波器（241点）算法,选用均匀大地和两层地电模型,讨论介质的介电常数和磁导率变化时的水平磁偶极子电磁场空间分布特征。计算结果表明,磁场的峨分量幅值与磁导率反相关;磁场H2和电场Ey分量随着磁导率的增大而增大;低频电磁测深中介电常数变化时对电磁场基本没有影响。     

A generalized derivative operator for potential field data‡

The enhancement of potential field data using filters based on horizontal and vertical derivatives is common. As well as the direct use of the gradients themselves they are used in filters such as sunshading, total horizontal derivative, analytic signal, horizontal and vertical tilt angles, the Theta map and other filters. These techniques are high‐pass filters of different types and so enhance noise as well as detail in the data. A new derivative operator is introduced in this paper, which generalizes the effects of some of the previously mentioned filters. This filter is a linear combination of the horizontal and vertical field derivatives, normalized by the analytic signal amplitude. The filter is demonstrated on aeromagnetic and gravity data from South Africa.     

瓜什则盆地的地应力场分析与模拟

基于水压致裂法地应力测试的原理,系统性研究分析了瓜什则盆地区域内的地应力场特征,这为成宁铁路控制性工程甘青特长隧道的选线设计提供了重要的参考依据。首先,根据区域内钻孔的地应力测试结果,采用有限元方法进行了区域地应力场的三维数值模拟,确定了该区域内的水平应力大于垂直应力的基本特征;其次,根据多个水平构造应力场的水平切片,分析了水平构造应力场随深度的变化,说明了该区域内水平构造应力场具有由东北向西南挤压的应力变化特征;然后,根据区域岩体抗拉强度分布,发现该区域内存在一个低抗拉强度带,它的地层结构相对较复杂,且其工程地质特征主要表现为地层破碎带或断层;最后,采用最小二乘法拟合了CK方案隧道横截面最大切向应力与最大水平主应力的线性回归模型,并预测了该方案隧道地层的岩爆发育程度。     

FREQUENCY ELECTROMAGNETIC SOUNDING USING A VERTICAL MAGNETIC DIPOLE*

A horizontal transmitter loop (vertical magnetic dipole) is used for frequency electromagnetic (FEM) soundings. The frequency ranges from approximately 6 Hz to about 4000 Hz. The vertical and radial magnetic field components are measured for 20 frequencies per decade several hundred meters from the transmitter loop. A very small bandwidth is selected for amplification using a reference signal. An Apple computer is used for data acquisition. A computer program for Marquardt inversion optimizes the parameters for the n-layer case: the resistivities and thicknesses of individual beds and a correction factor for the primary magnetic field. Interpretation of each component individually yields practically the same parameters. Examples from the field are given with interpretation; comparison with dc resistivity measurements is provided. The ratio of vertical/radial magnetic field components vs. frequency can be transformed simply into apparent resistivity vs. apparent depth. This can be done in the field to obtain an estimation of the depth of the layer boundaries. FEM results are compared with Schlumberger d.c. sounding obtained at the same site.     

Full-correlation analysis of turbulent scattering layers in the mesosphere observed by the MU radar

We have applied a full-correlation analysis technique to the echo power fluctuations observed by the MU radar (35°N, 136°E), and analyzed the horizontal structure of the scattering pattern in the mesosphere as well as their horizontal motions. The velocity of the scattering pattern did not agree with the background wind velocity, but was associated with the horizontal propagating direction of a saturated inertia gravity wave identified in the wind field. The length of the long axis of the characteristic ellipse of the scattering pattern was approximately 50 km, and the direction was almost perpendicular to the propagating direction of the wave. The correlation time of the scattering pattern was approximately 700 s, which is much longer than the lifetime of the isolated turbulence itself. This implies that the observed scattering pattern is associated with a region where the saturated inertia gravity wave generates turbulence.     

考虑航磁水平梯度变化的ΔT网格化方法研究

实测航磁横向水平梯度反映垂直于测线的磁场梯度,比传统航磁ΔT数据包含测线之间更多的磁场信息.针对航磁数据网格化问题,采用Hardwick提出的方法,利用航磁水平梯度与ΔT数据构建拟测线,并结合Akima插值法,开展了双方向测线型ΔT网格化方法研究,最终实现了考虑航磁水平梯度变化的ΔT网格化;针对网格化结果中的虚假异常采取了有效滤波方法.通过理论模型数据和实际数据网格化处理,表明该方法可以突出航磁测线之间的异常细节、更清晰地反映线性构造或磁性体走向,提高了网格化的精度和分辨率.     

Detecting co-seismic displacements in glaciated regions: An example from the great November 2002 Denali earthquake using SPOT horizontal offsets

We use SPOT image pairs to determine horizontal offsets associated with the Mw 7.9 November 2002 Denali earthquake in the vicinity of Slate Creek, AK. Field measurements and aerial photographs are used to further characterize the geometry of the surface rupture. Aerial photographs show that shear localization occurs where the rupture trace is linear, and distributed off-fault deformation is common at fault bends and step-overs, or at geologic contacts between rock, glacial sediments, and ice. The displacement field is generated using a sub-pixel cross correlation technique between SPOT images taken before and after the earthquake. We identify the effects of glacier motion in order to isolate the tectonic displacements associated with the Denali earthquake. The resulting horizontal displacement field shows an along-strike variation in dextral shear, with a maximum of approximately 7.5 m in the east near 144° 52′W, which decreases to about 5 m to the west near 145° 45′W. If the November 2002 earthquake represents the long-term behavior of the Denali fault, it implies a westward decrease in the long-term dextral slip rate. A possible mechanism to accommodate the westward decreasing slip on the Denali fault is to transfer fault slip to adjacent east-trending contractional structures in the western region of the central Alaskan Range.     

A direct method of interpreting gravity and magnetic anomalies: The case of a horizontal cylinder

Summary A new method of interpreting the gravity and magnetic anomalies is introduced with special reference to the magnetic anomalies of a horizontal cylinder. The method consists of calculating the functions of the anomaly and its distance from an arbitrary point. These form a simple linear equation with coefficients related to the parameters defining the body. Since each observation forms a separate linear equation, the required normal equations are formed by the method of least squares and solved for the coefficients and hence for the various parameters defining the target. The discussion here is confined to the vertical magnetic anomalies. The application of the method to horizontal and total field anomalies of two dimensional bodies is also outlined.     

Influence of water saturation on horizontal and vertical motion at a porous soil interface induced by incident SV wave

Recent field observations have indicated that water saturation of soils may strongly affect the vertical ground motion. A study is therefore carried out to investigate the effect of saturation on horizontal and vertical motion at an interface of porous soils with potential contributions directed to site evaluation based on field observations of both the horizontal and vertical motion. The problem described in this paper corresponds to an SV wave incident at the interface between the overlying soil and the underlying rock formation. The soils are modeled as partially water-saturated porous material with a small amount of air inclusions, while the rock are approximately regarded as ordinary one-phase solid. Theoretical formulation is developed for the computation of amplitudes of horizontal and vertical interface motion, which are expressed as functions of the degree of saturation, the angle of incidence as well as the frequency. Numerical results are given for a typical sand to illustrate the influence of saturation on the interface motion in two directions and their ratios. The present study demonstrates that the effect of water saturation may be substantial on both the horizontal and vertical motion as well as on their ratios, implying the importance of such effects in the interpretation of field observations.     

Estimating magnetic dike parameters using a non‐linear constrained inversion technique: an example from the Särna area,west central Sweden

In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated g_zz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of g_zz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.     

Characteristics of Late Quaternary Activity of the Gadê Segment in the Madoi-Gadê Fault Zone

Madoi-Gadê fault is an active fault in the Bayan Har block.According to field investigation,there is an earthquake surface rupture fairly well preserved on the Gadê segment of the Madoi-Gadê fault zone.The length of the rupture is approximately 50km,with a general strike of NW.The maximum horizontal sinistral displacement is about 7.6m and the maximum vertical displacement is about 4m.A large number of earthquake traces are to be found along the rupture zone,and the phenomena on the surface rupture are also...     

用滑动方向拟合法反演富蕴地震断裂带应力场

本文根据野外考察测得的1931年新疆富蕴8级地震形成的地震断裂带内42组断层水平和垂直位移数据,求得42个断层面上的滑动方向.利用滑动方向拟合法(使计算的断层面上的剪切力方向与观测的滑动方向之间夹角最小)求出了断裂带上由北向南四个区段的驱动断层滑动的平均应力场.结果表明,四个区段的最大和最小主压应力轴皆水平,最大主压应力轴方位大致在 N15E 至 N30E 之间.整个断裂带的应力场相当均一,这说明富蕴大地震的发震原因不象是局部因素,而可能是大范围构造运动的结果.      

VLF RESISTIVITY MAPPING AND VERTICALIZATION OF THE ELECTRIC FIELD1

For over 20 years, powerful VLF transmitters have been used as electromagnetic sources for subsurface investigations in mining exploration. Measurements initially concerned the vertical component of the magnetic field or the inclination of the field and were later extended to measurement of the horizontal electric field in the direction of the transmitter, to determine the resistivity of the terrain. Measurement of the electric field is usually performed with electric lines, grounded or not, with lengths of at least 5 m. This paper presents the concept of a VLF resistivity meter with a very short electric sensor (1 m) and the results obtained with it. This technique improves the measurement of the electric field, which is in principle a point value. It also permits a higher spatial sampling rate and, by closely linking the electric sensor with the magnetic sensor on a lightweight mount, makes it possible for the instrument to be used by a single operator. In addition, transformation of the electric field data, analogous to reduction to the pole in magnetism, is proposed to correct the horizontal deformation of the anomalies created by polarization of the primary field. Comparison with direct current electrical measurements shows highly satisfactory correlations. This transformation, validated for VLF, can be extended to any electrical or electromagnetic method using a uniform primary field, i.e. gradient array in direct current or low-frequency magnetotellurics. We call this verticalization of the electric field. Resistivity measurements and mapping using the VLF frequency range can be applied not only to mining but also to a wide range of shallow geophysical studies (hydrology, civil engineering, etc.) and are not limited to problems concerning the location of conductive targets     

基于实测值分段线性内插模型的InSAR视线向同震位错分解-以青藏高原昆仑山MemＳsub8.1地震为例x,auto,auto,415px);}stylehf=httprztl.com business week launch

利用现场GPS定位的实测值，在尝试过最小二乘拟合并取得初步成果的基础上，进而采用更符合破裂带形态的线性立方插值函数，在破裂带主断面上，建立起一种理论与实测相结合的InSAR视线向(LOS)变化量的分解方程，得到了InSAR视线向位移分解的具有唯一性的解析解. 本文的解析法吸纳了GPS定点实测值的精度优势，利用InSAR全天候、准实时获取连续形变场的技术特点，通过数学近似，最终获得了主破裂带上连续变化的水平位错及垂直位错同震形变曲线.      

Site-Specific Design for Dual Phase Recovery and Stabilization of Pooled DNAPL

Volume reduction and lowering of capillary pressure within a large DNAPL pool are utilized as objectives in the design of a large-scale dual phase recovery system at a chemical manufacturing facility in the United States. By reducing DNAPL pool height through mass removal, capillary pressure is lowered, resulting in a reduced potential for future vertical and horizontal mobilization of the chlorinated solvent DNAPL pool. The DNAPL pool extends over an approximately 200 m by 275 m area in low permeability fill deposits overlying a clay aquitard. A three-dimensional multiphase flow model was employed to arrive at a final design incorporating nine horizontal drains (total length 664 m) and a pulsed pumping system. The numerical model was calibrated to the results of a 42-day field pilot-test involving the removal of approximately 25,000 L of DNAPL from a single, 55 m long horizontal drain. Numerical simulation revealed that gravity drainage, as opposed to hydraulic gradients in the water phase, is the dominant recovery mechanism at this site. This stems from the relatively high density and the viscosity of the DNAPL, and the relatively low permeability of the formation deposits. The use of pulsed pumping is shown to reduce the volume of contaminated ground water recovered from the 9-drain system, without significant reduction of the total volume of DNAPL recovered.