黏滞声波地震波场模拟的SG-ONADM混合方法

交错网格方法（SG）和最优近似解析离散化方法（ONADM）是两类典型的地震波场数值模拟方法.两类方法各有其优势,相对于ONADM方法,SG方法在单个时间层内的计算更为简单;相对于SG方法,ONADM方法可以在较大空间步长条件下有效压制数值频散.结合两种方法的优势,本文提出了一种新的地震波场模拟方法（SG-ONADM）.该方法对控制方程中的一阶偏导数采用SG方法给出的一阶偏导数近似公式,高阶偏导数采用ONADM方法给出的高阶偏导数逼近公式.理论分析及数值算例表明,SG-ONADM方法保留了两种方法的优势,不仅能在较大空间步长条件下有效压制数值频散,同时具有较低的内存需求量;在对同一计算区域进行波场模拟时,SG-ONADM方法的计算效率要高于SG方法和ONADM方法.最后,我们使用SG-ONADM方法进行黏滞声波波场模拟,研究了黏滞声波在复杂介质中的传播.     

航磁异常地面快速查证方法组合研究

对于地形切割厉害、环境条件恶劣地区,采用传统勘查方法开展地面查证具有较大困难.结合工程应用实践,充分发挥放射性勘查方法快速分析、对地形要求较低的优势,代替部分传统物化探方法,提出了一种放射性物探方法与传统物探方法相结合的快速查证方法组合.本方法组合分为三个阶段完成地面查证,采用层层递进工作方式,前一阶段为后一阶段工作指明方向.各阶段采用不同的放射性勘查方法和传统物探方法,且勘查方法相互补充,弥补单一勘查方法的不足,完成地面查证工任务.在工程验证阶段,采用浅钻结合X荧光岩芯分析勘查技术,较大的降低了野外工作强度.实践证明,该方法组合是一种快速高效、经济的勘探方法组合,旨在为找矿勘探工作提供一种新的思路.     

改进的高精度拟谱微分方法及其在三维非均匀介质地震传播研究中的应用

改进了用于模拟地震波场的傅里叶拟谱微分方法,改进后的方法精度是常规拟谱方法的4倍,称为改进的傅里叶拟谱方法.在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用该方法和常规拟谱方法对Marmousi模型进行数值求解,结果表明,该方法的数值频散效应明显比常规拟谱方法弱.将该方法与有限元方法在各向异性介质中进行模拟比较,发现该方法的精度接近有限元方法,数值频散效应比有限元方法明显减小,而且可在较大空间网格间距下进行计算,从而提高计算效率.在3-D非均匀介质中的地震波传播数值模拟结果表明,该方法是一种研究复杂非均匀介质中地震波传播问题的高效方法.     

测井储层分类评价方法研究进展综述

测井储层分类方法多样,每种方法的原理、计算步骤及所需资料各不相同,且各方法的适用性及应用效果均存在很大差异,本文通过文献调研,将测井储层分类方法归纳为四大类:(1)基于交会图版法的半定量储层分类方法,此方法操作简单、应用范围广但适用性不强;(2)基于流动单元概念的测井储层分类方法,此方法基于岩心物性数据可以迅速达到储层分类的目的 ,具有明确的地质意义,但结果依赖于取心数据;(3)基于多元统计法及机器学习算法的测井储层分类方法,此类方法可以有效地避免人为因素的干扰,速度快、方法多样,可以实现储层的定量分类评价,是未来发展的趋势,但分类结果意义不明确;(4)基于测井新技术新方法的储层分类方法,该方法携带了大量的地质信息,与其他分类方法结合可以更有效、更准确地评价储层.最后比较了四种分类方法的优缺点并给出了相应的选择建议,该研究对测井储层分类方法的优选具有一定的参考意义.     

基于Born波路径的高斯束初至波波形反演

为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低.     

改进的特征量方法及其在震源动力学模拟中的应用

