1.一种基于井地电位法定量描述储层中含水饱和度的油水分布识别方法，其特征在于，该方法主要包括以下步骤：
1)首先收集应用油田的地质参数和流体参数，对这些参数进行统计分析，方便后续模型的输入计算；
2)建立井地电位方程，通过将渗流方程与电位微分方程相耦合，得到油水两相井地电位数学模型一般式；
3)建立模型差分方程，通过隐式求压力、显示求饱和度、隐式求电位三个步骤，按顺序建立差分方程组并进行求解；
4)使用Visual Basic程序设计语言完成井地电位模型的设计与计算，即已知地层含水饱和度、地层水矿化度、孔隙度等地层电阻率相关参数，来求解地表电位数据；
5)应用井地电位法反演，即通过正演得到地表电位数据与地下目标层电阻率相关关系，在求得目标层电阻率后，利用阿尔奇公式反演储层处含水饱和度。
2.如权利要求1所述的基于井地电位法定量描述储层中含水饱和度的油水分布识别方法，其特征在于，所述步骤1)中收集的应用油田的地质参数和流体参数包含井地电位法模拟所需的基本油藏参数，具体为地层原始压力、生产井井地流压、流体渗流面积、水的粘度、生产井与注入井之间的距离、地层渗透率、流体流量、原油粘度、地层孔隙度、束缚水饱和度、模拟设定的网格长度、模拟设定的网格数目、模拟设定的时差、模拟时间、地层水电导率、通电电流、地层初始电位等。
3.如权利要求1所述的基于井地电位法定量描述储层中含水饱和度的油 水分布识别方法，其特征在于，所述步骤2)中建立井地电位方程的过程为：
首先讨论油水两相渗流模型的建立，考虑到源汇项则油层条件下质量守恒形式的数学模型的一般式为：
式中，l表示相，l＝o,w分别表示油相或水相；qlv＝ql/ρlsc为源汇项，代表地面标准状况下单位体积岩石中注入或采出l相流体的体积流量,m3/s；ρlsc为地面标准状况下l相流体密度，kg/m3；B1为l相流体体积系数；k为油藏多孔介质绝对渗透率,m2；krl为l相流体相对渗透率,m2；μl为l相流体粘度,Pa·s；P1为l相流体压力，MPa；ρl为l相流体密度，kg/m3；g为重力加速度，m/s2；D为标高，由某一基准面算起的垂直方向方向深度(海拔)，m；Φ为孔隙度；Sl为l相流体饱和度；t为时间，s。
根据阿尔奇公式，地层电阻率与储层孔隙度、含水饱和度以及地层水电阻率的关系为：
式中，R为地层电阻率，Ω·m；a为比例系数，在0.6-1.5之间变化；m为孔隙度指数，取值范围为1.2-2.0；n为饱和度指数，这是一个由实验得到的因子，一般取值为2；Rw为充填于孔隙中水的电阻率，Ω·m。
根据电场强度、电位、电流以及电导率的相关关系(式3至式6)推导出电位微分方程(式
7)：


J＝σ·E             (5)


式中，E为电场强度，v/m；U为电位，V；j为电流密度，A/m2；Ω是被包含于曲面τ中的区域；
L为一线电流源；I为电流,A；σ为电导率，1/(Ω·m)；Φt为总电位，v。
式(2)和式(7)联立，建立油藏含水饱和度与地层电位的相关联系，式(1)与式(7)耦合得到油水两相井地电位数学模型一般式。将三维模型简化至一维并做如下假设：不考虑岩石的压缩性；不考虑重力作用的影响；不考虑流体体积的变化；油水粘度为常数；不考虑毛管力；油水互不相溶均符合达西渗流；供电电极电流均匀，没有奇异性。
于是一维井地电位基本模型为：
式中，σ为地层电导率，s/m；I为供电电流，A；U为电位，V。
模型中存在4个未知量分别为压力、含水饱和度、含油饱和度以及地层电位，因此还需1个辅助方程：
sw+so＝1                 (11)
初始条件为：0≤x≤l时，
p(x,0)＝pi                (12)
sw(x,0)＝swc               (13)
u(x,0)＝0                 (14)
边界条件为：t≥0时，
qv|x＝0＝qwv＝qv             (15)
qv|x＝l＝qwv+qov＝qv          (16)
I|x＝0＝Ii                  (17)
以上式(8)至式(17)组成了一维油水两相井地电位问题的完整数学模型。通过有限差分法数值求解，数值模拟计算结果中包含每个位置的压力、饱和度、电位等值，可以清晰地看到它们随时间的变化和分布。
4.如权利要求1所述的基于井地电位法定量描述储层中含水饱和度的油水分布识别方法，其特征在于，步骤3)中建立模型差分方程的过程为：在油藏渗流领域，有限差分法优点最为突出，实用，方便，简单，准确。因此本文亦选择此法求解。通过隐式求压力、显示求饱和度、隐式求电位三个步骤，按顺序建立差分方程组并进行求解。首先隐式求压力，即将式(8)至式(10)进行隐式差分。假设一维网格采用网格间距相同，大小相同的块状网格，并且注水井定义在i＝1处的网格，生产井定义在i＝n处的网格。
隐式差分结果为：
式中，λ表示总流动系数。
经过非线性系数线性化显式处理，系数取上游权等中间参数的处理，最 终可得三对角矩阵。通过托马斯算法求解矩阵得到每个网格n+1时刻的压力值。然后显示求饱和度，将已n+1
求得的n+1时刻的压力值代入式8至式10中便得到只含饱和度Swi 一个未知量的方程，因此下面只需将式8至式10进行两次中心差分，即：
由于每个方程只有一个待求量，因此很容易解得每个网格n+1时刻饱和度的值。
在求出Swin+1分布的基础上，再隐式求解电位，通过式(2)易得每一个网格处的电导率σin+1，然后隐式差分电位方程：
可得矩阵方程：
。