块状岩石矿物组分热红外光谱反演实验研究

针对地面研究用岩石的热红外光谱反演其矿物组分过程中,由于自然岩石内矿物组分分布不均匀,给确定观测区域内矿物组分的精确值造成困难的问题,该文采用块状矿物拼接的方法模拟自然岩石,使模拟岩石的矿物组分可以精确获取。进一步,针对模拟岩石制订旋转试样法观测方案,从理论上研究该方案存在的观测误差,并通过光谱解混验证其观测精度。     

Geostatistical approaches to refinement of digital elevation data

Data refinement refers to the processes by which a dataset’s resolution, in particular, the spatial one, is refined, and is thus synonymous to spatial downscaling. Spatial resolution indicates measurement scale and can be seen as an index for regular data support. As a type of change of scale, data refinement is useful for many scenarios where spatial scales of existing data, desired analyses, or specific applications need to be made commensurate and refined. As spatial data are related to certain data support, they can be conceived of as support-specific realizations of random fields, suggesting that multivariate geostatistics should be explored for refining datasets from their coarser-resolution versions to the finer-resolution ones. In this paper, geostatistical methods for downscaling are described, and were implemented using GTOPO30 data and sampled Shuttle Radar Topography Mission data at a site in northwest China, with the latter’s majority grid cells used as surrogate reference data. It was found that proper structural modeling is important for achieving increased accuracy in data refinement; here, structural modeling can be done through proper decomposition of elevation fields into trends and residuals and thereafter. It was confirmed that effects of semantic differences on data refinement can be reduced through properly estimating and incorporating biases in local means.     

Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

A determination of source properties of large intraplate earthquakes in Alaska

Historically, large and potentially hazardous earthquakes have occurred within the interior of Alaska. However, most have not been adequately studied using modern methods of waveform modeling. The 22 July 1937, 16 October 1947, and 7 April 1958 earthquakes are three of the largest events known to have occurred within central Alaska (M _s =7.3,M _s =7.2 andM _s =7.3, respectively). We analyzed teleseismic body waves to gain information about the focal parameters of these events. In order to deconvolve the source time functions from teleseismic records, we first attempted to improve upon the published focal mechanisms for each event. Synthetic seismograms were computed for different source parameters, using the reflectivity method. A search was completed which compared the hand-digitized data with a suite of synthetic traces covering the complete parameter space of strike, dip, and slip direction. In this way, the focal mechanism showing the maximum correlation between the observed and calculated traces was found. Source time functions, i.e., the moment release as a function of time, were then deconvolved from teleseismic records for the three historical earthquakes, using the focal mechanisms which best fit the data. From these deconvolutions, we also recovered the depth of the events and their seismic moments. The earthquakes were all found to have a shallow foci, with depths of less than 10 km.The 1937 earthquake occurred within a northeast-southwest band of seismicity termed the Salcha seismic zone (SSZ). We confirm the previously published focal mechanism, indicating strike-slip faulting, with one focal plane parallel to the SSZ which was interpreted as the fault plane. Assuming a unilateral fault model and a reasonable rupture velocity of between 2 and 3 km/s, the 21 second rupture duration for this event indicates that all of the 65 km long SSZ may have ruptured during this event. The 1947 event, located to the south of the northwest-southeast trending Fairbanks seismic zone, was found to have a duration of about 11 seconds, thus indicating a rupture length of up to 30 km. The rupture duration of the 1958 earthquake, which occurred near the town of Huslia, approximately 400 km ENE of Fairbanks, was found to be about 9 seconds. This gives a rupture length consistent with the observed damage, an area of 16 km by 64 km.     

基于方向tilt-Euler的三维磁数据快速反演

三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组（方向tilt-Euler法）,用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.

     

不同高度数据联合欧拉反褶积法研究

基于不同测量高度重力场及其梯度数据可同时对应同一场源并用于反演场源位置的分析原理,拓展不同高度场数据在欧拉反褶积法中的应用范围.首先,立足于对欧拉反褶积方法的理论研究基础,提出不同高度数据融合联合欧拉反演公式.其次,在理论模型上对多种高度数据联合反演做了测试分析计算,验证了不同高度场数据融合联合欧拉反褶积法能够改善位场解释中单一观测面数据计算带来的解的发散问题,收敛过程由此改善.最后,将本文方法应用于龙门山地区实际重力数据的解释,获得了研究区断裂分布特征.

     

A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.     

Nonstationary deconvolution using maximum kurtosis optimization

Deconvolution is an essential step for high-resolution imaging in seismic data processing. The frequency and phase of the seismic wavelet change through time during wave propagation as a consequence of seismic absorption. Therefore, wavelet estimation is the most vital step of deconvolution, which plays the main role in seismic processing and inversion. Gabor deconvolution is an effective method to eliminate attenuation effects. Since Gabor transform does not prepare the information about the phase, minimum-phase assumption is usually supposed to estimate the phase of the wavelet. This manner does not return the optimum response where the source wavelet would be dominantly a mixed phase. We used the kurtosis maximization algorithm to estimate the phase of the wavelet. First, we removed the attenuation effect in the Gabor domain and computed the amplitude spectrum of the source wavelet; then, we rotated the seismic trace with a constant phase to reach the maximum kurtosis. This procedure was repeated in moving windows to obtain the time-varying phase changes. After that, the propagating wavelet was generated to solve the inversion problem of the convolutional model. We showed that the assumption of minimum phase does not reflect a suitable response in the case of mixed-phase wavelets. Application of this algorithm on synthetic and real data shows that subtle reflectivity information could be recovered and vertical seismic resolution is significantly improved.     

塔里木盆地航磁场分析与磁源体结构

本文研究塔里木盆地区域磁异常图的反演及磁源体结构.由于众多异常的叠加和反演固有的多解性,区域磁异常图的准确解释是非常困难的.三维欧拉反褶积是一种确定地质体位置和埋藏深度的自动定量反演方法,比较适用于计算区域磁异常源的埋藏深度.由于大型克拉通沉积盆地地层具有上新下老的规律性,将磁异常源分解为三个深度层次,圈定它们各自的分布区域,便可将它们与形成的地质作用及时代联系起来,为准确解释区域磁异常图提供可靠的依据.本文应用三维欧拉反褶积反演方法,计算出的塔里木盆地深度为2~5 km、5~10 km、10~20 km三个等级的磁异常源,它们与形成的地质作用及时代分别为: 中生代构造运动,海西期玄武岩侵位和太古代结晶基底的变质作用;圈定了它们各自的分布区域.     

柯西约束盲反褶积技术在井间地震的应用

利用地面三维地震资料将二维井间地震资料推广到井间范围以外。对于提高井间地震的效益,加快其应用十分重要．为了将井间地震二维资料推广到井外三维,需要从地面地震低频信息提取层位、断面的几何信息,反褶积方法是重要的部分．本文给出了盲源反褶积方法的一种具体实现。并结合优化的预条件共轭梯度法以改善算法的稳定性,同时减少计算量．然后对经过高频恢复的地面地震数据与井间地震数据进行联合约束反演,有效地提高了地面地震的频带．并用实际资料的处理给予证实．