基于曲线网格有限差分方法,针对震源动力学模型中的边界条件的问题,本文提出了一种改进的特征量方法,即将牵引力镜像方法和特征量方法相结合,在保守形式的方程中通过镜像操作完成单边导数的计算,然后有约束地调整边界处的速度和应力值,以此保证满足模型中的边界条件.改进的特征量方法可用于处理曲线网格当中一般的边界条件,从而也可用于构建曲线网格有限差分法的分裂节点模型.我们将改进的特征量方法分别用于自由表面和断层面,并分别与前人计算结果进行了比较.计算结果表明,在相同的网格划分的情况下,改进的特征量方法能够得到与其他方法相似的计算结果,但改进方法的异常震荡更小,从而证明了本文提出的改进的特征量方法用于曲线网格有限差分法分裂节点模型的可行性.     

波动方程参数反演方法研究进展

根据波场信息重建波动方程系数项是地球物理反演中一个重要研究方向.波动方程参数反演方法按照求解方式的不同,可以分为直接反演方法、迭代反演方法、搜索类反演方法、混合优化反演方法和联合反演方法.目前发展的参数反演方法多数都处于理论计算阶段.研究能求解实际地震资料的参数反演方法具有重要意义.为取得可靠的反演结果,信息的转化开发...     

近震源破裂过程反演研究——Ⅰ. 方法和数字试验

利用非线性规划研究的最新成果, 设计了一种全新的震源破裂过程的反演方法, 同时反演震源破裂的时间图像和滑动在破裂面上的分布. 该方法与目前震源过程破裂反演计算所用的两类主要方法(线性矩阵方法和全局搜索寻优方法)相比, 在线性矩阵方法对初始模型依赖大、解不稳定和全局搜索寻优方法在高维解空间寻优效率不高、易出现伪解等方面有较显著改善. 数字试验的测试表明, 本方法计算效率高, 反演结果稳定、可靠, 可有效应用于实际地震震源过程的反演研究中.     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

基于八阶NAD算子的保辛分部Runge-Kutta方法及其波场模拟(英文)

基于声波方程扩充的哈密尔顿系统,本文给出了空间精度为八阶的近似解析离散化(NAD)保辛分部Runge-Kutta方法,简称八阶NSPRK方法。该方法采用八阶精度的近似解析离散算子近似空间高阶偏微分算子,并使用二阶精度的辛分部Runge-Kutta方法进行时间离散。我们从理论和数值计算两个方面研究了八阶NSPRK方法的稳定性条件和数值频散关系,并同四阶NSPRK方法、八阶Lax-Wendroff(LWC)方法和八阶交错网格(SG)方法进行了比较。结果表明八阶NSPRK方法压制数值频散的能力显著优于传统数值计算方法。与四阶NSPRK方法和传统四阶辛格式(SPRK)方法相比,八阶NSPRK方法具有最小的数值误差和最高的计算效率:在达到同样消除数值频散的前提下,八阶NSPRK方法的计算速度约为四阶NSPRK方法的2.5倍、为四阶SPRK方法的3.4倍;八阶NSPRK方法的存储量仅为四阶NSPRK方法的47.17%、为四阶SPRK方法的49.41%。在双层介质、非均匀介质和Marmousi等复杂速度模型中,八阶NSPRK方法模拟得到的波场快照非常清晰,无可见数值频散。这些结果表明,八阶NSPRK方法在粗网格条件下能有效地压制数值频散,从而能够极大地节省计算内存,提高计算速度。总体而言,八阶NSPRK方法是一种在地震探测领域和地震学研究中有着巨大应用潜力的数值计算方法。     

A new method for aquifer system identification and parameter estimation

The standard practice for assessing aquifer parameters is to match groundwater drawdown data obtained during pumping tests against theoretical well function curves specific to the aquifer system being tested. The shape of the curve derived from the logarithmic time derivative of the drawdown data is also very frequently used as a diagnostic tool to identify the aquifer system in which the pumping test is being conducted. The present study investigates the incremental area method (IAM) to serve as an alternative diagnostic tool for the aquifer system identification as well as a supplement to the aquifer parameter estimation procedure. The IAM based diagnostic curves for ideal confined, leaky, bounded and unconfined aquifers have been derived as part of this study, and individual features of the plots have been identified. These features were noted to be unique to each aquifer setting, which could be used for rapid evaluation of the aquifer system. The effectiveness of the IAM methodology was investigated by analyzing field data for various aquifer settings including leaky, unconfined, bounded and heterogeneous conditions. The results showed that the proposed approach is a viable method for use as a diagnostic tool to identify the aquifer system characteristics as well as to support the estimation of the hydraulic parameters obtained from standard curve matching procedures. Copyright © 2012 John Wiley & Sons, Ltd.     

