简介:地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况,为偏移和成像提供有效的依据。在弹性波波场数值模拟中,旋转交错网格数值模拟(RSM)修改了标准交错网格数值模拟(SSM)方法,将同类的参数定义在同样的节点上,拓宽了稳定性条件的约束,但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础,结合实际波场传播的情况进行误差计算,对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟(VRSM),即是籽变阶数方法应用到RSM中,它可以很好地解决RSM在低速区域的数值频散问题,以及减少不必要的时间损耗;同时讨论了旋转交错网格的理论频散特性,并基于波场分离的方法分析了实际波场传播的频散误差,将原方法的应用范围由声波推广到剪切波,由理论值推广到时变值。在数值模拟试验中,VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析,验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数,能保证计算的精确性,并合理控制计算的时间。
简介:储层物性参数的改变伴随着弹性参数某种程度的改变,而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因,致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足,以利用弹性参数反演储层物性参数为目的,依据贝叶斯反演框架,我们建立了新的储层物性参数目标反演函数。首先,采用兼具确定性与随机性特点的统计岩石物理模型,考虑到不同弹性参数间的精度存在差异,引入权重系数,建立起储层物性参数与弹性参数间的加权统计关系。其次,基于这种加权统计关系,结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后,建立解的快速求解准则,求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明,该方法具有较高的反演效率,应用前景较好。
简介:摘要弹性支承块式无砟轨道是一种新型的、少维修甚至免维修的道床结构,是今后我国长大隧道中普遍采用的轨道结构。弹性支承块式无砟轨道施工精度、施工质量要求高,本文以某铁路长隧道为例具体介绍了其施工技术要点。
简介:摘要:为了对工业用地的经济价值和现实作用进行强化,相关部门应在弹性出让和租赁制度的基础上全面分析工业用地,在标准化规章条例的基础上解决工业用地弹性出让和租赁中存在的问题,以彰显工业用地价值。
简介:波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。
简介:地震反演是当今最广泛应用于含油气储层预测的技术之一,取得了很多很好的预测效果,但也有失败的例子,达不到区分岩性和识别流体的目的。而本文介绍的建立在岩石物理建模和分析基础上的地震弹性反演,可将含油气储层预测由定性向(半)定量推进一步。根据岩石物理建模和正演扰动分析,可深刻理解岩石物性参数与地震弹性参数之间的内在关系,进而寻找储层岩性、物性和含油气性的敏感地震弹性参数,建立起理论岩石物理解释图版。岩石物理分析结果和所建立的岩石物理解释图版分别用以指导地震反演和反演结果的解释,实现油气储层分布预测和流体检测的目的。文中的含气砂岩分布预测实例研究应用结果表明,这种方法较叠后地震反演储层预测技术具有无可比拟的优越性,效果更佳,效率更高。
简介:基于Cauchy先验分布的贝叶斯AVO反射率反演,可以产生类似于稀疏脉冲反演的结果,增强了对大反射系数的识别能力,反演过程中正则项和参数矩阵的计算都需要用到由上一步迭代所估计出的模型参数信息,因而反问题呈现非线性,此外该过程依赖于模型参数之间的线性统计关系,而且反演结果为反射率而非弹性参数,不利于储层预测和流体检测。贝叶斯AVO波形反演通过将反射率改写为弹性参数差分形式,直接从叠前地震数据中提取弹性参数,无需对弹性参数之间关系做线性假设。本文综合考虑上述两种方法的优点,在研究过程中仍然采用Cauchy先验分布,同时对贝叶斯AVO反射率的反演过程进行了修改,实现了基于Cauchy先验分布的AVO弹性参数弱非线性波形反演。本方法的提出有效避免了模型参数之间的线性假设,同时也能够从地震数据中直接反演得到纵、横波速度和密度。理论合成数据实验证明,此方法可以直接从叠前地震数据中提取弹性参数,而且在含有噪音的情况下也能得到比较准确的反演结果。
简介:采用高阶交错网格有限差分法算子构建了各向异性介质多分量弹性波动方程的高精度正演模拟和叠前逆时深度偏移的数值离散方程,在正演过程中采用最大绝对振幅能量法实现各向异性介质多分量初至旅行时的计算,并以此作为多分量双程弹性波叠前逆时成像条件,同时给出计算所需的稳定性条件、多分量叠前逆时成像的基本原理、吸收边界条件等,采用凹陷模型合成了共炮点道集,并进行多分量双程各向异性弹性波叠前逆时成像研究,并与对正演炮集进行多分量各向同性叠前逆时成像结果进行对比。计算结果表明,考虑了各向异性弹性波场的矢量特性,能够更为准确地实现叠前多分量弹性波场的成像问题,使地质层位中的断层、断点等复杂目标成像更加清晰准确,并且偏移成像精度较高,因此开展各向异性叠前逆时成像可为当前高精度岩性地震勘探提供方法指导。
简介:弹性参数在甜点区预测和页岩气的开发过程中扮演着重要的角色,因此研究等效弹性参数随页岩气储层属性的变化是一项很有意义的工作。研究中我们用x射线CT扫描技术获得了较为精确的页岩样品微观结构图像。从这些图像中,我们可以获得孔隙度和矿物的详细情况,据此,我们构建了三维数字岩心,并应用有限元法对弹性参数进行了数值模拟,其间深入考察了子样选取、网格划分、求解器类型以及边界条件等,该方法易于区别不同的矿物及其百分含量。本文重点研究孔隙度和干酪根含量对弹性参数的影响,计算结果表明,孔隙度和干酪根含量对弹性性质有较大的影响,当孔隙度和干酪根含量增加时,弹性模量降低,且当孔隙度小于0.75%左右、干酪根含量大于3%左右时弹性参数减小速率较缓。因为孔隙度仅仅为4.5%,孔隙中填充油或气对弹性参数的影响甚微。不同岩心样本具有不同的孔隙度和干酪根含量,传统岩石物理实验不仅昂贵而且费时,而数值模拟是基于数字岩心来计算弹性参数,更加经济、方便。本研究证实了将页岩样品的微观结构图像与弹性模量的计算相结合来预测页岩弹性参数的可行性。
简介:摘要:弹性层状体系理论是目前我国采用的柔性路面设计的理论模型。逐条分析讨论弹性层状体系理论基本假定的理论依据,探究弹性层状体系理论基本假定与实际道路的差异。
简介:有效应力〈σij〉的精确表达式,特别是导致有孔隙流体材料的弹性应变的应力〈P〉(有效应力)是基于假设之上的,仅对Hook定律有效,即〈σij〉=σij-αPδij和〈P〉=Pc-Pp,这里α=1-(K/Ks),Pc和Pp是围压和孔隙压力,K和Ks分别是干燥岩石(排水状态下)和岩石基质(岩石固体部分)的体积模量。Gccrtsma(1957年)和Skempton(1960年)在实验的基础上首次提出关于〈P〉的方程。该表达式虽然不直接依赖于孔隙度,但是当有效应力〈P〉等于围压时孔隙消失,因此K=Ks。如果使用〈σij〉方程中的有效应力定理,那么可根据没有孔隙压力的固体弹性模量确定一个具有孔隙压力的多孔固体的应变。有效应力表达式非常准确地描述了砂岩和花岗岩样品在围压和孔隙压力达到2.5kb(250MPa)时的应变。结果表明,该有效应力定理不适用于非弹性过程(如断裂)。
简介:地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层(PML)边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件(HABC)具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。
简介:线性弹性断裂力学(LEFM)被广泛应用于水力压裂设计和评价。有很多压裂作业的结果与基于LEFM的模型模拟结果并不一致,对这类水力压裂作业中的非线性和塑性效应已开展过多项研究。大部分针对LEFM在岩石和准脆性材料(如混凝土)中的应用研究仅限于LEFM在测量断裂韧性的室内试验中的应用。这些研究发现,LEFM不是描述准脆性材料断裂特性的有效工具。运用已发表文献中有关原地应力状态和巴奈特页岩力争性质的大量数据,本文重点研究了LEFM在水力压裂研究中的适用性。LEFM已成功应用于金属材料,有关其应用的局限性已开展过研究,而且这些研究结果已反映在ASTM标准中。在本文中,我们运用相同的方法研究了LEFM在巴奈特页岩水力压裂处理中的适用性。首先,为了有一个总体的概念,我们研究了不同的断裂方式,以找出静水孔隙压力和超高压储层条件下的边界应力。其次,运用近裂纹尖端解和全应力分布解定义和描述了过程区域(processzone)和奇异性主导区域(singularity—dominatedzone)。然后,根据全应力场分布,计算剪切过程区域和拉伸过程区域,并将其与奇异性主导区域的大小进行比较。最后,基于奇异性主导区域和过程区域的相对大小得出结论:LEFM不能精确描述巴奈特页岩的水力压裂特性。