海上地震虚反射相位剔除法反褶积

在分析了虚反射的形成机制及其滤波器频谱特性的基础上，提出了一种在频率域剔除地震鬼波的新方法，该方法直接剔除的是虚反射滤波器相位谱的影响，而对于振幅谱，则是通过同时消除虚反射与地震子波两者共同的影响来拓宽地震记录的有效频谱，处理过程自动避免了鬼波滤波器的“零值点”问题，其最终结果是以获取反射生活费数为期望输出，因此，该方法在消除鬼波的同时可进一步提高地震剖而后 分辨率。实验表明，本文所提出的反褶积方法具有较好的容噪能力。且在鬼波延迟时以及虚反射界面反射系数存在一定误差的情况下仍能获得令人满意的处理结果。     

带乘性噪声系统在加性噪声相关时的最优固定域平滑及反褶积算法

最优平滑及最优反褶积在石油地震勘探、通讯工程、语音处理等应用领域都具有十分重要的意义。以往的带乘性噪声系统的平滑及反褶积大都对系统模型的噪声特性有着较强的限制条件 ,要求动态噪声及观测噪声互相独立或只能在同时刻相关。该文给出了一种在动态噪声为有色噪声及动态噪声和量测噪声在有限时间段上相关的情形下带乘性噪声系统的固定域平滑及反褶积算法 ,该算法在线性最小方差意义下是最优的。通过仿真计算 ,说明了该算法的有效性。     

Deconvolution of squared velocity waveform as applied to the study of a noncoherent short-period radiator in the earthquake source

