专利名称：一种岩性数据采集处理方法
技术领域：
本发明涉及地球物理勘探开发方法，具体为一种岩性数据采集处理方法，属于测井海量信息综合解释评价范畴及涉及地层岩性参数的数据采集测试分析领域，尤其适用于基于岩性数据来识别地下水层或油层或含可动水油层。

背景技术：
快速准确地评价砂岩油藏油水层一直是油田勘探开发科技人员要解决的技术难题。目前油水层评价主要有五类方法分别为取心法、测试法、录井法、图版法和综合法。
众所周知，利用取心或测试技术是评价油水层最直接，也是最有效的评价技术，但不能口口井都进行取心或每个产层都进行测试，这既不经济也不现实。后三种方法特别是图版法评价油水层经济高效，应用更为普遍。需要指出的是测井资料是应用任何方法评价油水层时都不能抛开的最普遍、最有价值的地层信息。
图版法通常是采用测井技术结合取心和测试成果，依据Archie方程制定油水层识别标准，但由于钻井周期长，泥浆对地层(包括油层和水层)污染严重，破坏了地层原有的电性特征，同时由于井与井之间的泥浆性能不统一，从而无法制定统一的油水层识别标准，造成油水层识别率低，给油田开发带来很大的损失，为解决由于泥浆污染对油水层评价造成的不利影响，开展了深浅侧向电阻率测井受泥浆污染后的电阻率评价技术研究。
泥浆侵入原状地层后由于含油饱和度发生变化，使地层中的导电流体性质发生改变。如何客观准确地描述地层流体性质的变化是解决泥浆侵入原状地层后客观求取地层电阻率和识别流体特征的关键。
一大批专家学者在此方面做了大量的研究工作，但多少都存在一些不足，如专利号为CN1243958A的“储层流体性质测量方法及设备”，该专利是通过记录设备发出的两种以上频率的复合电流在储层中产生对应频率的可测量信号，通过对记录信号的对比，认为幅度差值比例小的井段储层流体是水，幅度差值比例大的井段储层流体是油或气。该专利对于受污染油层或水层判别会出现较大偏差，当砂岩地层中的原状流体被泥浆部分或全部驱替时，原状地层(油层或水层)的流体性质均发生变化，应用该专利技术存在油层或水层误判现象。
再如，专利号为CN 101363315A的“定量化薄层地层的电阻率和含烃饱和度的方法”，该专利是根据各层的束缚水饱和度和总孔隙度及估算的各层的束缚水总体积，计算地层水平电阻率和垂直电阻率，将估算值与测定的水平电阻率和垂直电阻率进行比较，调整所估算的各层的束缚水饱和度并反复估算该值，直到估算值与所测定的垂直电阻率之间的差值降到所选的阀门值，由经过调整的各层的束缚水饱和度估算烃体积。该专利采用反复调整后的束缚水饱和度实现估算电阻率与实测电阻率的一致性，并用调整后的束缚水饱和度估算含烃体积，该专利把由于泥浆侵入造成的地层电阻率升高或降低都归结到束缚水饱和度上，从而直接造成采用调整后的束缚水饱和度估算的烃体积可信度降低。


发明内容
为了解决上面背景技术中存在的问题，本发明提出一种岩性数据采集处理方法，属于测井海量信息综合解释评价范畴及涉及地层岩性参数的数据采集测试分析领域，其基于电阻率的计算来识别地层岩性或岩层中流体属性。
依据本发明所述的技术方案，提供一种岩性数据采集处理方法，其包括以下步骤 1)利用地质勘探设备，可以采集地层岩性数据； 2)基于采集到的地层岩性数据，通过岩性数据分析设备获得地层中的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度； 3)基于采集到的地层岩性数据，通过岩性数据分析设备获得包括可动地层水电阻率Rwf和地层束缚水电阻率Rwi的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率； 4)根据在步骤2)中获得的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度，以及根据在步骤3)中获得的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率，来获得地层流体替换率和泥浆滤液分配系数；其中地层流体替换率是指在井眼周围的指定范围内(指侧向电阻率测井的探测深度)，侵入渗透性砂岩地层中的侵入液体积与原状砂岩地层中可动流体体积之比，侵入液体积是指泥浆和泥浆滤液体积之和；泥浆滤液分配系数是指侵入渗透性砂岩地层中的泥浆滤液与侵入液体积之比。
地层流体替换率表达式 VF＝Vmud+Wmf(2) 式中VF-侵入指定地层(侧向电阻率测井探测深度)中的泥浆和泥浆滤液总体积，m3(立方米)；Vmud-侵入指定地层(侧向电阻率测井探测深度)中的泥浆体积，m3(立方米)；Vmf-侵入指定地层(侧向电阻率测井探测深度)中的泥浆滤液体积，m3(立方米)；V-指定原状地层(侧向电阻率测井探测深)中可动流体体积，m3(立方米)；uf-地层流体替换率，0≤uf≤1，无因次； 泥浆滤液分配系数表达式 式中泥浆滤液分配系数取值范围0≤vf＜1。
5)根据地层流体替换率和泥浆滤液分配系数，可以判定出砂岩地层的属性，即地层是否是纯水层或纯油层或含可动水油层。
其中，泥浆和泥浆滤液侵入纯水层后，将孔隙中的可动水全部或部分驱替，此时地层中的可动水饱和度为Swc＝1-Swi，深浅侧向电阻率测井在纯水层中的视电阻率解释式为 式中RLLso-100％含水砂岩浅侧向视电阻率，Ω·m；RLLdo-100％含水砂岩深侧向视电阻率，Ω·m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向视电阻率测井表现出水层为“增阻”(递增或倒增)或“半增阻”泥浆侵入特征，即Ro＜RLLdo＜RLLso(附

图1为淡水泥浆侵入水层梯形趋势的第一示意图A)，或Ro＜RLLso＜RLLdo(附图2为淡水泥浆侵入水层梯形趋势的第二示意图B)，或RLLso＜Ro＜RLLdo(附图3为淡水泥浆侵入水层梯形趋势的第三示意图C)；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出水层为“减阻”泥浆侵入特征，即Ro＞RLLdo＞RLLso(附图4为盐水泥浆侵入水层梯形趋势的第四示意图D)，在满足上述四种情形之一的同时还要满足

