利用改进的粒子群算法反演视震源时间函数

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

基于交替方向乘子算法的二维磁异常稀疏反演

磁异常反演是获取地下场源磁化率分布的重要手段之一,在地球勘探中扮演着重要角色.在磁异常反演中,对比光滑反演,稀疏反演的结果具有边界分明,物性参数分布集中的特点,更符合实际情况.针对稀疏反演,本文首先构建了具有代表性的基于L₁范数目标函数,利用交替方向乘子算法可分离凸函数的特点,将极小化L₁范数的优化问题分解为一系列的子问题,通过对子问题求解获得原问题的解;为了增强交替方向乘子算法的适应性,本文结合广义软阈值函数将交替方向乘子法推广于L_p（0 < p < 1）范数的反演中.为了验证本文提出的算法的有效性,采用了三种常规模型进行模拟实验.与基于L₂范数的反演算法进行实验对比,结果表明,本文算法得到了边界清晰,磁化率分布更集中的反演结果.最后,将基于交替方向乘子算法的L₁和L_p（0 < p < 1）范数的反演应用到青海省尕林格铁矿保护区获得的实际磁异常数据中,获得了较为符合实际地质情况的稀疏反演结果.

     

正则化偏移成像的全局优化快速算法

目前,偏移后的地震剖面往往只是一个地质构造图像,还不能为后续的岩性分析和油气储层属性的提取提供更精确的信息.为了得到高分辨率真振幅的图像,建议采用正则化偏移成像方法.针对本问题数据规模大和正演算子矩阵稀疏的特点,提出采用一种新的算法——无记忆拟牛顿-模拟退火法对偏移算子方程进行求解.该方法综合了无记忆拟牛顿法优良的局部搜索能力以及模拟退火法的全局最优性质.用该方法得到的全局最优解不仅成像位置正确,还能提供更准确的振幅及其他属性信息.通过理论模型试算和实际资料处理,表明在正演算子准确的情况下,该方法能减弱Kirchhoff偏移的假象,得出接近真实的反射系数分布.因此,该方法在地震成像方面是有实际应用前景的.     

位场各阶垂向导数换算的新正则化方法

位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.     

基于辛格式奇异核褶积微分算子的地震标量波场模拟

本文在对地震波场进行模拟时,采用辛差分格式对波动方程进行时间离散,采用奇异核褶积微分算子对波动方程进行空间离散.该方法尽管增加了一些计算量,但提高了计算精度和稳定性;相对于其他非辛算法,它是全局保结构的,并且具有较强的长时间跟踪能力.该方法为解决大尺度、长时程地震波场的高精度模拟问题提供了一种新的、有效的选择.     

2007年宁洱6.4级地震序列精定位

Five mobile digital seismic stations were set up by the Earthquake Administration of Yunnan Province near the epicenter of the main shock after the Ning’er M6.4 earthquake on June 3, 2007. In this paper, the aftershock sequence of the Ning’er M6.4 earthquake is relocated by using the double difference earthquake location method. The data is from the 5 mobile digital seismic stations and the permanent Simao seismic station. The results show that the length of the aftershock sequence is 40km and the width is 30km, concentrated obviously at the lateral displacement area between the Pu’er fault and the NNE-trending faults, with the majority occurring on the Pu’er fault around the main shock. The depths of aftershocks are from 2km to 12km, and the predominant distribution is in the depth of 8～10km. The mean depth is 7.9km. The seismic fault dips to the northwest revealed from the profile parallel to this aftershock sequence, which is identical to the dip of the secondary fault of the NE-trending Menglian-Mojiang fault in the earthquake area. There are more earthquakes concentrated in the northwest segment than in the southeast segment, which is perhaps related to the underground medium and faults. The depth profile of the earthquake sequence shows that the relocated earthquakes are mainly located near the Pu’er fault and the seismic faults dip to the southwest, consistent with the dip of the west branch of the Pu’er fault. In all, the fault strike revealed by earthquake relocations matches well with the strike in the focal mechanism solutions. The main shock is in the top of the aftershock sequence and the aftershocks are symmetrically distributed, showing that faulting was complete in both the NE and SW directions.     