We consider an inverse problem of determination of short-period (high-frequency) radiator in an extended earthquake source. This radiator is assumed to be noncoherent (i.e., random), it can be described by its power flux or brightness (which depends on time and location over the extended source). To decide about this radiator we try to use temporal intensity function (TIF) of a seismic waveform at a given receiver point. It is defined as (time-varying) mean elastic wave energy flux through unit area. We suggest estimating it empirically from the velocity seismogram by its squaring and smoothing. We refer to this function as observed TIF. We believe that one can represent TIF produced by an extended radiator and recorded at some receiver point in the earth as convolution of the two components: (1) ideal intensity function (ITIF) which would be recorded in the ideal nonscattering earth from the same radiator; and (2) intensity function which would be recorded in the real earth from unit point instant radiator (intensity Green's function, IGF). This representation enables us to attempt to estimate an ITIF of a large earthquake by inverse filtering or deconvolution of the observed TIF of this event, using the observed TIF of a small event (actually, fore-or aftershock) as the empirical IGF. Therefore, the effect of scattering is stripped off. Examples of the application of this procedure to real data are given. We also show that if one can determine far-field ITIF for enough rays, one can extract from them the information on space-time structure of the radiator (that is, of brightness function). We apply this theoretical approach to short-periodP-wave records of the 1978 Miyagi-oki earthquake (M=7.6). Spatial and temporal centroids of a short-period radiator are estimated.     

A blind deconvolution method for ground based telescopes and Fizeau interferometers

In the case of ground-based telescopes equipped with adaptive optics systems, the point spread function (PSF) is only poorly known or completely unknown. Moreover, an accurate modeling of the PSF is in general not available. Therefore in several imaging situations the so-called blind deconvolution methods, aiming at estimating both the scientific target and the PSF from the detected image, can be useful. A blind deconvolution problem is severely ill-posed and, in order to reduce the extremely large number of possible solutions, it is necessary to introduce sensible constraints on both the scientific target and the PSF.In a previous paper we proposed a sound mathematical approach based on a suitable inexact alternating minimization strategy for minimizing the generalized Kullback–Leibler divergence, assuring global convergence. In the framework of this method we showed that an important constraint on the PSF is the upper bound which can be derived from the knowledge of its Strehl ratio. The efficacy of the approach was demonstrated by means of numerical simulations.In this paper, besides improving the previous approach by the use of a further constraint on the unknown scientific target, we extend it to the case of multiple images of the same target obtained with different PSFs. The main application we have in mind is to Fizeau interferometry. As it is known this is a special feature of the Large Binocular Telescope (LBT). Of the two expected interferometers for LBT, one, LINC-NIRVANA, is forthcoming while the other, LBTI, is already operating and has provided the first Fizeau images, demonstrating the possibility of reaching the resolution of a 22.8 m telescope. Therefore the extension of our blind method to this imaging modality seems to be timely.The method is applied to realistic simulations of imaging both by single mirrors and Fizeau interferometers. Successes and failures of the method in the imaging of stellar fields are demonstrated in simple cases. These preliminary results look promising at least in specific situations. The IDL code of the proposed method is available on request and will be included in the forthcoming version of the Software Package AIRY (v.6.1).     

鄂尔多斯南缘地区重力变化场源特征

以鄂尔多斯南缘地区布格重力异常数据及2014—2017年相对重力观测数据为基础,采用欧拉反褶积方法,对引起重力变化的场源深度进行反演,并对空间分布规律予以解释。通过构建理论模型,反演得到最优模型参数,并对实际数据进行计算和分析。为了减弱和消除欧拉解的发散性,利用水平梯度滤波法优化反演结果。结果表明:①构造指数为1时,适合对相对重力数据进行反演,当滑动窗口选择4—8倍测点间距时,可以获得较为可靠的场源参数;②在主要断裂附近,不同时间尺度的场源呈较好的一致性特征;③引起重力变化的场源深度集中在10—30 km,与鄂尔多斯南缘地区的震源深度基本一致;④不同时间尺度的场源位置相对分散,表明引起重力变化的物质流动具有随机性和波动性。     

由斑点图重建天文图像的方法

首先介绍了斑点图形成机制以及获取方法,并借助计算机模拟了斑点图的生成。然后重点描述斑点图的图像重建方法,包括频率域和空间域重建方法。在频率域重建方法部分依次详细叙述了用于复原模信息的斑点干涉法和用于复原相位信息的由模复原相位法、Knox-Thompson法、斑点掩模法;空间域重建方法部分介绍了典型的LWH法和迭代位移叠加法;另外,盲消卷积法也为斑点图重建提供了新的思路,文章对此方法也做了介绍;特别地,还对国内在该领域的相关工作做了介绍;关于斑点图重建方法在天文学、天体物理学等领域的应用,文章做了有针对性的介绍;最后一章是斑点图图像重建方法的总结。     

The enhanced spectrum of the symbiotic variable star CH Cygni on November 2, 2006

The symbiotic variable star CH Cygni went on the deepest minimum phase in December 2006. A high-resolution spectrum of CH Cygni obtained at the Bohyunsan Optical Astronomy Observatory (BOAO) on November 2, 2006, is reported. The spectrum of CH Cygni obtained at the BOAO was calibrated on a scale of an absolute-flux density. Emission lines of H I, [O III], [N II], and [S II] were strong, with several components, and Ca II, H, and K lines had P-Cygni profiles. Removal of the spectrum of the underlying M-type star from the CH Cygni spectrum was completed to obtain the true emission profiles. The radial velocities of all emission and absorption lines were measured. Some of the observed lines were also deconvoluted using multi-Gaussian functions to find the relationship between line shapes and the corresponding sources of these lines. The results obtained are discussed in terms of an accretion disk around the hot star, and in non-spherically extended nebular regions.     

Structural features of South-Adamawa (Cameroon) inferred from magnetic anomalies: Hydrogeological implications

Structural directions of the southern part of the Adamawa plateau in Cameroon are highlighted by filters applied on the magnetic anomaly map. Determination of the magnitude maxima of (1) the horizontal gradient of the field reduced to the pole, and (2) the analytic signal, allowed a map to be produced showing various lineaments interpreted as the fault system of the studied area. In particular, the faults on the shaded relief map of the horizontal gradient magnitude suggest a slide to the left of the Foumban shear zone in the studied area. Interpretation of the observed anomaly along a profile directed SSE-NNW shows the presence of a magnetized body considered as an intrusion into the granitic basement of probably basaltic volcanic rocks, put in place during the reactivation of the Foumban shear zone. In addition, the observed correlation between the mapped structural directions, the sites of thermo-mineral springs and the hydrographical network shows that the results of this study will be helpful for further hydrogeological research in the studied area.     

海洋天然气水合物勘探地震处理的最小平方反滤波设计及其应用

地震勘探的BSR识别技术是发现海洋天然气水合物的最经济、快捷、方便、有效的方法.在地震处理识别上,对子波的精确处理是水合物地震资料处理中最关键的一个环节,采用最小平方误差准则,即利用实际输出与期望输出的误差平方和为最小的条件,来确定反滤波因子,因此又称为最小平方子波整形.在地震处理程序中引入三种期望输出,即俞氏子波、雷克子波、Buttworth子波,对子波零相位化有较好的效果.根据上面的原理,我们开发了一套最小平方反滤波地震处理软件,对天然气水合物地震勘探资料进行试处理的结果表明,该软件在提高分辨率的同时保持了较高的信噪比.