或

的数值关系，需判定的砂岩地层是纯水层。
其中Ro为砂岩水层真电阻率；Rts为浅侧向电阻率测井实测电阻率；Rtd为深侧向电阻率测井实测电阻率。
其中，当泥浆和泥浆滤液侵入纯油层后，将孔隙中的石油部分或全部驱替，由于砂岩孔隙中无原始可动水存在，即Swf＝0，得到深浅侧向视电阻率测井在砂岩地层中的视电阻率解释式为 式中RLLs-砂岩油层浅侧向视电阻率，Ω·m；RLLd-砂岩油层深侧向视电阻率，Ω·m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
根据RLLs、RLLd、RLLso、RLLdo、Rt、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是纯油层。
进一步地，采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向电阻率测井表现出砂岩油层为“减阻”(递减)泥浆侵入特征，即存在Rt＞RLLd＞RLLs的数值关系时，附图5为淡水泥浆侵入油层梯形趋势的第五示意图E，需判定的砂岩地层是纯油层；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出油层为“减阻”(速减)泥浆侵入特征，即存在Rt＞RLLd＞＞RLLs(附图6为盐水泥浆侵入油层梯形趋势的第六示意图F)，在满足上述两种情形之一时，还要满足

或

的数值关系，需判定的砂岩地层是纯油层。
其中Rt为砂岩油层真电阻率；Rts为浅侧向电阻率测井实测电阻率；Rtd为深侧向电阻率测井实测电阻率。
另外地，对于含可动水油层，深浅侧向电阻率测井在砂岩地层中的视电阻率解释式为 式中R′LLs-含可动水油层浅侧向视电阻率，Ω·m；R′LLd-含可动水油层深侧向视电阻率，Ω·m；Swf-砂岩地层可动水饱和度，0≤Swf＜1；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
根据RLLd、RLLs、RLLso、RLLdo、R′LLs、R′LLd、Rt、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是含可动水油层。
进一步地，当RLLd＞R′LLd＞RLLdo或RLLs＞R′LLs＞RLLso，同时存在

或

的数值关系时，要判定的砂岩地层是含可动水油层。
优选地，根据岩心分析孔隙度和实测声波时差得到砂岩总孔隙度。
优选地，采用不受井眼环境影响的压汞资料确定束缚水饱和度。
优选地，地层水密度是地层水矿化度、地层温度和地层压力的函数。
优选地，可动地层水电阻率Rwf可通过地层水矿化度(S·C)转换为等效NaCl矿化度后，通过总矿化度与对应离子含量的等效系数关系图版查取(附图7为总矿化度与对应离子含量的等效系数关系图、附图8为地层水等效NaCl溶液电阻率与温度交会关系的示意图)。
优选地，纯水层束缚水饱和度和可动水饱和度之间的关系 1-Swi＝Swc 式中Swc-100％含水砂岩中可动水饱和度，0≤Swc＜1；Swi-砂岩束缚水饱和度，0＜Swi≤1。
使用本发明所述的基于一种岩性数据采集处理方法来识别地下纯水层或纯油层或含可动水油层，并对射孔井段进行选择，剔除纯水层，选择纯油层对要生产的井段实施射孔，这样操作可有效地提高油藏的单井开采效率，降低开发成本，有利于油藏经济高效地开发。本发明评估了泥浆对井眼周围砂岩地层的污染程度，避免了油井在投入开发射孔时，无意中误射了纯水层或高含可动水的油层所造成的浪费，减少了对异常的误判。
附图简要说明 图1为淡水泥浆侵入水层梯形趋势的第一示意图A； 图2为淡水泥浆侵入水层梯形趋势的第二示意图B； 图3为淡水泥浆侵入水层梯形趋势的第三示意图C； 图4为盐水泥浆侵入水层梯形趋势的第四示意图D； 图5为淡水泥浆侵入油层梯形趋势的第五示意图E； 图6为盐水泥浆侵入油层梯形趋势的第六示意图F； 图7为总矿化度与对应离子含量的等效系数关系图； 图8为地层水等效NaCl溶液电阻率与温度交会关系的示意图； 图9为含油饱和度与油藏高度关系的示意图； 图10为原始溶解气油比与地层原油密度关系的示意图； 图11为应用本发明技术方案的葵东103井解释成果示意图； 图12为应用本发明技术方案的锦315井解释成果示意图。
图13为应用本发明技术方案的葵东102井解释成果示意图。

