一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

航空重力垂直梯度探测潜艇方法研究     

张志强  郑晗  崔银锋 《海洋测绘》2019,39(4):6-9

航空重力梯度测量属于被动探测,抗干扰能力强,如果能和其他探潜手段相配合将极大地提高航空搜潜的效率。针对航空重力梯度测量是否能够用于探测潜艇的问题,依据俄亥俄级弹道导弹核潜艇的结构特点,建立了适用于重力梯度计算的潜艇模型,分别给出了潜艇外壳、内部质量亏损产生的重力垂直梯度的计算方法,并对不同精度重力梯度仪可探测的潜艇重力垂直梯度值进行了计算,从航空反潜的角度对探测潜艇效果进行了分析。经计算表明,如航空重力梯度仪精度达到10~(-2) E,将具备一定的实际探潜效能;如精度达到10~(-4) E,反潜机搜索宽度可与现有航空磁性探潜相当。  相似文献   

井间电磁法综述     

刘四新  倪建福 《地球物理学进展》2020,(1):153-165

地面电磁法在实际工作中在很大程度上受到了探测深度和分辨率的限制.为了克服这些缺点,出现了井间电磁法.井间电磁法指的是在两个(或多个)钻孔中分别发射或接收电磁波信号,利用电磁波信号进行成像并探测井间物理性质的地球物理方法.由于发射机和接收机可以分别放置在很深钻孔中,其具有大透距、大探测深度的特点,因而广泛应用于工程环境物探、矿产勘查、石油勘探等中.针对不同的应用,产生了各具特点的一些特殊方法,包括井间无线电波成像、跨孔雷达、井间电磁成像.井间无线电波成像仪目前只测量电场强度数据,工作频率低,一般是单频的电磁波,频率范围通常在1 kHz至10 MHz.由于缺少走时数据修正射线路径,井间无线电波成像主要是进行基于直射线追踪的衰减层析成像.井间无线电波成像既可用于工程与环境地球物理也可用于找矿.跨孔雷达是钻孔雷达的一种探测方式,用高频电磁脉冲探测两个井孔间介电常数和电导率的变化.跨孔雷达层析成像也叫地质雷达CT,既可进行走时成像,还可进行衰减成像.一般来说,地质雷达CT的电磁波工作频率较高,中心频率通常在10 MHz和1 GHz之间,因此在分辨率指标上占有优势,跨孔雷达主要用于工程与环境地球物理.井间电磁成像采用更低的频率,测量复电磁信号,适合油气储集层监测,是一种地球物理前沿技术.经过在多个试验区初步试验表明,井间电磁成像是油藏研究的有效手段,可用于分析剩余油分布,寻找油气富集区,进而达到提高钻探高效井成功率和提高采收率的目的.本文详细介绍对比了这三种方法在理论和实践中的一些特点,并对未来的发展进行了展望.  相似文献   

非连续变形分析中块体大转动问题研究     

巩师林  凌道盛  胡成宝  钮家军 《岩土力学》2020,41(11):3810-3822

原有非连续变形分析（DDA）采用一阶近似后的位移增量表达式更新块体构形，推导相关子矩阵，且对不同时步计算出的应变增量直接叠加，当模拟的块体发生大转动时往往会产生较大误差。为考虑块体转动与变形的耦合作用，引入先变形、后转动的块体位移增量表达式。重新推导了惯性力子矩阵，将块体转动时的离心力与科氏力加到荷载矩阵中。计算时对应变分量及其相关变量进行坐标变换与修正，并采用新引入的位移增量表达式计算块体顶点位移，进行后接触修正与更新块体构形。数值算例表明，改进后的程序能够消除转动带来的误差，自动考虑了块体转动时离心力和科氏力引起的变形，应变计算精度更高。改进方法克服了块体体积自由膨胀、应变场畸变等问题，给出了合理的块体应变。  相似文献   

基于高压压汞技术和分形理论的致密砂岩储层分级评价标准     

程泽虎  李文浩  薛海涛  卢双舫  谭昭昭 《东北石油大学学报》2019,43(1)

利用铸体薄片与场发射扫描电镜观察,分析南襄盆地泌阳凹陷核桃园组三段致密砂岩储层储集空间;基于高压压汞技术和分形理论,对致密砂岩储层进行分类与分级。结果表明:研究区致密砂岩储层孔隙系统可划分为小孔(孔隙直径<0.1μm)、过渡孔(孔隙直径为0.1～1.0μm)、中孔(孔隙直径为1.0～3.0μm)和大孔(孔隙直径>3.0μm)。根据储层中不同类型的孔隙所占比例,可将致密砂岩储层分为4类,Ⅰ类(以大孔为主)、Ⅱ类(以中孔为主)、Ⅲ类(以过渡孔为主)和Ⅳ类(以小孔为主)。根据不同微观孔喉参数与储层物性的相关关系,可将致密砂岩储层分为4级,Ⅰ级储层孔隙度>10.0%,渗透率>1.000×10-3μm2;Ⅱ级储层孔隙度为5.0%～10.0%,渗透率为(0.200～1.000)×10-3μm2;Ⅲ级储层孔隙度为2.5%～5.0%,渗透率为(0.030～0.200)×10-3μm2;Ⅳ级储层的孔隙度<2.5%,渗透率<0.030×10-3μm2。利用储层分级评价标准选取最优质储层,为致密砂岩储层的油气勘探与开发提供指导。  相似文献   

不同非感应起电及感应起电参数化方案对青海东部一次雷暴云电荷结构影响的数值模拟研究          下载免费PDF全文