2013年洱源M_S5.5地震序列M_L≥3.4地震的震源机制解研究

利用云南省地震台网的宽频带数据,采用CAP方法计算了2013年3月3日洱源MS5.5地震前后共10次ML≥3.4地震的震源机制解,分析了洱源MS5.5地震前后震源机制解的一致性变化特征,结合精定位结果、震后现场灾害评估报告及震区地质构造情况,判定维西—巍山断裂以西的苍山西麓炼铁盆地东缘主边界断裂带为洱源MS5.5地震的发震构造,震源机制解节面Ⅱ(走向160°、倾角48°、滑动解-99°)为其发震断层面。     

华北地区中小地震重新定位和地震活动特征研究

选取了华北地区1978—2011年289个台站记录到的24753个地震事件, 使用双差定位方法对这些地震进行了重定位, 获得了20512个地震的重定位结果, 平均定位误差在EW、 NS、 UD方向分别为0.54 km, 0.57 km和0.57 km。 相比于初始结果和绝对定位结果, 双差定位结果在平面上分布更集中, 沿断裂分布特征更加明显, 震源深度也更加合理。 通过剖面分析, 张北地震带、 唐山地震带、 邢台地震带均为高倾角断裂, 南部剖面6震源深度分布有随时间变浅的趋势。 通过与前人的研究结果进行对比, 本次研究使用的时间较长, 地震数量较大, 得到了更加可信的结果。     

CEC比值法在凝灰质砂岩储层测井评价中的应用——以海-塔盆地X凹陷凝灰质砂岩储层为例

凝灰质砂岩储层一直是测井评价中的难点.凝灰质的存在使储层的孔隙结构和物性变化都很大,其含量直接影响储层参数的求取,对储层饱和度的影响不容忽视.本文以海-塔盆地X凹陷的凝灰质砂岩储层为例进行测井评价研究.基于泥质和凝灰质的测井响应差异利用粒子群和细菌觅食的混合优化算法计算储层中两者的含量,然后通过阳离子交换量(CEC)的实验数据验证了凝灰质具有导电性,进而利用CEC与电阻率之间的关系得到凝灰质电阻率的计算方法,并且应用到饱和度的计算,最终得到了一种新的计算凝灰质砂岩储层饱和度的方法——CEC比值法,并取得了良好的应用效果.     

70 kV超低管电压低对比剂用量冠状动脉CTA研究

目的：探讨采用基于AI的新迭代重建算法与低管电压在降低肥胖患者冠状动脉CTA辐射剂量的可行性。方法：选取行冠状动脉CTA检查的肥胖患者40例,随机分为两组。常规剂量组（A组）和低管电压剂量组（B组）各20例。常规剂量组采用管电压120 kV,管电流200 mAS扫描,低剂量组采用管电压80 kV,管电流200 mAS。采集数据传输到后处理工作站后,分别采用常规迭代重建算法和基于AI的新迭代重建算法的图像重建,对比分析两组的主动脉平均CT值、主动脉及左冠状动脉的信噪比、对比噪声比和右冠状动脉的噪声、信噪比。结果：B组主动脉的平均CT值、主动脉及左冠状动脉的信噪比、对比噪声比和右冠状动脉的噪声、信噪比、对比噪声比与A组无明显统计学差异;两组冠状动脉主观图像质量评分无明显差异;但B组的辐射剂量较A组降低68.2％。结论：新迭代重建算法能够显著性提高重建的图像质量,在新迭代重建算法的加持下,用80 kVp的管电压可以明显降低辐射剂量,同时获得与常规剂量图像质量相当的冠状动脉CTA图像。