混沌控制反演系统构造及算法逻辑设计

给出了适用于非线性反演的混沌控制反演系统构造方法及求解控制矩阵的算法逻辑设计。该方法在迭代反演控制参数和迭代反演输出结果之间建立耦合关系，通过时时修改控制参数保证迭代稳定收敛到预期的解空间。在求解过程中，应用数据结构方法分析了混沌控制理论中计算控制矩阵的数据结构，采用树形结构表述出控制矩阵计算过程中下标的取值逻辑，利用可以复读栈中元素的中序遍历来遍历每一棵树，寻找各元素的组合序列，给出的数值算例说明了混沌控制反演方法及其算法的有效性。     

中-意气候变化与气候变率学术研讨会总结

1 基本情况 2006年11月28～29日，由意大利欧洲一地中海气候研究中心（国家火山-地球物理研究所（CMCC—INGV）与中国气象局国家气候中心（NCC／BCC）联合举办的“中-意气候变化与气候变率学术研讨会”在意大利博洛尼亚市举行。中国气象局副局长王守荣研究员、预测减灾司前司长章国材研究员、国家气候中心任国玉研究员等一行9人出席了本次研讨会。     

基于稀疏离散τ-p变换的叠后地震道外推

这里研究的是一种基于稀疏离散τ-p变换的地震道缺失重建方法。该方法根据叠后零偏移距剖面,在局部时窗内可以看作是一系列线性同相轴的组合,使用稀疏离散τ-p和预条件双共轭梯度算法进行地震道外推,使空间方向的缺失道得到恢复。τ-p变换由于信息不足,有限的孔径和离散等因素,导致在τ-p域的结果存在假象、不准确,而在τ-p域进行稀疏对于重建缺失的信息是非常有用和必要的。同时,还可使用迭代的预条件双共轭梯度算法进行计算。     

地震信号广义时频分析及其数值实现

基于分数阶Fourier变换的信号时频分析称为广义时频分析,将其引入地震信号处理中,是对现有基于传统变换时频分析方法的广义化拓展。这里着重数值实现,探索基于分数阶Fourier变换的时频分析,并应用了其在物探技术中新思路和比较的优势。在分数阶Fourier变换基本定义的基础上,数值实现离散分数阶Fourier变换;研究三类广义时频分析方法及其数值实现;针对地震信号,以伪Wigner分布为例,展示广义时频分析效果。经研究表明:实现的数值算法正确有效;地震信号新的广义时频分布能量更加集中,时频分析的聚焦度和信息特征更加清晰。由此可见,广义时频分析为地震信息特征分析、预测和提取提供新的视角和工具,值得进一步研究和推广。     

利用太阳光度计反演渤海湾西岸大气柱水汽总量

利用CE-318型太阳光度计936nm水汽吸收通道的太阳辐射观测值和太阳光度计在该通道的透过率与水汽量关系，采用瞬态法反演了渤海湾西岸大气柱水汽总量。结果表明：利用太阳光度计936nm通道可以反演晴空大气柱水汽总量，其局限性是要在晴空下使用；渤海湾西岸大气柱水汽含量时间分布极不均匀；不同季节晴空日的水汽含量日变化有所不同。     

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

Toeplitz方程组的近似计算

从通常所使用的求解Toeplitx方程组的Fourier变换方法出发，结合预条件共乾梯度法(Preconditioned Conjugate Gradient Method,记为PCGM)，在对Toeplitz矩阵系统中的系数矩阵作ω循环延拓后再对其进行求解。理论和实际数值计算表明，该方法化于传统的采用简单循环的普通Fourier变换方法，所得结果具有较好的精度。     

波动方程叠前深度偏移及并行计算在南黄海地区的应用

近年来发展起来的微分解型叠前深度偏移成像技术，在复杂油气勘探的实践中不断完善，已经成为复杂地质目标偏移成像必不可缺的工具．同时，随着计算机技术的飞速发展，使得并行算法能够极大地提高数值计算的效率．本文以南黄海地区为例，描述了运用二维波动方程叠前深度偏移辛几何算法及并行运算的过程并给出结果．     

实现稀疏反褶积的预条件双共轭梯度法

地震勘探稀疏反褶积计算一般要导出一个Toeplitz矩阵的线性系统，通常可以用矩阵求逆、Levison递推及共轭梯度等方法直接求解。当Toeplitz矩阵的条件数很大时，数值稳定性差，甚至无法求解。使用共轭梯度法，在矩阵的对角元素上加入规则化因子，可以改善这种情况，但不能彻底解决数值稳定性和精度问题。若求解最小二乘问题的原始问题，结果会好些。线性系统形式的细微改变，将导致不同的数值计算特性。在规则化策略基础上，可巧妙地构造稀疏反褶积的问题原型，引入预条件，采用双共轭梯度法求解，从而实现稀疏反褶积，获得较好结果。数值算例表明，预条件双共轭梯度法比直接稀疏反褶积方法收敛快、精度高。     

基于小波变换的地震相干体算法研究

基于计算相干体算法,提出了用模拟地震子波的小波函数(或高分辨导数小波函数)的小波变换,得到分频瞬时相位,再计算相干体的相干体算法1和用小波变换得到的实、虚部(相当Hilbert变换)计算相干体的相干体算法2. 在油田构造解释中,为了突出小断层特征,用分频计算出的相干体进行重构. 实际资料计算表明,相干体算法2较K. J. Marfurt相干体算法抗噪声能力强;相干体算法1较相干体算法2在实际应用中效果更明显.