首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 162 毫秒
1.
微震源反演数学模型、反演算法等对象是微震源定位误差的重要来源。针对P波速度易衰减、首波到时数据拾取误差较大的情况,基于到时差模型(ATD)和P波到时方程组,建立到时差比值法反演数学模型(ATDRM),从目标函数中消除了P波速度和发震时刻等参数;同时,引入范数概念,得到改进的L1、L2、L3、L4范数模型,设计一系列反演模型对比实验和反演算法对比实验综合分析研究这些模型的优越性。工程数据显示:ATDRM模型的定位效果整体比经典ATD模型更好,ATDRM-L1模型震源误差可减少13.437 0 m;反演模型的定位误差受反演算法影响较大,单形替换法相比遗传算法(GA)、非线性最小二乘法(NLSM)、模拟退火算法(SA)等方法定位精度更高,稳定性更好;在应用单形替换法求解模型时,反演误差会随着范数参数N的增加(L1→L4)呈现先下降再上升的现象;ATDRM模型及其范式形式对提高微震源定位精度有一定价值。  相似文献   

2.
利用水平与竖向谱比(HVSR)方法反演场地速度结构是国际上迅速发展的研究领域.HVSR反演计算实质是一个土层场地模型空间搜索的全局优化问题,当模型搜索空间的复杂程度增大时,目前常用的搜索算法收敛速度慢,计算效率较低.本文实现了一种结合遗传和模拟退火方法优点的混合全局优化HVSR反演算法,通过理论模型和竖向台阵实测数据的检验,表明该算法能获得很好的反演效果,较好地解决了蒙特卡罗方法收敛速度慢,遗传算法收敛早熟和模拟退火算法搜索效率低的问题.本文在此基础上讨论了单台加速度S波记录用于场地速度结构HVSR反演的适用性,为基于单个地震台的地震观测记录反演浅层速度结构提供了一种高效且较为准确的反演方法.  相似文献   

3.
基于FFT-MA谱模拟的快速随机反演方法研究   总被引:3,自引:2,他引:1       下载免费PDF全文
虽然基于地质统计学的随机反演方法能够有效融合测井资料中的高频信息,但计算效率低,占用内存大,限制了它在实际资料中的应用.本文在保留传统随机反演方法优点的基础上,创造性地引入傅里叶滑动平均(Fast Fourier Transform-Moving Average,FFT-MA)谱模拟进行频率域的地质统计模拟,并利用逐步变形算法(Gradual Deformation Method,GDM)确保模拟结果与实际地震数据的匹配,构建了基于FFT-MA谱模拟的新的快速随机反演方法.与常规随机反演相比,新方法不仅分辨率高,而且能够使反演解得到快速收敛,有效提高计算效率,减少内存占用.模型试算获得了与理论模型吻合度较好的高分辨率反演结果.实际资料分析也表明新方法所得到的高分辨率反演结果能够对薄互储层进行良好的展示,为薄储层的识别提供高效可靠的技术支持.  相似文献   

4.
柳旭峰  许才军 《地震学报》2013,35(2):151-159
视震源时间函数的提取是研究震源参数的重要途径. 本文提出了利用改进的粒子群(PSO)算法反演视震源时间函数的方法, 以水平线方法得到的结果作为PSO算法的初值, 并对PSO算法的惯性因子和学习因子进行改进, 提高计算效率. 采用改进的PSO算法对模拟数据进行了反演计算, 并与映射Landweber反褶积(PLD)方法和遗传算法(GA)进行了对比分析. 结果表明, 相对于PLD方法, 改进的PSO算法反演结果与真实结果误差更小; 相对于遗传算法, 改进的PSO算法计算效率提高了5倍以上. 最后, 利用改进的算法对2005年10月8日巴基斯坦克什米尔MW7.6地震的P波视震源时间函数进行了提取, 结果表明此次地震P波视震源时间函数在25 s之内, 震源沿西北向破裂. 该结果与张勇等的结果一致.   相似文献   

