发明名称 一种基于电导探针阵列传感器的水平井多参数估计方法
摘要 本发明涉及一种基于电导探针阵列传感器的水平井多参数估计方法,可用于对低产液水平井产出剖面的动态监测。其特征是在使用电导探针阵列传感器对水平井内的多个参数进行测量之前,首先对电导探针阵列传感器进行标定,获得标定数据并建立标定数据库;在进行参数估计时,根据现场得到的电导探针阵列传感器的实际响应计算距离矩阵;最后根据距离矩阵中最小元素所处的位置找出标定数据库中的匹配向量,进而对电导探针阵列传感器方位角、水平井中的水相电导率和油水界面高度进行估计。本发明的有益效果是当采用电导探针阵列传感器进行水平井产出剖面测井时,可以对井内的油水界面高度、水相电导率以及电导探针阵列传感器的方位角同时在线测量,测量下限低、精度高、可靠性好。
申请公布号 CN103277084B 申请公布日期 2015.07.15
申请号 CN201310193498.3 申请日期 2013.05.23
申请人 北京航空航天大学 发明人 徐立军;徐文峰
分类号 E21B47/00(2012.01)I;E21B47/047(2012.01)I 主分类号 E21B47/00(2012.01)I
代理机构 代理人
主权项 一种基于电导探针阵列传感器的水平井多参数估计方法,其特征在于包括下列步骤:(1)在使用电导探针阵列传感器对水平井内的多个参数进行测量之前,首先对电导探针阵列传感器进行标定,具体步骤如下:A.进行标定之前,在电导探针阵列传感器方位角的取值区域上等间隔取m个点,分别记为θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>m</sub>,在油水界面高度的取值区域上等间隔取n个点,分别记为h<sub>1</sub>,h<sub>2</sub>,…,h<sub>n</sub>;B.在标定过程中,将电导探针阵列传感器的方位角调整为θ<sub>i</sub>,将油水界面高度调整为h<sub>j</sub>,其中i=1,2,…,m,j=1,2,…,n,测量此刻电导探针阵列传感器的响应,记为<img file="FDA00003234737500011.GIF" wi="122" he="90" />在标定过程中,保持水相电导率为常数,记为σ<sup>c</sup>,反复执行步骤B,测量每种条件下电导探针阵列传感器的响应,得到标定数据库,用矩阵表示为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>G</mi><mi>c</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>G</mi><mn>1,1</mn><mi>c</mi></msubsup></mtd><mtd><msubsup><mi>G</mi><mn>1,2</mn><mi>c</mi></msubsup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msubsup><mi>G</mi><mrow><mn>1</mn><mo>,</mo><mi>n</mi></mrow><mi>c</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>G</mi><mn>2,1</mn><mi>c</mi></msubsup></mtd><mtd><msubsup><mi>G</mi><mn>2,2</mn><mi>c</mi></msubsup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msubsup><mi>G</mi><mrow><mn>2</mn><mo>,</mo><mi>n</mi></mrow><mi>c</mi></msubsup></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msubsup><mi>G</mi><mrow><mi>m</mi><mo>,</mo><mn>1</mn></mrow><mi>c</mi></msubsup></mtd><mtd><msubsup><mi>G</mi><mrow><mi>m</mi><mo>,</mo><mn>2</mn></mrow><mi>c</mi></msubsup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msubsup><mi>G</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mi>c</mi></msubsup></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA00003234737500012.GIF" wi="551" he="349" /></maths>(2)对油水界面高度、电导探针阵列传感器方位角和水相电导率进行参数估计,具体步骤如下:A.在测量现场获取电导探针阵列传感器的实际响应,记为G<sup>m</sup>;B.计算距离矩阵D,表示为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>d</mi><mn>1,1</mn></msub></mtd><mtd><msub><mi>d</mi><mn>1,2</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>d</mi><mrow><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mn>2,1</mn></msub></mtd><mtd><msub><mi>d</mi><mn>2,2</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>d</mi><mrow><mn>2</mn><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>d</mi><mrow><mi>m</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>m</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>d</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00003234737500013.GIF" wi="514" he="335" /></maths>其中矩阵D中元素d<sub>i,j</sub>的计算公式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mrow><mo>|</mo><mo>|</mo><mfrac><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>c</mi></msubsup><msub><mrow><mo>|</mo><mo>|</mo><msubsup><mi>G</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>c</mi></msubsup><mo>|</mo><mo>|</mo></mrow><mo>&infin;</mo></msub></mfrac><mo>-</mo><mfrac><msup><mi>G</mi><mi>m</mi></msup><msub><mrow><mo>|</mo><mo>|</mo><msup><mi>G</mi><mi>m</mi></msup><mo>|</mo><mo>|</mo></mrow><mo>&infin;</mo></msub></mfrac><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA00003234737500014.GIF" wi="578" he="200" /></maths>C.确定矩阵D中最小元素的位置,记为(r,s),其中r表示矩阵D中最小元素所在的行号,s表示矩阵D中最小元素所在的列号,然后在标定数据库G<sup>c</sup>中找出(r,s)位置上的元素,定义为匹配向量,记为<img file="FDA00003234737500015.GIF" wi="123" he="85" />D.利用公式<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>&theta;</mi><mo>~</mo></mover><mo>=</mo><msub><mi>&theta;</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><mover><mi>&sigma;</mi><mo>~</mo></mover><mo>=</mo><mfrac><msub><mrow><mo>|</mo><mo>|</mo><msubsup><mi>G</mi><mrow><mi>r</mi><mo>,</mo><mi>s</mi></mrow><mi>c</mi></msubsup><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><msub><mrow><mo>|</mo><mo>|</mo><msup><mi>G</mi><mi>m</mi></msup><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mfrac><mo>&CenterDot;</mo><msup><mi>&sigma;</mi><mi>c</mi></msup></mtd></mtr><mtr><mtd><mover><mi>h</mi><mo>~</mo></mover><mo>=</mo><msub><mi>h</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00003234737500021.GIF" wi="331" he="376" /></maths>确定各参数估计值,其中<img file="FDA00003234737500022.GIF" wi="140" he="79" />和<img file="FDA00003234737500023.GIF" wi="46" he="78" />分别为电导探针阵列传感器方位角、水平井中的水相电导率和油水界面高度的估计值。
地址 100191 北京市海淀区学院路37号
您可能感兴趣的专利