具体实施例方式 依据本发明所述的岩性数据采集处理方法来识别地下纯水层、纯油层或含可动水油层，具有识别准确、误判率低的优点。该岩性数据采集处理方法主要基于以下事实进行研究来获得的。
在侧向电阻率测井探测范围内，泥浆和泥浆滤液侵入原状砂岩地层，将孔隙中的地层水或石油部分或全部驱替，孔隙中地层水性质会因泥浆和泥浆滤液的侵入而发生变化。为了准确描述泥浆侵入砂岩地层后，砂岩地层中原状流体的变化，本发明提出了“地层流体替换率”和“泥浆滤液分配系数”的思想。“地层流体替换率”是指在一定范围内(指侧向电阻率测井的探测深度)，侵入渗透性砂岩地层中的侵入液体积与原状砂岩地层中可动流体体积之比，侵入液体积是指侵入渗透性砂岩地层中的泥浆和泥浆滤液体积之和；“泥浆滤液分配系数”是指侵入渗透性砂岩地层中的泥浆滤液与侵入液体积之比。
本发明采用如下技术方案来识别地下纯水层、纯油层或含可动水油层，即一种岩性数据采集处理方法，其包括以下步骤 1)利用地质勘探设备，可以采集地层岩性数据； 2)基于采集到的地层岩性数据，通过岩性数据分析设备获得地层中的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度； 3)基于采集到的地层岩性数据，通过岩性数据分析设备获得包括可动地层水电阻率Rwf和地层束缚水电阻率Rwi的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率； 4)根据在步骤1)中获得的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度，以及根据在步骤2)中获得的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率，来获得地层流体替换率和泥浆滤液分配系数；其中地层流体替换率是指在井眼周围的指定范围内(指侧向电阻率测井的探测深度)，侵入渗透性砂岩地层中的侵入液体积与原状砂岩地层中可动流体体积之比，侵入液体积是指泥浆和泥浆滤液体积之和；泥浆滤液分配系数是指侵入砂岩地层中的泥浆滤液与侵入液体积之比。
地层流体替换率表达式 VF＝Vmud+Fmf(2) 式中VF-侵入指定地层(侧向电阻率测井探测深度)中的泥浆和泥浆滤液总体积，m3(立方米)；Vmud-侵入指定地层(侧向电阻率测井探测深度)中的泥浆体积，m3(立方米)；Vmg-侵入指定地层(侧向电阻率测井探测深度)中的泥浆滤液体积，m3(立方米)；V-指定原状地层(侧向电阻率测井探测深)中可动流体体积，m3(立方米)；uf-地层流体替换率，0≤uf≤1，无因次。
泥浆滤液分配系数表达式 式中泥浆滤液分配系数的极限取值范围为0≤vf＜1。
5)根据地层流体替换率和泥浆滤液分配系数，可以判定出砂岩地层的属性，即地层是否是纯水层或纯油层或含可动水油层。
其中，泥浆和泥浆滤液侵入水层后，将孔隙中的可动水全部或部分驱替，此时地层中的可动水饱和度为Swf＝1-Swi，深浅侧向电阻率测井在水层中的视电阻率解释式为 式中RLLso-100％含水砂岩浅侧向视电阻率，Ω·m；RLLdo-100％含水砂岩深侧向视电阻率，Ω·m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
当采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向视电阻率测井表现出水层为“增阻”(递增或倒增)或“半增阻”泥浆侵入特征，即Ro＜RLLdo＜RLLso，或Ro＜RLLso＜RLLdo，或RLLso＜Ro＜RLLdo，同时存在

或

的数值关系；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出纯水层为“减阻”泥浆侵入特征，即Ro＞RLLdo＞RLLso，同时存在

或

的数值关系，需判定的砂岩地层是纯水层。
其中Ro为砂岩水层真电阻率；Rts为浅侧向电阻率测井实测电阻率；Rtd为深侧向电阻率测井实测电阻率。
其中，当泥浆和泥浆滤液侵入纯油层后，将孔隙中的石油全部或部分驱替，由于砂岩孔隙中无原始可动水存在，即Swf＝0，得到深浅侧向视电阻率测井在砂岩地层中的视电阻率解释式为 式中RLLs-砂岩油层浅侧向视电阻率，Ω·m；RLLd-砂岩油层深侧向视电阻率，Ω·m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
其中Rt为砂岩油层真电阻率，Rts为浅侧向电阻率测井的实测电阻率，Rtd为深侧向电阻率测井的实测电阻率。
根据RLLs、RLLd、RLLso、RLLdo、Rt、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是纯油层。
进一步地，采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向电阻率测井表现出砂岩油层为“减阻”(递减)泥浆侵入特征，即存在Rt＞RLLd＞RLLs的数值关系，同时存在

或

的数值关系，需判定的砂岩地层是纯油层。
另外地，对于含可动水油层，深浅侧向电阻率测井在砂岩地层中的视电阻率解释式为 式中R′LLs-含可动水油层浅侧向视电阻率，Ω·m；R′LLd-含可动水油层深侧向视电阻率，Ω·m；Swf-砂岩地层可动水饱和度，0≤Swf＜1；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
根据RLLs、RLLd、RLLso、RLLdo、R′LLs、R′LLd、Rt、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是含可动水油层。
进一步地，当深浅侧向电阻率存在RLLd＞R′LLd＞RLLdo或RLLs＞R′LLs＞RLLso的数值关系，同时存在

或

