Nonstationary deconvolution using maximum kurtosis optimization     

Javad Jamali  Abdorrahim Javaherian 《Geophysical Prospecting》2020,68(5):1443-1455

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.  相似文献   

基于欧拉反褶积法的无人机航磁找矿应用     

乔中坤  马国庆  于平  周文纳  张志厚  焦健  周帅  孟兆海  唐水亮 《吉林大学学报(地球科学版)》2021,51(2):552-560

无人机航磁测量以其高效、安全、低成本等优点逐渐应用到中小面积大比例尺资源勘探领域。本文针对安徽芜湖市三山区22 km²区域构造刻画和找矿前景预测,开展1∶20 000高精度无人机航磁测量工作;推导了欧拉反演求解过程,并分析对比构造指数选取,使用构造指数为0的欧拉反褶积对航磁数据进行反演求解,获得地下磁异常体构造边界解集信息;将其应用到研究区域航磁数据反演,根据反演得到的边界位置和深度信息圈定了18处断裂,并结合异常特征给出了8.7 km²岩浆岩分布区,根据宁芜区域铁矿成矿特点,圈定岩浆分布区及边界为成矿详查区。  相似文献   

金融领域的随机建模与基于软件R的Monte Carlo模拟(4):随机微分方程模型     

毛学荣  李晓月 《南京气象学院学报》2015,7(4):313-322

主要研究线性随机微分方程模型,为此定义Itô随机微积分,建立Itô公式.鉴于研究的重点是利用R软件进行数值模拟,所以详细讨论了过去10多年来随机微分方程数值解的研究.  相似文献   

应用改进欧拉算法解算磁性目标空间位置参数     

卞光浪  翟国君  刘雁春  黄谟涛  欧阳永忠 《测绘学报》2011,40(3):386-392

指出目前磁性目标探测中常规欧拉方法存在的局限性,提出在欧拉窗口内视地磁场和其他目标干扰异常联合影响为线性变化,将构造指数作为动态变化参数代入欧拉方程一并求解,并利用多元线性回归方法解决了欧拉方程非线性问题.在质量控制方案中,根据构造指数和目标深度变化规律,给出离散欧拉解滤波措施;采用层次聚类分析方法对滤波后欧拉解进行分...  相似文献   

基于SIFT特征配准的多帧迭代盲解卷积算法   总被引：1，自引：0，他引：1  

宋向  耿则勋  王洛飞  魏小峰 《测绘科学》2011,36(6):109-111

在未配准的情况下,多帧图像解卷积算法不能发挥其利用多帧图像信息互补的优势,针对这种情况,本文将基于SIFT特征配准的算法引入多帧图像解卷积算法中,这一做法克服了各帧影像之间相同位置不同细节像素之间的影响,有效地提高了复原图像的质量。实际观测数据的复原结果证明了本文思想的正确性,配准后图像的复原结果更加清晰,且体现出观测目标更多的细节信息。  相似文献   

基于单位四元数的单独像对相对定向     

柳建锋  盛庆红  王惠南 《测绘工程》2015,(6)

目前对单位四元数描述空间旋转的应用,基本上是将其转换成旋转矩阵的形式,而忽略了单位四元数描述空间向量旋转具有描述几何关系的直观性以及运算简洁的优势。基于此提出一种基于单位四元数的单独像对相对定向法,方法首先采用向量描述同名光线,然后利用单位四元数描述立体像对中左右影像同名光线的变换关系,建立单位四元数共面条件方程模型,再按照带有约束条件的间接平差原理进行最小二乘平差,确定相对定向元素。针对RC30航空影像和低空无人机影像采用单位四元数法分别进行单独像对相对定向实验,结果表明单位四元数相对定向方法计算精度较高,为解决摄影测量中坐标变换问题提供更多选择。  相似文献   

不同高度数据联合欧拉反褶积法研究          下载免费PDF全文

周文月  陈昌昕  侯振隆  秦朋波  王泰涵 《地球物理学报》2018,61(8):3400-3409

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

基于非结构有限元的时间域海洋电磁三维反演          下载免费PDF全文

惠哲剑  殷长春  刘云鹤  张博  任秀艳  王聪 《地球物理学报》2020,63(8):3167-3179

海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性.  相似文献   

Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step     

Songhao Shang 《水文研究》2012,26(22):3338-3343

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

中国地区水平地壳形变模型的建立     

朱新慧  王刃  孟键 《测绘科学》2010,35(1):13-15

我国大陆板块是欧亚大陆板块的重要组成部分,其板块运动剧烈,地质构造复杂,是全球研究板块运动的重要地区。本文通过对各个块体分别建模和解算,得出结论:中国地区除东北块体相对于欧亚板块不动外,其他块体均有不同程度的形变运动,东部形变较小,西部形变较大;华北块体总体上以大约6mm/a的速率向东运动;华南块体以约5~8mm/a的速率向东南方向运动;新疆块体以约10~15mm/a的速率向北东方向运动;西藏块体也向北东方向运动,但运动速率从南到北逐步减少。  相似文献