Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

New insights into the P- and S-wave velocity structure of the D" discontinuity beneath the Cocos plate

Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary.     

气枪震源激发模式及应用

气枪震源是一种重要的人工地震震源。气枪阵列理论的提出,使得气枪阵列设计技术日趋成熟,并能够在石油勘探和地球物理探测中得到更加广泛的运用。气枪震源在不同领域中应用时,需要不同的组合和激发模式,以适应不同的探测要求。加强主脉冲和加强气泡脉冲,是目前两种主要的气枪激发模式。通过比较研究两种气枪激发模式,讨论各种激发模式在激发时间、气枪间距、频率、分辨率等方面的差异,为气枪震源的广泛运用提供依据。     

甫沙-曲吕西地区电磁阵列剖面法研究

通过对塔里木盆地西南缘甫沙-曲吕西地区电法勘探(CEMP),以电性剖面为主,综合地震、重力、地面地质等资料,运用前陆盆地构造分析学理论,对该区构造式样及局部构造进行了研究,摸索了一套研究复杂构造的方法,弥补了地震勘探的不足.     

Broadband synthetic aperture borehole radar interferometry

Trials in mines have established that wideband VHF borehole radars (BHR), working in the 10–100-MHz band, can be used to probe the rockmass between boreholes over ranges from <5 m to as much as 150 m with submeter resolution. There is evidence that ore bodies reflect these radar signals both specularly and diffusely, much as the ground/air interface does when overflown by synthetic aperture radar (SAR). In both SAR and BHR, multiple flight lines, together with diffuse reflections admit the possibility of developing interferometric 3D images of the object. This paper examines the possibility of imaging buried objects in three dimensions by interferometrically combining broadband VHF borehole radar profiles shot in adjacent pairs of boreholes. Broadbanding in BHR has the advantage of releasing the image from 2nπ phase ambiguities, but practically, interferometric borehole radar (InBHR) needs high signal-to-noise ratios (SNR) to avoid noise capture. This means that 3D InBHR is limited to ranges in wavelengths which are less than the rock's attenuation factor Q. Interferometric methods are developed which are capable of mapping ore bodies and other structures in three dimensions. Tangent plane migration methods are developed here in order to reconstruct surfaces that lie in the near-field of sparse interferometric arrays.     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

New approach to milling circuit control—robust inverse Nyquist array design

The aim of this paper is to present an application of robust inverse Nyquist array control design method to a milling circuit in Majdanpek Copper Mine. An existent model of a milling circuit was adopted to respond to robust inverse Nyquist array methodology demand by introducing model uncertainty, and optimal robust controller of milling circuit is designed. The proposed controller is robust for reference signals. The robustness of multivariable feedback systems with designed controller is tested and compared with different, previously designed multivariable controllers. Simulation results, in time domain, demonstrate the preferable performance of the designed robust controller in setpoint tracking. Implementation of the new controller without any investments in new equipment significantly improves performance of the milling process.     

卫星摄影三线阵CCD影像的EEP法空中三角测量（三）

4 平差方案及模拟数据平差实验4.1 自由网+控制点空中三角测量由前方交会(2.1)式,后方交会(2.2)式以及外方位元素平滑制约条件(2.3)式,可以构成类似经典的光束法空中三角测量,但由于EFP像点坐标是推算出来的,所以控制点不宜直接参与平差过程.于是平差要分成自由网平差及利用控制点作三维线性变换两个步骤.另一方面,卫星摄影中起始角元素大约在±0.5°左右,对于经典空中三角测量,角元素起始近似值均可按零处理.但在EFP平差中,φ角起始值φ0对空中三角测量,像元素起始近似值均可按零处理.在EFP平差中,φ角起始值φ0对整条航线的几何状态影响很大,必须采用特殊的程序预先加以确定.我们采用不断步进φ角值,比较由前方交会及(1.2)式计算的Y视差的均方根值最小者,即作为φ的最佳起始值.平差框图见图4.