基于多面元约束反演各向异性参数的方法和装置制造方法及图纸

本发明专利技术公开了基于多面元约束反演各向异性参数的方法和装置。该方法包括：获得地震叠前方位道集数据体，地震叠前方位道集数据体中的数据表示为dataiii,jjj,kkk，其中下标iii表示面元序号，jjj表示每个面元内的道序号，kkk表示每道内的采样点序号；从所获得的地震叠前方位道集数据体中抽取多面元数据；基于所抽取的多面元数据的约束，反演得到各向异性参数。应用本公开能够在不牺牲反演精度的情况下提高反演得到的各向异性参数的稳定性。

【技术实现步骤摘要】

本专利技术涉及油气储层地震勘探领域，更具体地，涉及一种基于多面元约束反演各向异性参数的方法和一种基于多面元约束反演各向异性参数的装置。
技术介绍
裂缝型储层作为一种重要的油藏资源，已经成为地震勘探研究的热门课题。在碳酸盐岩、页岩等储层中，裂缝的发育都对油藏勘探开发起着重要作用。裂缝的形成受多种因素控制，其物理属性复杂，横向、纵向变化大，表现出很强的各向异性特征。因此获得稳定的各向异性信息对于检测裂缝特征非常重要。现有技术中，提高基于叠前数据反演得到的各向异性参数的稳定性的方法主要有两种，一种是在处理阶段，通过优化叠前方位道集处理流程，进而获得更稳定的反演结果；另一种是通过优化反演策略，进而获得更稳定的反演结果。现有的反演策略优化大多是通过多道叠加、降低反演公式阶数等手段，以适当牺牲反演精度的方式来提高反演结果的稳定性。
技术实现思路
本专利技术的目的是提出了一种新的反演策略优化方法，其能够在保持反演精度的前提下提高反演得到的各向异性参数的稳定性。本专利技术还提供了相应的装置。根据本公开的一个方面，提供了一种基于多面元约束反演各向异性参数的方法，该方法包括：获得地震叠前方位道集数据体，地本文档来自技高网...

【技术保护点】
一种基于多面元约束反演各向异性参数的方法，该方法包括：获得地震叠前方位道集数据体，地震叠前方位道集数据体中的数据表示为dataiii,jjj,kkk，其中下标iii表示面元序号，jjj表示每个面元内的道序号，kkk表示每道内的采样点序号；从所获得的地震叠前方位道集数据体中抽取多面元数据；基于所抽取的多面元数据的约束反演得到各向异性参数。