的数值关系，要判定的砂岩地层是含可动水油层。
优选地，根据岩心分析孔隙度和实测声波时差得到砂岩总孔隙度。
优选地，采用不受井眼环境影响的压汞资料确定束缚水饱和度。
优选地，地层水密度是地层水矿化度、地层温度和地层压力的函数。
优选地，可动地层水电阻率Rwf可通过地层水矿化度(S·C)转换为等效NaCl矿化度后，通过总矿化度与对应离子含量的等效系数关系图版查取(附图7为总矿化度与对应离子含量的等效系数关系图、附图8为地层水等效NaCl溶液电阻率与温度交会关系的示意图)。
优选地，纯水层束缚水饱和度和可动水饱和度之间的关系 1-Swi＝Swc 式中Swc-100％含水砂岩中可动水饱和度，0≤Swc＜1；Swi-砂岩束缚水饱和度，0＜Swi≤1。
为了更进一步详细说明本发明的技术方案，下面描述如何得到本发明中采用各种岩性参数。但是本领域技术人员应当清楚，下面对各种岩性参数的解释和说明仅仅是示例性的；本领域技术人员也可以通过其他途径来获得砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度，以及地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率等。
1岩性参数 1.1砂岩总孔隙度 1986年J.P.Martin等人在Raymer-Hunt工作的基础上，提出了声波地层因素公式 式中φt-砂岩总孔隙度，f；Δt-砂岩声波时差，μs/m；Δtma-砂岩骨架时差，μs/m；x-岩性系数。
由于(12)式与地层电阻率因素

很相似，故把