人工智能体积密度法判断肺亚实性结节的浸润性研究

目的:探讨人工智能（AI）体积密度法判断肺亚实性结节（SSNs）浸润性的价值。方法:回顾性分析106例患者的108枚SSNs的CT和病理结果,将结节分为腺体前驱病变组和腺癌组。通过肺结节AI软件测量并比较两组的最大CT值、最小CT值、平均CT值、峰度、偏度、Perc.25％、Perc.50％、Perc.75％、Perc.95％、结节体积、结节平均径等CT定量参数。使用Medcalc软件得出受试者工作特征曲线（ROC）,评价诊断SSNs浸润性的敏感度、特异度、阳性预测值及阴性预测值,用逻辑回归分析评估他们的诊断性能。结果:SSNs的多数CT定量参数差异存在统计学意义,其中,诊断效能最高的是Perc.25％,AUC达0.797;其次为Perc.50％ 和平均CT值,AUC均为0.787。Logistic回归分析显示,将诊断效能最高的Perc.25％ 分别与Perc.50％ 和平均CT值两两建立联合诊断模型1,其中Perc.25％ 与平均CT值的模型诊断效能最高,且联合诊断模型诊断效能高于Perc.25％ 与平均CT值单独的诊断效能。Medcalc软件分析显示,Perc.25％≥-578 HU和平均CT值≥-468 HU的SSNs病理表现为腺癌的可能性大。将Perc.25％ 与结节平均径结合,可获得对判断SSNs浸润性非常有价值的联合诊断模型2。结论:AI体积密度法对SSNs的浸润性有较高的诊断价值,联合使用Perc.25％ 与平均CT值比单独使用更能准确地判断浸润性。      

双源CT对冠脉搭桥术后桥血管评价的Meta分析

目的:利用循证医学方法,对双源CT评价冠状动脉搭桥术(CABG)后桥血管狭窄程度的价值进行Meta分析。方法:检索CNKI、万方、维普数据库以及Cochrance Library、PubMed、OVID、Embase等外文数据库,搜索双源CT诊断桥血管狭窄程度的诊断性试验。对纳入资料使用Meta-disc1.4软件进行分析计算,计算合并敏感度、特异度和诊断比值比。结果:最终入选5篇文献,其合并敏感度、特异度和诊断比值比分别为99%、99%、4984,受试者工作特征曲线下面积为99.84%。结论:双源CT对于CABG后桥血管的评价具有很高的诊断效能和临床意义。      

急性颈段脊柱外伤的MRI检查价值

目的：分析MRI在颈段脊柱急性外伤中的诊断价值，方法：回顾性分析53例急性颈部脊柱外伤患者的MRI表观特点，并与CT及X线平片比较，结果：椎体骨折41例，椎体脱位3例，脊髓损伤30例，椎间盘损伤7例，结论：MRI可以对颈椎损伤形态，及脊髓韧带椎间盘等软组织损伤作出明确的诊断，是急性颈椎损伤患者最佳的检查方法。     

探讨震兆信息及其判别依据

本文以山东省境内近十几年电感地应力测量素材为样本,提出曲线的正常、干扰、震兆信息的判别依据。并以此,对震兆信息进行判别收到了一定效果。     

3种不同肠道准备方法在CT小肠造影中的应用研究

目的:探讨3种不同肠道准备方法在CT小肠造影（CTE）中的应用价值。方法:对156例临床怀疑为小肠病变而行CTE的病例进行前瞻性研究。156例随机分为A、B和C组各52例,肠道造影剂均采用口服法。A和B两组在扫描前16h进食半流质饮食,并均以番泻叶6g导泄,扫描前12h禁食,扫描前1h起每隔15min A组尽量口服等渗甘露醇、B组同法尽量口服引用水。C组扫描前6h禁食,扫描前1h起每隔15min尽量口服饮用水。3组均在口服对比剂10min后进行CTE检查。统计比较3组小肠的充盈效果及对小肠病变的诊断价值。结果:小肠充盈效果比较A组最佳,B组次之,C组最差;对小肠病变的诊断结果显示A组对小肠各段病变的诊断价值均较显著;B组对十二指肠、空肠病变诊断价值较大,对回肠的细小病变诊断效果较差;C组对十二指肠及空肠上段病变诊断效果较好,但不利于空肠上段以远细小病变的诊断。结论:充分的肠道准备是CTE的重要技术因素,分次口服等渗甘露醇尽量充盈小肠有利于小肠病变的显示,是一种安全可靠、经济方便的肠道准备方法。      