【技术特征摘要】
1.一种基于多面元约束反演各向异性参数的方法，该方法包括：获得地震叠前方位道集数据体，地震叠前方位道集数据体中的数据表示为dataiii,jjj,kkk，其中下标iii表示面元序号，jjj表示每个面元内的道序号，kkk表示每道内的采样点序号；从所获得的地震叠前方位道集数据体中抽取多面元数据；基于所抽取的多面元数据的约束反演得到各向异性参数。2.根据权利要求1所述的方法，其中，从所获得的地震叠前方位道集数据体中抽取多面元数据包括：通过滚动式抽取、每次抽取R条地震测线的数据以形成反演数组，其中，除第一次抽取外，每次抽取都只抽取一条新的地震测线的数据来替换当前反演数组中地震测线序号最小的地震测线的数据，其中R表示宏面元半径。3.根据权利要求1所述的方法，其中，抽取多面元数据以及基于多面元数据的约束反演得到各向异性参数包括：步骤201，初始化令i＝0，readbool＝1，linecdp＝linemin，其中linecdp表示当前计算的地震测线在地震叠前方位道集数据体中的序号，linemin表示最小地震测线序号，i用于指示所计算的地震测线标号，分别标记为第0、1、2、……条地震测线，realbool表示读数标记，realbool取0或1；步骤202，如果i＝0，直接跳转到步骤204；步骤203，初始化令bool＝0，以及搜索是否存在t∈[ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise0-1]满足公式1：orgn_linet＝linecdp，(公式1)其中，bool表示读数标记，bool取0或1，ise0用于标记读数的开始，ntrcdp表示每
\t次读入的面元的数量，有ntrcdp＝ncdp*(R*2+1)，ncdp表示每条地震测线包含的面元数量，orgn_linet表示面元标号为t的面元所在的地震测线的序号；如果存在满足公式1的t，则设置bool＝1，readbool＝1；如果不存在满足公式1的t，则设置i＝i-1，readbool＝0；步骤204，如果i≤R，则设置readbool＝0；步骤205，如果readbool＝0，直接跳转至步骤208；步骤206，搜索t∈[ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise0-1]中满足公式2的值：orgn_linek≠orgn_linet，其中如果k≠ntrcdp则k＝t+1，如果k＝ntrcdp则k＝0，(公式2)然后记录满足公式2的k，记录为ise1＝k，ise1用于标记读数的结束；步骤207，如果ise0≥ise1并且k∈[ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise1-1,ise1]，或者如果ise0＜ise1并且k∈[ise0,ise0+1,ise0+2,...,ise1-1,ise1]，则设置orgn_linek＝-1，orgn_cdpk＝-1，其中，orgn_cdpk表示面元标号为k的面元在地震叠前方位道集数据体中的面元序号；步骤208，如果i≤R，则设置ise0＝0，ise1＝0；步骤209，如果i＞0且i≤R，并且readbool＝0，则直接跳转至步骤211；步骤210，从地震叠前方位道集数据体dataiii,jjj,kkk中读取数据，如果ise0≥ise1则读取ntrcdp-ise0+ise1+1个面元的数据，如果ise0＜ise1则读取ise1-ise0+1个面元的数据，写到数组data_inii,jjj,kkk中，ii表示抽取后的地震数据的面元标号，如果ise0≥ise1则ii＝ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise1-1,ise1，如果ise0＜ise1则ii＝ise0,ise0+1,ise0+2,...,ise1-1,ise1；记录所读取的面元标号为ii的面元所在的地震测线序号orgn_lineii、面元标号为ii的面元在地震叠前方位道集数据体中的面元序号orgn_cdpii；步骤211，如果i＝0，初始化令bool＝0，以及搜索是否存在t∈[ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise0-2,ise0-1]满足上述公式1；如果存在，则令bool＝1，readbool＝1；如果不存在，则令i＝i-1，readbool＝0；步骤212，初始化令j＝cdpmin，其中cdpmin表示最小面元序号；搜索是否存在t∈[ise0,ise0+1,ise0+2,...,ntrcdp-1,ntrcdp,0,1,ise0-2,ise0-1]满足公式3：j-R≤orgn_cdpt≤j+R(公式3)linecdp-R≤orgn_linet≤linecdp+R；记录满足公式3的t值：data_cdprjj＝t，其中jj＝0,1,2,...,temp_int-1，用于表示满足公式3的面元标号，temp_int表示满足公式3的面元的数量；把地震叠前方位道集数据体中与data_cdprjj相对应的全部数据记录到数组中；当orgn_cdpt＝j且orgn_linet＝linecdp时，记录该t值：cal_i＝t；步骤213，在多面元数据的约束下进行各向异性参数反演；步骤214，输出反演得到的各向异性参数data_outl,kkk，其中l＝0,1,...,xn，xn表示经反演输出的数据体个数；步骤215，设置ise0＝ise1；如果j<cdpmax，则令j＝j+1，并返回步骤212，其中cdpmax表示最大面元序号；步骤216，如果i<linemax，则令i＝i+1，返回步骤202，其中linemax表示最大地震测线序号。4.根据权利要求3所述的方法，其中，在步骤213中，在多面元数据的约束下进行各向异性参数反演包括：采用预条件共轭梯度反演，其中，在预条件共轭梯度反演中的每次外部迭代前，计算多面元数据中每个面元数据的多道约束值矩阵，以及在针对每个面元进行预条件共轭梯度反演时，将相应的多道约束值矩阵加入到其反演过程中。5.根据权利要求4所述的方法，其中，计算多面元数据中每个面元数据的多道约束值矩阵包括：A，初始化令deltamoyjj,l,kkk＝0,deltamozjj,l,kkk＝0，其中jj＝0,1,2,...,temp_int-1，l＝0,1,...,xn-1，kk＝0,1,2,...,ns-1，deltamoyjj,l,kkk表示水平方向约束项，deltamozjj,l,kkk表示垂直方向约束项；初始化令p＝0.0001，p表示无量纲最小量；B，计算deltamoyjj,l,kkk＝RIjj,l,kkk-RIcal_i,l,kkk，其中RIjj,l,kkk为预条件共轭梯度反演的中间结果，被初始化为输入的初始模型，通过预条件共轭梯度反演的每次外部迭代来更新其结果，RIcal_i,l,kkk表示反演结果；C，计算deltamoyjj,l,kkk=(deltamoyjj,l,kkk)2(deltamoyjj,l,kkk)2+p2;]]>D，计算deltamozjj,l,kkk＝RIjj,l,kkk-RIcal_i,l,kkk-1，其中kkk＝1,2,...,ns-1；E，计算deltamozjj,l,kkk=(deltamozjj,l,...

【专利技术属性】
技术研发人员：肖鹏飞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11