称为“声波地层因素”。
(12)式可转换为 根据岩心分析总孔隙度(φc)和实测声波时差(Δt)，由上式变形后的公式形式回归求得骨架时差(Δtma)和岩性系数(x)，最终得到孔隙度解释方程。
logΔt＝logΔtma-xlog(1-φc)(14) 式中φc-实验室分析的砂岩总孔隙度，f。
1.2束缚水饱和度 利用J(Sw)函数对毛管压力曲线进行分类和平均。
式中σ，θ和Pc-分别为实验室内的界面张力、接触角和毛管压力；K-砂岩渗透率，mD；φ-砂岩有效孔隙度(实验室分析)，f。
若毛管压力试验为水银-空气系统，公式可转化为 J(Sw)函数确定含油饱和度的具体步骤 第一步，计算每块样品的C值。
对某一样品的毛管压力曲线而言，σ，cosθ，K和φ均为常数。为计算方便，令(16)式中的常数项为C，即 第二步，计算岩样每个测压点的J(Sw)函数值。
J(Sw)＝CPc (18) 第三步，统计不同饱和度区间的平均汞饱和度、平均J(Sw)函数和平均压力数据。
第四步，以纵坐标表示J(Sw)函数，横坐标表示汞饱和度(SHg)，点上J(Sw)点子。
如果点子集中，说明这些样品同属于一种孔隙结构类型，可以综合为一条代表该类储层的J(Sw)函数曲线。
第五步，求取该类储层的平均毛管压力曲线。
已知该类储层的平均渗透率和平均孔隙度，利用(17)式算出平均C，而后将J(Sw)函数曲线上任一含汞饱和度对应的J(Sw)值乘以平均C值的倒数，即可得到该点对应的平均毛管压力值(Pc)，进而得出该类储层的平均毛管压力曲线。
第六步，将实验室条件下的平均毛管压力曲线换算为油藏条件下的毛管压力曲线。
在实验室内测量的毛管压力曲线，虽然使用的样品是油藏储层实际岩心，但非润湿相和润湿相流体不能直接采用油藏条件下的流体，而不同组合的流体界面张力和润湿接触角不同，测得的毛管压力也各不相同。因此无论采用油-水或水银-空气那个系统测得的毛管压力曲线，在计算原始含水饱和度时，必须校正为油藏条件下的平均毛管压力曲线。
实验室毛管压力表达式和油藏毛管压力表达式分别为 式中σ，θ和Pc-分别为实验室内的界面张力、接触角和毛管压力； σR，θR和PcR-分别为油藏条件下的界面张力、接触角和毛管压力。
由(20)、(21)式联立得 实验室流体的界面张力和润湿接触角是已知的，而由于油藏流体处在地层较高温度和地层压力下，且油和水中有溶解气，因此油藏流体的界面张力很难求取。测量油藏条件下润湿接触角十分困难，在实际使用中，将公式(22)简化为 将实验室分析的参数(表1为界面张力、接触角实验室分析成果)代入(22)式得 表1界面张力、接触角实验室分析成果
第七步，将油藏条件下的毛管压力换算为油柱高度。
油藏的毛管压力由油水的重力差来平衡，而毛管压力可表示为 PcR＝H(ρwf-ρof)g(25) 式中H-油藏的自由水面以上高度，m；g-重力加速度，9.80m/s2。
将(25)式转化为SI制实用单位，求解含油高度为 式中PcR-油藏毛管压力，MPa；ρwf、ρof-分别为油藏条件下地层水密度和原油密度，g/cm3。
第八步，求取含油饱和度。
将(24)式代入(26)式得 由(27)式和油藏平均毛管压力曲线即可得到含油饱和度与油藏高度关系，附图9为含油饱和度与油藏高度关系的示意图，进而求得不同油藏高度下的含油饱和度和束缚水饱和度。
1.3地层水密度 地层水密度是地层水矿化度、地层温度和地层压力的函数，确定地层水密度采用以下方法。
ρw2＝1-1.87×10-5t-1.063×10-6t2 (30) ρw3＝1-1.40×10-5t-2.40×10-6Pi (31) 式中的地层水矿化度(S·C)为实测值，地层温度(t)是根据研究地区的实测地层温度与垂直深度(D/m)建立的统计方程计算 式中地层压力(Pi/MPa)是根据实测地层压力与垂直深度(D/m)建立的统计方程计算 Pi＝0.10+0.01D (33) 1.4地层原油密度 在地层条件下，由于原油中不同程度地溶有天然气，使地层原油密度的变化与溶解气含量密切相关。根据实验室提供的实际样品分析数据，建立了原始溶解气油比与地层原油密度相关统计方程(299组分析数据)，附图10为原始溶解气油比与地层原油密度关系的示意图，相关系数为0.9090。
ρof＝1.30736-0.29186logRsi (34) Rsi＝321(Boi-1)1.02 (35) 式中ρo-地面原油密度(实际分析值)，g/cm3；Boi-原油体积系数，无因次；Rsi-原始溶解气油比(由PVT分析或由公式求取)，m3/m3；D-垂直深度，m。
1.5地层水电阻率 综合地层水电阻率(Rw*)和可动地层水电阻率(Rwf)可由实验室分析获取 式中Ro-100％含水砂岩真电阻率，Ω·m；φt-砂岩总孔隙度，f。
可动地层水电阻率(Rwf) 地层束缚水电阻率(Rwi)求取方法 另外，可动地层水电阻率(Rwf)还可通过地层水矿化度(S·C)转换为等效NaCl矿化度后，通过总矿化度与对应离子含量的等效系数关系图获得，附图7为总矿化度与对应离子含量的等效系数关系图、附图8为地层水等效NaCl溶液电阻率与温度交会关系的示意图。
1.6泥浆及泥浆滤液电阻率 泥浆和泥浆滤液电阻率是根据实测泥浆电阻率(R′m)、泥浆密度(ρm)、泥浆温度(t′)等参数，采用D.W.Hilchie在1984年提出的公式和Schlumberger公司提供的公式计算求取 式中Rm-地层温度下泥浆电阻率，Ω·m；t′-井口实测温度，℃；t-地层温度，℃；R′m-井口实测泥浆电阻率，Ω·m；ρm-井口泥浆密度，g/cm3；Rmf-地层温度下泥浆滤液电阻率，Ω·m；z-转换系数。
1.7砂岩地层真电阻率 (1)纯水层真电阻率 纯水层束缚水饱和度和可动水饱和度之间的关系 1-Swi＝Swc (43) 式中Swc-100％含水砂岩地层中可动水饱和度，f；Swi-砂岩束缚水饱和度，f。纯水层电阻率计算方程 式中N-单位立方体砂岩的切片数，个。
(44)式化简后得 式中Ro-100％含水砂岩地层真电阻率，Ω·m。
(2)含可动水油层真电阻率 含可动水油层的束缚水饱和度、可动水饱和度及含油饱和度之间的关系 1-Swi-Sof＝Swf (46) 根据(46)式，得含可动水油层真电阻率方程 式中R′t-含水油层真电阻率，Ω·m。
(3)纯油层真电阻率 1-Swi＝Soi (48) 根据(48)式，得纯油层真电阻率方程 式中Rt-砂岩油层真电阻率，Ω·m。
1.8地层流体替换率和泥浆滤液分配系数 (1)地层流体替换率 VF＝Vmud+Vmf 式中VF-侵入指定地层(指侧向电阻率测井探测深度)的泥浆和泥浆滤液总体积，m 3；Vmud-侵入指定地层(指侧向电阻率测井探测深度)中的泥浆体积，m3；Vmf-侵入指定地层(指侧向电阻率测井探测深度)中的泥浆滤液体积，m3；V-指定原状地层(指侧向电阻率测井探测深度)中可动流体体积，m3；uf-地层流体替换率，0≤uf≤1，无因次。
地层中泥浆和泥浆滤液的侵入量(VF)为钻井过程中滤失泥浆体积(VL)与钻井岩屑中可动流体体积(Vsc)之差 VF＝VL-Vsc (50) 钻井岩屑中可动流体体积(Vsc) 式中Vsc-钻井岩屑中可动流体体积，m3；d-钻头直径，一般采用直径9.5in的钻头，即9.5×0.0254＝0.2413m；h-渗透性地层厚度，m；φt-渗透性地层厚度权衡孔隙度，f；Swi-渗透性地层孔隙厚度权衡束缚水饱和度，f。
钻井过程中滤失泥浆体积(VL) 由于泥浆与岩屑、地层孔隙流体交换的体积量(VL)相等，则混合(composition)了地层水的泥浆密度可表示为 式中ρmc-混合了地层流体的泥浆密度(实际分析)，g/cm3；ρm-钻井泥浆密度(现场分析)，g/cm3；ρf-在地面条件下，原状地层中被驱替的孔隙流体密度(实际分析)，g/cm3；Vm-理论计算的泥浆体积，m3；VL-钻井过程中滤失的泥浆体积，m3；α-钻井过程中泥浆实际使用量与理论计算泥浆体积的比值(在实际钻井过程中所使用的泥浆体积量为理论计算值的2.0～3.0倍)，无因次。
由(52)式转换得 上式中的三种密度(钻井泥浆密度、原状地层中被驱替的孔隙流体密度、混合了地层流体的泥浆密度)均可由现场取样分析及时获取，而钻井过程所需要的泥浆体积量可通过理论计算。
式中π-圆周率，取值3.1415；L-钻井轨迹长度(对直井而言为钻井深度D)，m。
深浅侧向电阻率测井探测范围内的泥浆侵入量浅侧向电阻率测井探测深度为0.35m，深侧向电阻率测井探测深度为1.15m。
浅侧向电阻率测井探测范围内的可动流体体积 式中rs-浅侧向电阻率测井探测半径(等于探测深度加钻头半径)，m；Vs-浅侧向电阻率测井探测范围内的可动流体体积，m3。
深侧向电阻率测井探测范围内的可动流体体积(含浅侧向电阻率测井探测范围内的可动流体体积) 式中rd-深侧向电阻率测井探测半径(等于探测深度加钻头半径)，m；Vd-深侧向电阻率测井探测范围内的可动流体体积，m3。
在井筒中滤失泥浆总量的百分比ω分布在浅侧向电阻率探测范围内，余下的1-ω分布在深侧向电阻率探测范围内，即 式中ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；ω-侵入深浅侧向探测范围内的侵入液比例(分析获取)，0＜ω≤1。
(2)泥浆滤液分配系数表达式 式中vf-泥浆滤液分配系数，0≤vf＜1。
泥浆滤液分配系数是由实验室根据岩心分析获取。
1.9纯水层视电阻率方程 泥浆和泥浆滤液侵入砂岩地层后，将孔隙中的可动水部分或全部驱替，此时地层中的可动水饱和度为Swc＝1-Swi，则由(44)式可导出深浅侧向电阻率测井在砂岩地层中的视电阻率解释方程 式中侵入泥浆与泥浆滤液混合电阻率(Rmfs)是根据实验室分析的泥浆滤液分配系数计算求取。
1.10含可动水油层视电阻率方程 泥浆和泥浆滤液侵入含可动水油层后，将孔隙中的可动水和石油部分或全部驱替，由(44)式可导出深浅侧向电阻率测井在砂岩地层中的视电阻率解释方程 1.11纯油层视电阻率方程 泥浆和泥浆滤液侵入纯油层后，将孔隙中的石油部分或全部驱替，由于砂岩孔隙中无原始可动水存在，即Swf＝0，根据(44)式可导出深浅侧向电阻率测井在砂岩地层中的视电阻率解释方程 2油水层评价 通过对上述系列电阻率解释方程的分析，采用淡水泥浆钻井时(侵入液电阻率大于地层水电阻率)，泥浆和泥浆滤液侵入水层后，会使砂岩水层的电阻率增大；泥浆侵入油层会使砂岩油层的电阻率降低。而采用盐水泥浆钻井时(侵入液电阻率小于地层水电阻率)，泥浆侵入水层后，会使砂岩水层的电阻率减小；侵入油层会使砂岩油层的电阻率大幅度降低。
2.1水层评价方法 采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向电阻率测井表现出水层为“增阻”(递增或倒增)或“半增阻”泥浆侵入特征，即Ro＜RLLdo＜RLLso，附图1为淡水泥浆侵入水层梯形趋势的第一示意图A，或Ro＜RLLso＜RLLdo，附图2为淡水泥浆侵入水层梯形趋势的第二示意图B，或RLLso＜Ro＜RLLDo，附图3为淡水泥浆侵入水层梯形趋势的第三示意图C；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出水层为“减阻”泥浆侵入特征，即Ro＞RLLdo＞RLLso，附图4为盐水泥浆侵入水层梯形趋势的第四示意图D；在满足上述四种情形之一的同时，还要满足