5.
地震全波形反演理论与技术虽已得到了广泛研究,但周波跳跃等问题的存在严重制约了常规全波形反演方法的实用化进程.基于遗传算法的全波形反演方法能够在一定程度上较好地缓解常规全波形反演面临的初始模型依赖性问题,但是当前方法仍存在收敛性和巨大的计算量问题.本文提出一种混合自适应遗传算法(HAGA),并提出基于HAGA的稳健全波形反演方法,该方法将基于HAGA的反演与基于共轭梯度法的常规全波形反演交替迭代进行,其可兼顾反演计算效率与精度.数值测试结果表明,局部与全局优化交替迭代的全波形反演方法集合了局部优化反演的高效与全局优化反演的稳定的优点,大大降低了全波形反演对初始模型的依赖性,能够有效的缓解常规全波形反演的周波跳跃问题.  相似文献   

6.
双侧向电阻率测井反演算法研究   总被引:12,自引:3,他引:9       下载免费PDF全文
双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难.本文利用变形玻恩(Born)迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演.在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林(Green)函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量.通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小.  相似文献   

7.
收集了1978~2003年位于渤海湾盆地及其邻近地区(大体位于东经112°~122°,北纬34°~44°之间)所记录到的的总计10534个地震或人工爆破事件的Pn、Pg以及P波到时数据,运用GT5准则,从中挑选出相对精确和可信度较高,且具有代表性的的10个地震事件,使用Messy GA算法,分别反演了这10个地震事件所对应的一维地壳P波速度模型.反演的局部地壳模型与已知结果对应较好,基本反映了震中区的地壳速度结构特征.该方法原理简单,计算方便,可用于确定局部和区域的地壳速度模型.  相似文献   

8.
上世纪中叶以来,微震探查法等基于面波相速度频散测定以及其反演的地下构造勘探方法,在地壳研究、工程地质勘探等方面得到了广泛的应用.本研究分析了遗传算法(GA)应用在面波频散反演当中出现的问题,采用遗传算法和Levenberg-Marquardt 算法(LM法)的联合运用法——以遗传算法反演出来的大局的搜索结果为初始点,再进行LM法搜索,由此克服两者的缺点,实现高精度的反演.  相似文献   

9.
基于不等式约束的三维电阻率探测混合反演方法   总被引:2,自引:2,他引:0       下载免费PDF全文
三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.  相似文献   

10.
随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明:该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。  相似文献   

11.
一种全局优化算法:遗传算法-单纯形法   总被引:20,自引:0,他引:20  
针对遗传算法搜索效不高和早熟的缺点,本文提出了将遗传算法和单纯形法组合的一种全局优化算法,并用七个典型检验函数验证了此方法的有效性。  相似文献   

12.
通常认为,成层介质中P波和SV波互相转化、耦合交织而不可以解耦得到单独的P和SV波入射场.针对钻井台阵的两处三分量加速度记录,根据反向方位角将两水平分量分解为径向水平分量和横向水平分量,假设横向水平分量由SH波系组成,竖向分量和径向分水平量由P-SV波系组成.本文提出在频域内解联立方程组的方法,并结合遗传算法-单纯形法混合优化策略,对P-SV耦合场包括SH波系进行了联合剥离反演,最终得到了基岩斜入射P、SV和SH波.最后,通过场地条件差别较大的另一附近钻井台站的联合反演结果,初步验证了其合理性.该项研究对于断层特征和震相特性研究、近地表“软”基岩或土层介质参数研究、场地效应的细致评价具有一定意义.  相似文献   

13.
ABSTRACT

