一种用井轴电场分布求取地层真电导率的方法与流程

文档序号:11152384阅读:来源:国知局

技术特征:

1.一种用井轴电场分布求取地层真电导率的方法,其特征在于,包括如下步骤,

步骤1,在井眼中沿井轴测量电场强度的分布;测井仪器在井眼中每移动一个场点对应测量一个电场分布;

步骤2,将测量得到的电场强度分布转换为视电导率函数;

步骤3,对视电导率函数进行井眼校正;

步骤4,从经井眼校正后的视电导率函数值中选取正确的函数值,得到地层真电导率。

2.根据权利要求1所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,步骤1中,由如下公式得到场点对应电场分布中的电场;

式中,E为电场强度;j为复数中虚数单位;ω为圆频率;μ为介质的磁导率;π为圆周率;e为自然常数;x为场点位置矢量;x1为源点的位置矢量;J为电流密度;k0(x)为场点处的传播系数;是源点处的传播系数;▽为梯度算子。

3.根据权利要求2所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,当传播的介质为轴对称介质时,感应测井电场如下所示,

式中,I为发射电流,I=I0ejωt;nT为发射线圈匝数;S0为发射线圈面积;r为场点径向坐标。

4.根据权利要求2所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,步骤2中,将测量得到的电场强度分布转换为视电导率函数σA(x)时,设发射线圈位于点(0,0,h),则由如下公式得到视电导率函数;

式中,

<mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>re</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <msub> <mi>E</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>re</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mi>a</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

其中,σa(x)——井眼以外区域等效电导率;σA(x)——全空间等效视电导率,等于感应测井视电导率。

5.根据权利要求4所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,步骤3中,根据如下公式得到井眼以外所有地层的等效电导率σa1(x),作为在井眼以外空间等效的电导率,就是井眼电导率为σa1(x)时的视电导率;完成对视电导率函数进行井眼校正;

<mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>E</mi> <mn>01</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </msubsup> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <msub> <mi>r</mi> <mn>0</mn> </msub> </msubsup> <msub> <mi>F</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>&sigma;</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>dr</mi> <mn>1</mn> </msub> <msub> <mi>dz</mi> <mn>1</mn> </msub> </mrow>

式中,

<mrow> <msub> <mi>E</mi> <mn>01</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>ae</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <msub> <mi>F</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>j</mi> <mi>&omega;</mi> <mi>&mu;</mi> <mi>a</mi> </mrow> <mn>4</mn> </mfrac> <mfrac> <mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

其中,V0为井眼区域,r0为井眼半径;σ0为泥浆电导率;a为场点与井轴的距离;σa1(x)为井眼以外所有地层的等效电导率。

6.根据权利要求5所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,步骤4中,根据如下公式,从经井眼校正后的视电导率函数值中选取正确的函数值,得到地层真电导率;

<mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>E</mi> <mn>02</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mo>&Integral;</mo> <msub> <mi>h</mi> <mn>1</mn> </msub> <msub> <mi>h</mi> <mn>2</mn> </msub> </msubsup> <msubsup> <mo>&Integral;</mo> <msub> <mi>r</mi> <mn>0</mn> </msub> <mi>&infin;</mi> </msubsup> <msub> <mi>F</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>&sigma;</mi> <mo>(</mo> <msub> <mi>r</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>-</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mo>(</mo> <mi>x</mi> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>dr</mi> <mn>1</mn> </msub> <msub> <mi>dz</mi> <mn>1</mn> </msub> <mo>;</mo> </mrow>

式中,

<mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>re</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <msub> <mi>E</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <msub> <mi>E</mi> <mn>02</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>j&omega;&mu;In</mi> <mi>T</mi> </msub> <msub> <mi>S</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>re</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mi>x</mi> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mi>x</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

<mrow> <msub> <mi>F</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>j</mi> <mi>&omega;</mi> <mi>&mu;</mi> <mi>a</mi> </mrow> <mn>4</mn> </mfrac> <mfrac> <mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>j</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> </mrow> </msup> </mrow> <msup> <mrow> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> </mrow> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>j</mi> </mrow> <mo>)</mo> <msqrt> <mrow> <mn>0.5</mn> <msub> <mi>&omega;&mu;&sigma;</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <mo>|</mo> <mrow> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> </mrow> <mo>|</mo> <mo>)</mo> </mrow> <mo>;</mo> </mrow>

其中,h1为目的层的下界面纵坐标;h2为上界面纵坐标;H为地层厚度;σ(r1)为目的层电导率即真电导率,;σa2(x)为目的层以外所有地层的等效电导率。

7.根据权利要求6所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,令σa2(x)等于σa1(x),对无侵地层,有一个视电导率值等于真电导率,则得到相应的场点坐标。

8.根据权利要求6所述的一种用井轴电场分布求取地层真电导率的方法,其特征在于,令σa2(x)等于σa1(x),对有侵地层,有一个视电导率等于目的层侵入带电导率和真电导率的加权平均,则得到相应的场点坐标。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1