发明名称 趋于低频趋势的多井测井曲线标准化方法
摘要 本发明公开了一种趋于低频趋势的多井测井曲线标准化方法,该方法包括以下步骤:1)在已知油气田工区内,对要标准化的多口测井进行井震标定,获得需要测井的纵波曲线、密度曲线和地震叠加数据;对各口井要进行标准化的测井曲线采用FIR数字滤波器进行滤波处理,只保留测井曲线的低频信息,得到低频趋势线;对多条测井的低频趋势线进行平均,求取油田工区内所有井公共的低频趋势线;以公共的低频趋势线为标准,对各口井要标准化的测井曲线进行校正,得到标准化曲线。本发明提供了一种不需要利用测井曲线的标准井和标准层就能完成的趋于低频趋势的多井测井曲线标准化方法。
申请公布号 CN104849763A 申请公布日期 2015.08.19
申请号 CN201510157358.X 申请日期 2015.04.03
申请人 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 发明人 黄麟云;郁飞
分类号 G01V1/48(2006.01)I 主分类号 G01V1/48(2006.01)I
代理机构 湖北武汉永嘉专利代理有限公司 42102 代理人 胡建平
主权项 趋于低频趋势的多井测井曲线标准化方法,其特征在于,包括以下步骤:1)在已知油气田工区内,对要标准化的多口测井进行井震标定,获得需要测井的纵波曲线、密度曲线和地震叠加数据;其中,井震标定公式如下:X(t)=R(t)*W(t);式中,X(t)为合成记录,R(t)为反射系数、W(t)为子波;其中R(t)由测井的纵波曲线和密度得到,W(t)从井旁道的地震中提取;2)根据各口井要进行标准化的测井曲线确定滤波器技术指标和FIR滤波设计,采用FIR数字滤波器进行滤波处理,只保留测井曲线的低频信息,得到低频趋势线;所述低频信息的低频范围为0~10Hz;3)对多条测井的低频趋势线进行平均,求取油田工区内所有井公共的低频趋势线;计算公共低频趋势线公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub></mrow><msub><mi>n</mi><mn>2</mn></msub></munderover><msub><mi>x</mi><mi>ij</mi></msub></mrow>]]></math><img file="FDA0000694179580000011.GIF" wi="274" he="146" /></maths>式中,j为第j条测井曲线(j=1,2…,N),i为某一条测井曲线的第i个采样点,x<sub>ij</sub>为第j条测井曲线的第i个采样点,<img file="FDA0000694179580000012.GIF" wi="59" he="72" />为公共低频趋势线;4)以公共的低频趋势线为标准,对各口井要标准化的测井曲线进行校正;具体如下:测井曲线标准化矩阵如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>M</mi><mi>c</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>3</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>iN</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>iN</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>&Sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>iN</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>N</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>iN</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>N</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000694179580000021.GIF" wi="1488" he="283" /></maths>求出矩阵特征值λ,找出最大的λ值,再求出λ最大值对应的特征向量h,h的各个分量即为每条测井曲线的加权因子;根据加权因子对各条曲线加权滤波,得到标准化测井曲线。
地址 430223 湖北省武汉市东湖高新区大学园路18号(华师科技园)