或

的数值关系，这是识别纯水层的典型特征，附表2为流体性质识别模式。
2.2油层评价方法 采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向电阻率测井表现出砂岩油层为“减阻”(递减)泥浆侵入特征，即Rt＞RLLd＞RLLs，附图5为淡水泥浆侵入油层梯形趋势的第五示意图E；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出油层为“减阻”(速减)泥浆侵入特征，即Rt＞RLLd＞＞RLLs，附图6为盐水泥浆侵入油层梯形趋势的第六示意图F；在满足上述两种情形之一的同时，还要满足

或

的数值关系，这是识别纯油层的典型特征，附表2为流体性质识别模式。
2.3含可动水油层评价方法 采用淡水或或盐水泥浆钻井时，侵入液侵入含可动水油层后，出现RLLd＞R′LLd＞RLLdo或RLLs＞R′LLs＞RLLso的数值关系，同时存在

或

的数值关系，这是识别含可动水油层的典型特征，附表2为流体性质识别模式。
表2流体性质识别模式
本发明提出的砂岩地层电阻率计算及孔隙流体识别技术在辽河油田十余个区块应用后，效果显著，使油水层测井解释精度由原来的75％～80％提高到90％～95％，解释符合率提高了十五个百分点。下面列举实例加以说明。
实例1 葵东103井完钻井深1813.4m，井口温度26.9℃，井口泥浆密度1.160g/cm3，混合泥浆密度1.145g/cm3，属盐水泥浆。42号层埋深对应地层温度48.38℃，泥浆使用量207m3，地层温度下泥浆电阻率0.187Ω·m，地层温度下泥浆滤液电阻率0.410Ω·m，钻头直径9.5in(0.2413m)，总孔隙度31.1％，束缚水饱和度40.0％，综合地层水电阻率1.771Ω·m，束缚水电阻率2.833Ω·m，可动地层水电阻率1.417Ω·m，100％纯水层电阻率5.70Ω·m，泥浆滤液分配系数0.80，浅侧向地层流体替换率15.1％，深侧向地层流体替换率1.2％，侵入液65％分布在浅侧向探测范围，由此计算该层含油时(未受泥浆污染)真电阻率为22.80Ω·m，计算受泥浆污染后的水层深浅侧向视电阻率RLLso为4.14Ω·m，RLLdo为5.58Ω·m，此时存在4.14(RLLso)＜5.58(RLLdo)＜5.70(Ro)的数值关系，同时存在