CT能谱成像对肝脏微小转移瘤的诊断价值

目的:探讨肝脏微小转移瘤患者应用CT能谱成像诊断的效果及对患者肝功能的影响。方法:回顾性分析2015年5月至2017年5月我院收治的肝脏微小转移瘤患者80例的临床资料,根据诊断方法不同分成对照组和研究组两组。对照组常规的CT诊断,研究组应用CT能谱成像的诊断方法,对比诊断结果。结果:对比诊断结果,研究组诊断出的病灶个数以及检出率显著高于对照组（P<0.05）,同时漏检率以及误检率均低于对照组（P<0.05）;研究组应用诊断方法的诊断敏感度、阳性预测值以及特异性都高于对照组（P<0.05）;对比治疗效果,研究组的疗效显著优于对照组（P<0.05）。结论:应用CT能谱成像诊断肝脏微小转移瘤患者,有助于疾病早发现、早治疗,有更高的诊断价值,利于患者改善肝功能,值得在临床上推广应用。      

地震数字钟诊断程序设计分析

根据地震台网数据处理系统时间码出错的原因,提出并采用计算机程序诊断替代人工诊断的方法,着重分析了时间码变化规律,编写了地震数字钟诊断程序。这个诊断工具软件,对维护地震台网数据处理系统正常运行,确保地震时间准确无误,具有很强的实用性。     

多排螺旋CT冠状动脉造影结合超声心动评价冠心病

目的评价多排探测器螺旋CT(multi-detect row computed tomography MDCT)冠状动脉造影结合超声心动图对冠心病的诊断价值.材料方法收集2002年-2003年间54例临床怀疑或确诊的冠心病患者.所有患者进行MDCT(美国GE公司四排探测器螺旋CT)冠状动脉造影及超声心动(美国惠普公司5500超声)检查,其中30例患者在CT和超声检查前后2周范围内进行了选择冠状动脉血管造影检查(selected coronary angiography SCA),并以此为金标准,对MDCTA结合超声心动图诊断冠心病的价值进行评估.结果 54例患者中,42例患者显示56支冠状动脉病变;超声心动发现左室舒张功能障碍24例.30例SCA显示25例患者38支冠状动脉病变(左冠32支,右冠6支).MDCT冠脉造影结合超声心动法与SCA法检查结果一致性为:92.86%,X2检验:X2＝0.5;0.5>P>0.25,无显著性差异.结论 MDCT冠脉造影扫描可以作为冠状动脉狭窄性冠心病的初步筛选手段,结合超声心动图有利于提高诊断准确性.     

多层螺旋CT血管成像对急性肠系膜缺血性疾病的诊断价值

目的：评价多层螺旋cT血管成像（MSCTA）在急性肠系膜缺血性疾病中的诊断价值。方法：搜集2010年1月至2013年1月临床拟诊急性肠系膜缺血性疾病患者30例,分别行MSCTA检查,动脉期扫描延迟时间采用手动触发扫描,动脉期扫描结束后立即进行静脉期扫描,对比剂60mL,注射速率4mL／s。将获得薄层数据传至AW4．4工作站,分别进行多平面重组、曲面重组、容积再现和最大密度投影。结果：30例急性肠系膜缺血性疾病均经临床证实,其中20例行数字减影血管造影证实,MSCTA检查结果与DSA结果一致;肠系膜上动脉栓塞14例,合并动脉夹层2例,肠系膜下动脉栓塞6例,肠系膜上静脉8例,肠系膜下静脉2例。结论：MSCTA检查安全、快速、准确及价廉,是诊断急性肠系膜缺血性疾病的首选检查方法,对临床早诊断和早治疗有重要的指导价值。