In this study, a data-driven streamflow forecasting model is developed, in which appropriate model inputs are selected using a binary genetic algorithm (GA). The process involves using a combination of a GA input selection method and two adaptive neuro-fuzzy inference systems (ANFIS): subtractive (Sub)-ANFIS and fuzzy C-means (FCM)-ANFIS. Moreover, the application of wavelet transforms coupled with these models is tested. Long-term data for the Lighvan and Ajichai basins in Iran are used to develop the models. The results indicate considerable improvements when GA selection and wavelet methods are used in both models. For example, the Nash-Sutcliffe efficiency (NSE) coefficient for Lighvan using FCM-ANFIS is 0.74. However, when GA selection is applied, the NSE is improved to 0.85. Moreover, when the wavelet method is added, the performance of new hybrid models shows noticeable enhancements. The NSE value of wavelet-FCM-ANFIS is improved to 0.97 for Lighvan basin.
Editor D. Koutsoyiannis Associate editor E. Toth  相似文献   

14.
地震叠前逆时偏移算法的CPU/GPU实施对策   总被引:16,自引:8,他引:8       下载免费PDF全文
相较于单程波偏移算法而言,逆时偏移成像方法以其物理基础为依托优势,几十年来一直备受国内外地球物理学家的青睐.目前的逆时偏移(RTM)若直接采用双程波动方程进行延拓,尽管可以回避上下行波的分离处理,然就已有算法而言,其计算量和I/O(输入/输出)量却是最大的.针对此问题,本文在分析现行逆时偏移的多种算法基础上,提出利用CPU/GPU(中央处理器/图形处理器)作为数值计算核心,建立随机边界模型,从而克服存储I/O难题和提高计算效率.在实际的数据测试中,本文的方法可以大幅度的提高计算效率和减少存储单元,从而促使其高效地应用于生产实际.  相似文献   

15.
基于牛顿迭代法和遗传算法的CSAMT近场校正   总被引:1,自引:0,他引:1       下载免费PDF全文
栾晓东  底青云  雷达 《地球物理学报》2018,61(10):4148-4159
由于可控源音频大地电磁法(CSAMT)采用人工场源在大大增加信号强度的同时也带来了在近区产生非平面波效应的问题,其表现特征为在双对数坐标中卡尼亚电阻率曲线呈45°上升,即所谓的近场效应.本文首先提出了基于牛顿迭代法求解视电阻率方程的近场效应校正方法,通过对二层、三层理论模型试算验证了该方法的有效性;为进一步增强对噪声的压制能力,本文进一步引入了遗传反演的思想,将全频段误差最小作为目标函数,采用自适应正则化方法引入模型参数梯度最小作为稳定因子,由遗传算法求解得到校正后视电阻率,通过对理论模拟数据加随机噪声的校正结果,表明该方法在实现近场校正的同时能剔除噪声干扰,得到的结果与理论模拟曲线吻合度较高.通过对内蒙古曹四夭钼矿区实测CSAMT数据处理结果表明,牛顿法和遗传算法均能明显校正CSAMT因近场效应引起的假高阻异常,能更好地反映地电结构特征.  相似文献   

16.
Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity.  相似文献   

17.
有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点.  相似文献   

18.
角度域成像道集是叠前深度偏移的重要输出结果,它是偏移速度分析、各向异性分析和AVA分析等研究工作的基础.目前存在的角度域成像道集的生成方法受计算效率或角度分辨率的影响,仍然满足不了实际生产的要求.角度域成像道集的生成方法可以大致分为直接法和间接法两大类.在直接方法中,波矢量方向计算和局部平面波分解是两个最重要的内容,它们共同决定角度域成像条件的实现效率和角度域成像道集的质量.为了完善现有的角度域成像道集生成方法,本文提出一种新的波矢量计算方法和局部平面波分解方法.本文先用波动方程任意时刻的柯西条件构造一个只含原波场负频率成分的柯西波场,然后根据柯西波场在时间波数域的振幅谱计算波场的波矢量方向.该方法仅在需要计算波矢量方向的时刻合成柯西波场,不需要增加额外的数据读写操作,是一种高效的波矢量计算方法.本文还以柯西波场为基础提出一种高效的局部平面波分解方法,保证角度域成像条件的实现效率.结合柯西波场和局部平面波分解方法,本文最后给出一种新的角度域成像方法.文中最后的数值实验证明该方法得到的角度域成像道集具有理想的角度分辨率,可以反映地下构造的角度照明情况.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号