或

的数值关系，即为盐水泥浆侵入水层的“减阻”泥浆侵入特征，42号层应用本发明解释为水层。经试油证实是水层日产水86.1方，附图11为葵东103井解释成果示意图。
实例2 锦315井完钻井深1250m，井口温度20℃，井口泥浆密度1.100g/cm3，混合泥浆密度1.095g/cm3，属偏盐水泥浆。17号层埋深对应地层温度27.46℃，泥浆使用量114m3，地层温度下泥浆电阻率0.827Ω·m，地层温度下泥浆滤液电阻率1.824Ω·m，钻头直径9.5in(0.2413m)，总孔隙度31.0％，束缚水饱和度40.0％，综合地层水电阻率2.002Ω·m，束缚水电阻率3.204Ω·m，可动地层水电阻率1.602Ω·m，100％纯水层电阻率6.45Ω·m，泥浆滤液分配系数0.75，浅侧向地层流体替换率5.0％，深侧向地层流体替换率1.7％，侵入液30％分布在浅侧向探测范围，由此计算该层含油时(未受泥浆污染)地层真电阻率为25.80Ω·m，计算污染后地层深浅侧向视电阻率RLLs为22.06Ω·m，RLLd为24.68Ω·m，此时存在22.06(RLLs)＜24.68(RLLd)＜25.80(Rt)的数值关系，同时存在

或

的数值关系，即污染后油层“减阻”泥浆侵入特征，17号层原解释为水层，采用本发明解释为油层，经试油证实锦315井17号层是油层累产油61.632方，附图12为锦315井解释成果示意图。
实例3 葵东102井完钻井深2306m，井口温度18℃，井口泥浆密度1.190g/cm3，混合泥浆密度1.181g/cm3，属盐水泥浆。50号层垂深对应地层温度60.93℃，泥浆使用量263m3，地层温度下泥浆电阻率0.087Ω·m，地层温度下泥浆滤液电阻率0.190Ω·m，钻头直径9.5in(0.2413m)，总孔隙度28.7％，束缚水饱和度45.0％，综合地层水电阻率2.239Ω·m，束缚水电阻率3.471Ω·m，可动地层水电阻率1.736Ω·m，100％纯水层电阻率7.80Ω·m，泥浆滤液分配系数0.90，浅侧向地层流体替换率6.9％，深侧向地层流体替换率1.5％，侵入液40％分布在浅侧向探测范围，由此计算该层含油时(未受泥浆污染)真电阻率为26.87Ω·m，受污染后计算地层深浅侧向视电阻率RLLso为5.38Ω·m，RLLdo为7.16Ω·m，此时存在5.38(RLLso)＜7.16(RLLdo)＜7.80(Ro)的数值关系，同时存在

或