李江林  余晔  李万莉  李亚珺 《地球物理学报》2019,62(7):2366-2381

本研究利用加入起电、放电参数化方案的数值模式（Weather Research and Forecasting Model（Version 3.7.1），WRF3.7.1_ELEC），通过设计五组不同非感应起电及感应起电参数化方案敏感性试验，对发生在青藏高原东北部青海大通地区的一次雷暴过程进行模拟研究，对比分析了不同非感应起电机制及感应起电机制对雷暴云电荷结构的影响.结果表明：在雷暴云发展旺盛阶段，Saunders（S91）、Riming Rate（RR）、和Saunders和Peck（SP98）三种非感应起电方案模拟的雷暴云最低层均为负电荷区，而混合方案（Brooks and SP98，BSP）模拟的雷暴云最低层为正电荷区，主电荷区自下而上为"+-+-"排列的四层电荷结构.与甚高频辐射源定位法推算的结果对比，BSP方案模拟的本次高原雷暴云电荷结构更接近实际情况；几种不同非感应起电方案模拟的主电荷区外围与主电荷区电荷结构不同，说明在雷暴发展的不同阶段雷暴云的电荷结构是不同的；几种非感应起电方案模拟的电荷结构不尽相同，主要是由于霰、冰和雪粒子在不同高度所带电荷的极性及电量的大小不同，霰粒子的电荷密度对低层的影响较大，冰粒子和雪粒子的电荷密度对中上层的影响较大；加入感应起电机制后，雷暴云电荷结构分布几乎没有变化，但能使雷暴云发展旺盛阶段低层和中层的正负电荷区电荷密度有所加强.  相似文献   

地层格架正则化约束下的二维立体层析反演          下载免费PDF全文

张力起  杨锴  邢逢源  李振伟  倪瑶 《地球物理学报》2019,62(2):634-647

通过把地层格架信息作用于立体层析Fréchet导数矩阵,使得更新后的速度模型呈现出符合地质规律的块状特征.地层格架信息基于立体层析反演中得到的反射点位置进行非规则B样条插值拟合得到,因此在反演中它将会随着反射点位置的更新自然得到更新.与前人提出的保边缘层析算法或多层立体层析算法相比,本文提出的地层格架正则化无需引入混合正则化项或定义某种复杂的混合速度格式,更为直接也更容易实现.理论和实际数据算例证实了该正则化技巧的稳健性和可靠性,能够得到与实际地质构造特征更为一致的地质一致性反演结果.  相似文献   

非迁移潮作用下的电离层赤道异常区单峰现象及其经度分布          下载免费PDF全文

田耀宇  郝永强  张东和  肖佐 《地球物理学报》2019,62(11):4067-4081

当喷泉效应较弱而双峰结构发展不充分的时候，可能在赤道异常区仅能够观测到一个电子密度的峰值，称之为单峰现象.本文利用CHAllenging Minisatellite Payload卫星在2001-2010年的电子密度数据给出了单峰的发生规律，单峰在地方时早上08：00-10：00和下午16：00-19：00发生率高，发生位置在经度上呈现多波数分布，尤其在10：00-18：00明显：在分季时多呈现四波，而在冬至季时以三波为主.单峰发生多的经度，正好对应着双峰的结构特征较弱之处.究其原因，是非迁移潮的DE2和DE3分量调制了背景风场和大气发电机电场，在电场和喷泉效应减弱的经度，双峰结构难以形成时，就会表现为单峰结构.本文扩展了对单峰现象的地方时、季节和经度分布等规律的了解，明确了非迁移潮在其中施加的影响，由此，单峰同双峰现象一样可以用于研究非迁移潮对热层-电离层的作用.  相似文献   

2019年长宁M_S6.0地震邻区分钟采样地电场分析     

王玮铭  谭大诚 《地震》2020,40(3):41-51

2019年6月17日四川宜宾市长宁县发生M_S6.0地震, 分析震中周边600 km内7个地电场台站观测分钟值数据, 获知震前2个月左右红格等6个台站地电场均存在不同方位之间、 同方位长短极距之间的相关系数持续降低, 6月17日这些相关系数基本降至最低, 反映出这6个场地的地电场在此次地震前后发生了非均匀变化; 同时, 这6个台站的地电场优势方位角出现了约45°或90°跃变, 并伴有显著发散或收窄现象。分析表明长宁M_S6.0地震前后, 这6个场地的多种地电场异常现象具有局部场地效应、 时间准同步性。  相似文献   

曲靖非相干散射雷达电离层F区日间电子温度变化特征初步分析   总被引：1，自引：0，他引：1       下载免费PDF全文

丁宗华  代连东  杨嵩  唐志美  吴健  许正文 《地球物理学报》2020,63(4):1282-1293

利用曲靖非相干散射雷达2017-2018年春夏季观测数据首次分析了电离层日间150~450km电子温度的地方时与高度变化特征及其与电子密度的相关性.发现hmF2及以上的电子温度在日出日落时具有两个峰值,在11∶00-16∶00LT之间变化较小,高度越高午后上升的时间越早;从150km开始迅速增加,在约220km达到最大值,然后开始降低,在约300~350km达到最小值,最后单调上升;200km以下电子温度与电子密度成正相关(主要由热传导控制),200~450km之间存在明显的反相关(光电离过程占主导),电子-离子温度差与电子密度对数之间存在近似线性关系,电子温度逐日变化与光电离因子的变化趋势相似,这种相关性在中午与午后更明显;以上结果与其他非相干散射雷达观测和电离层模型计算结果基本一致,但也存在一些差别,需要结合更多数据深入分析.  相似文献