的数值关系，属盐水泥浆侵入水层的“减阻”泥浆侵入特征，50号层原解释为油层，采用本发明解释是水层，经单层试油证实50号层是水层测液面731.3m，日产水97.80方，附图13为葵东102井解释成果示意图。
如上述，已经清楚详细地描述了本发明提出的岩性数据采集处理方法来识别地下水层、油层或含可动水油层。尽管本发明的优选实施例详细描述并解释了本发明，但是本领域普通的技术人员可以理解，在不背离所附权利要求定义的本发明的精神和范围的情况下，可以在形式和细节中做出多种修改。
权利要求
1.一种岩性数据采集处理方法，其包括以下步骤
1)利用地质勘探设备，可以采集地层岩性数据；
2)基于采集到的地层岩性数据，通过岩性数据分析设备获得地层中的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度；
3)基于采集到的地层岩性数据，通过岩性数据分析设备获得包括可动地层水电阻率Rwf和地层束缚水电阻率Rwi的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率；
4)根据在步骤2)中获得的砂岩总孔隙度、束缚水饱和度、地层水密度和地层原油密度，以及根据在步骤3)中获得的地层水电阻率、泥浆及泥浆滤液电阻率和砂岩地层真电阻率，来获得地层流体替换率和泥浆滤液分配系数；其中地层流体替换率是指在井眼周围的指定范围内(指侧向电阻率测井的探测深度)，侵入渗透性砂岩地层中的侵入液体积与原状砂岩地层中可动流体体积之比，侵入液体积是指侵入的泥浆和泥浆滤液体积之和；泥浆滤液分配系数是指侵入渗透性砂岩地层中的泥浆滤液与侵入液体积之比。
地层流体替换率表达式
VF＝VMud+Vmf(2)
式中VF-侵入指定地层(侧向电阻率测井探测深度)中的泥浆和泥浆滤液总体积，m3(立方米)；Vmud-侵入指定地层(侧向电阻率测井探测深度)中的泥浆体积，m3(立方米)；Vmf-侵入指定地层(侧向电阻率测井探测深度)中的泥浆滤液体积，m3(立方米)；V-指定原状地层(侧向电阻率测井探测深)中可动流体体积，m3(立方米)；uf-地层流体替换率，0≤uf≤1，无因次；
泥浆滤液分配系数表达式
5)根据地层流体替换率和泥浆滤液分配系数，可以判定出砂岩地层的属性，即地层是否是纯水层或纯油层或含可动水油层。
2.根据权利要求1中所述方法，其中泥浆和泥浆滤液侵入纯水层后，将孔隙中的可动水全部或部分驱替，此时地层中的可动水饱和度为Swc＝1-Swi，深浅侧向电阻率测井在纯水层中的视电阻率解释式为
式中RLLso-100％含水砂岩浅侧向视电阻率，Ω·m；RLLdo-100％含水砂岩深侧向视电阻率，Ω·m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi ≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
当采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向视电阻率测井表现出水层为“增阻”(递增或倒增)或“半增阻”泥浆侵入特征，即Ro＜RLLdo＜RLLso，或Ro＜RLLso＜RLLdo，或RLLso＜Ro＜RLLdo，同时存在
或
的数值关系；采用盐水泥浆钻井时，侵入液电阻率小于地层水电阻率，深浅侧向电阻率测井表现出水层为“减阻”泥浆侵入特征，即Ro＞RLLdo＞RLLso，同时存在
或
的数值关系，需判断的砂岩地层是纯水层。
其中Ro为砂岩水层真电阻率，Rts为浅侧向电阻率测井的实测电阻率，Rtd为深侧向电阻率测井的实测电阻率。
3.根据权利要求1中所述方法，其中当泥浆和泥浆滤液侵入纯油层后，将孔隙中的石油部分或全部驱替，由于砂岩孔隙中无原始可动水存在，即Swf＝0，得到深浅侧向电阻率测井在砂岩地层中的视电阻率解释式为
式中RLLs-砂岩油层浅侧向视电阻率，Ω·m；RLLd-砂岩油层深侧向视电阻率，Ω.m；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
根据RLLs、RLLd、RLLso、RLLdo、Rt、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是纯油层。
4.根据权利要求3中所述方法，其中采用淡水泥浆钻井时，侵入液电阻率大于地层水电阻率，深浅侧向电阻率测井表现出砂岩油层为“减阻”(递减或速减)泥浆侵入特征，即存在Rt＞RLLd＞RLLs或Rt＞RLLd＞＞RLLs的数值关系，同时存在
或
的数值关系，需判定的砂岩地层是纯油层。
5.根据权利要求1中所述方法，其中对于含可动水油层，深浅侧向电阻率测井在砂岩地层中的视电阻率解释式为
式中R′LLs-含可动水油层浅侧向视电阻率，Ω·m；R′LLd-含可动水油层深侧向视电阻率，Ω·m；Swf-砂岩地层可动水饱和度，0≤Swf＜1；Rwi-地层束缚水电阻率，Ω·m；Rwf-可动地层水电阻率，Ω·m；Rmf-地层温度下泥浆滤液电阻率，Ω·m；Rmfs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的侵入液综合电阻率，Ω·m；Swi-砂岩束缚水饱和度，0＜Swi≤1；ufs-在浅侧向电阻率测井探测范围内(探测深度0.35m)的地层流体替换率，0≤ufs≤1；ufd-在深侧向电阻率测井探测范围内(探测深度1.15m)的地层流体替换率，0≤ufd≤1；vf-泥浆滤液分配系数，0≤vf＜1；φt-砂岩总孔隙度，0＜φt≤47.64％。
根据RLLd、RLLs、RLLso、RLLdo、R′LLs、R′LLd、Rts和Rtd之间的数值关系，可以判定出砂岩地层是否是含可动水油层。
6.根据权利要求5中所述方法，当砂岩地层存在RLLd＞R′LLd＞RLLdo或RLLs＞R′LLs＞RLLso的数值关系，同时存在
或
的数值关系，要判定的砂岩地层是含可动水油层。
7.根据权利要求1-6中之任一所述方法，其中根据岩心分析孔隙度和实测声波时差得到砂岩地层总孔隙度。
8.根据权利要求1-6中之任一所述方法，其中采用不受井眼环境影响的压汞资料确定束缚水饱和度。
9.根据权利要求1-6中之任一所述方法，其中地层水密度是地层水矿化度、地层温度和地层压力的函数。
10.根据权利要求1-6中之任一所述方法，其中可动地层水电阻率Rwf可通过地层水矿化度(S·C)转换为等效NaCl矿化度后，通过总矿化度与对应离子含量的等效系数关系图版查取(附图7为总矿化度与对应离子含量的等效系数关系图、附图8为地层水等效NaCl溶液电阻率与温度交会关系的示意图)。
11.根据前述权利要求中之任一所述方法，其中纯水层束缚水饱和度和可动水饱和度之间的关系1-Swi＝Swc
式中Swc-100％含水砂岩中可动水饱和度，0≤Swc＜1；Swi-砂岩束缚水饱和度，0＜Swi≤1。
全文摘要
本发明涉及地球物理勘探开发的岩性数据采集处理方法，涉及地层岩性参数的数据采集测试分析领域，适用于基于岩性数据来识别地下水层或油层或含可动水油层。包括利用地质勘探开发设备采集地层岩性数据；通过分析设备对采集的地层岩性数据进行分析获得计算砂岩地层视电阻率所需的基础参数及地层流体替换率和泥浆滤液分配系数，通过比较受污染的砂岩地层深浅侧向视电阻率和真电阻率及实测电阻率之间的数值关系，判定砂岩地层是水层、是油层还是含可动水油层。本发明通过对井眼周围砂岩地层污染程度分析，避免了油井在投入开发射孔时误射水层所造成的浪费，减少了对异常层的误判，有效地提高油藏的单井开采效率，有利于油藏经济高效地开发。
文档编号G01V3/38GK101776771SQ20101010711
公开日2010年7月14日 申请日期2010年2月9日 优先权日2010年2月9日
发明者康志勇, 闫家宁, 李铁军, 孔令福, 郭晓军 申请人:康志勇