CN107605477A - 一种确定井下泵抽取样油气突破和含水率稳定时间的方法 - Google Patents
一种确定井下泵抽取样油气突破和含水率稳定时间的方法 Download PDFInfo
- Publication number
- CN107605477A CN107605477A CN201710704338.9A CN201710704338A CN107605477A CN 107605477 A CN107605477 A CN 107605477A CN 201710704338 A CN201710704338 A CN 201710704338A CN 107605477 A CN107605477 A CN 107605477A
- Authority
- CN
- China
- Prior art keywords
- pumping
- moisture content
- oil gas
- time
- cycle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Sampling And Sample Adjustment (AREA)
Abstract
本发明提供了一种确定井下泵抽取样油气突破和含水率稳定时间的方法,基本原理为在电缆地层测试器对目的层流体进行泵抽取样的过程中,利用特别设计的井下流体监测模块和泵抽模块对泵抽的流体性质进行分析,提取出瞬时含水率和泵向状态信息,然后在matlab中编写程序对原始瞬时含水率进行周期化处理获得周期平均含水率。最后,利用相关数学算法编程实现了泵抽过程中油气突破时间和含水率稳定时间的自动拾取。通过采用经特别设计的井下流体监测和泵抽仪器提供的数据编程实现了泵抽油气突破时间和含水率稳定时间的自动精确确定,弥补了常规方法只能通过人为拾取油气突破时间和含水率稳定时间的不足,提高了工作效率,同时大大方便了后续的分析研究工作。
Description
技术领域
本发明涉及油气勘探技术领域,尤其涉及一种确定井下泵抽取样油气突破和含水率稳定时间的方法。
背景技术
随着电缆地层测试技术的不断完善和发展,在对海上油田进行油气资源评价时,电缆地层测试能快速获得储层的物性、压力数据和流体样品,节约大量的测试成本,降低作业风险,在油田勘探开发过程中的作用越来越大。目前油田研究人员和相关服务公司应用较多的是测压压力、流度和取样样品化验分析结果等,缺少对泵抽过程中井下监测的流动压力、流量、泵抽体积、含水率等动态资料的分析和应用。尤其在低渗储层泵抽过程中泥浆滤液驱替效率低,泵抽作业需要的时间长,在海上油田意味着占用更多的平台时间,耗费更高的作业费用。
发明内容
本发明提供了一种确定井下泵抽取样油气突破和含水率稳定时间的方法,将泵抽过程中的动态信息加以分析利用,提高资料应用价值。该方法将杂乱的瞬时含水率转化为有明显趋势的周期平均含水率数据,精确确定泵抽过程中油气突破时间和含水率稳定时间,减少人为拾取误差,方便后续分析和经验总结,提高现场作业效率,降低作业成本,提高测井解释评价的精度。
为解决上述技术问题,本发明提供了一种确定井下泵抽取样油气突破和含水率稳定时间的方法,包括以下步骤:
S1:通过常规测井曲线确定电缆地层测试工具泵抽取样的作业点深度,根据该点的测井渗透率并信息优化测前设计,通过电缆绞车将电缆地层测试工具下至设计的泵抽取样作业深度,仪器坐封开展预测试作业分析地层渗透性,调整泵抽初始泵速,开展连续泵抽作业,利用井下流体监测设备分析泵抽流体性质,获取含水率原始数据;
S2:将电缆地层测试工具泵抽模块和流体监测模块的原始数据进行解编,挑选出泵抽时间、泵速、泵向状态参数和瞬时含水率的原始数据;
S3:对原始数据进行规范化处理:将仪器记录的瞬时含水率中的非值统一赋值为0.0001或者0.9999;
S4:利用泵向状态参数确定泵抽周期的开始和结束时间:若泵向状态参数为111000111000…,则将状态参数由0变为1的时刻作为一个周期的开始时刻,同时也是上一个周期的结束时刻;
S5:将瞬时含水率、泵抽时间等数据按泵向状态参数进行周期化平均处理:将相同周期内的瞬时含水率、泵抽时间进行累加求和,然后获取周期平均值;
S6:确定油气突破时间:将同一泵抽周期内的瞬时含水率最大值和最小值进行对比,若二者之差大于0.1且连续3个周期都出现这一情况,则确定相应3个周期中第一个周期的时间点为油气突破时间;
S7:将周期平均含水率数组的最后5个元素互相做差,若差值大于0.05的点小于3个,则初步判断泵抽结束时含水率近稳定,若出现3个或更多差值大于0.05,则泵抽结束时含水率还未稳定;
S8:求取近稳定泵抽点的周期平均含水率数组最后5个元素的平均值,将平均含水率数组的元素从前往后逐个与平均值对比,如果连续3个点满足二者之差小于0.02,则确定相应3个点中第一个点对应的时间为泵抽过程中含水率开始稳定时间;
S9:在安装内存为4G,处理器为酷睿I7,64位操作系统或以上配置的电脑上安装matlab7.0版本;
S10:利用matlab软件按上述方法原理编译算法实现油气突破时间和含水率稳定时间的自动拾取,同时实现多泵抽点数据的批量处理,批量输出各点拾取的结果。
本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
将泵抽过程中的动态信息加以分析利用,提高资料应用价值。该方法将杂乱的瞬时含水率转化为有明显趋势的周期平均含水率数据,精确确定泵抽过程中油气突破时间和含水率稳定时间,减少人为拾取误差,方便后续分析和经验总结,提高现场作业效率,降低作业成本,提高测井解释评价的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中井下流体监测含水率随泵抽作业时间变化的关系图形。
具体实施方式
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
如图1所示,本实施例所述的一种确定井下泵抽取样油气突破和含水率稳定时间的方法,包括以下步骤:
S1:通过常规测井曲线(伽玛、电阻率、中子、密度等)确定电缆地层测试工具泵抽取样的作业点深度。根据该点的测井渗透率等信息优化测前设计。通过电缆绞车将电缆地层测试工具下至设计的泵抽取样作业深度,仪器坐封开展预测试作业分析地层渗透性,调整泵抽初始泵速。开展连续泵抽作业,利用井下流体监测设备分析泵抽流体性质,获取含水率原始数据;
S2:将电缆地层测试工具泵抽模块和流体监测模块的原始数据进行解编,挑选出泵抽时间、泵速、泵向状态参数和瞬时含水率的原始数据;
S3:原始资料规范化处理。将仪器记录的瞬时含水率中的非值统一赋值为0.0001或者0.9999;
S4:利用泵向状态参数确定泵抽周期的开始和结束时间。如泵向状态参数为111000111000…,则将状态参数由0变为1的时刻作为一个周期的开始时刻,同时也是上一个周期的结束时刻;
S5:将瞬时含水率、泵抽时间等数据按泵向状态参数进行周期化平均处理。将相同周期内的瞬时含水率、泵抽时间进行累加求和,然后获取周期平均值;
S6:确定油气突破时间。如果同一泵抽周期内的瞬时含水率最大值和最小值之差大于0.1(瞬时含水率单位为小数)且3个周期连续(避免个别非值点引起的误判)出现,则确定相应3个周期中第一个周期的时间点为油气突破时间;
S7:将周期平均含水率数组的最后5个元素互相做差,如果差值大于0.05的点小于3个,则初步判断泵抽结束时含水率近稳定,如果出现3个或更多差值大于0.05,则泵抽结束时含水率还未稳定;
S8:求取近稳定泵抽点的周期平均含水率数组最后5个元素的平均值(剔除异常点后),将平均含水率数组的元素从前往后逐个与平均值对比。如果连续3个点满足二者之差小于0.02,则确定相应3个点中第一个点对应的时间为泵抽过程中含水率开始稳定时间。
S9:在安装内存为4G,处理器为酷睿I7,64位操作系统或以上配置的电脑上安装matlab7.0版本;
S10:利用matlab软件按上述方法原理编译算法实现油气突破时间和含水率稳定时间的自动拾取,同时实现多泵抽点数据的批量处理,批量输出各点拾取的结果。
本发明的基本原理为在电缆地层测试器对目的层地层流体进行泵抽取样的过程中,利用专门设计的井下流体监测模块分析泵抽的流体性质,同时对仪器泵抽模块的泵向状态参数进行分析,提取出周期信息,然后利用自编的程序对仪器监测的原始瞬时含水率数据进行周期化处理获得周期平均含水率。最后,对周期平均含水率数组的元素进行数学分析,利用相关数学算法编程实现了泵抽过程中油气突破时间和含水率稳定时间的自动拾取。
在WW14井2656.2米深度进行电缆地层测试泵抽取样作业,预测试流度为5.5mD/cP。将流体监测模块的瞬时含水率数据(蓝色实线)按照仪器泵抽模块的泵向状态参数进行自动识别计算出每个泵抽周期的平均含水率(绿色三角形),然后将泵抽时间和瞬时含水率、周期平均含水率进行绘图,得到随泵抽时间变化的两种含水率变化图。利用瞬时含水率的变化按上述原则确定油气突破时间,利用周期平均含水率的变化按上述原则确定含水率基本稳定的时间。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (1)
1.一种确定井下泵抽取样油气突破和含水率稳定时间的方法,其特征在于,包括以下步骤:
S1:通过常规测井曲线确定电缆地层测试工具泵抽取样的作业点深度,根据该点的测井渗透率并信息优化测前设计,通过电缆绞车将电缆地层测试工具下至设计的泵抽取样作业深度,仪器坐封开展预测试作业分析地层渗透性,调整泵抽初始泵速,开展连续泵抽作业,利用井下流体监测设备分析泵抽流体性质,获取含水率原始数据;
S2:将电缆地层测试工具泵抽模块和流体监测模块的原始数据进行解编,挑选出泵抽时间、泵速、泵向状态参数和瞬时含水率的原始数据;
S3:对原始数据进行规范化处理:将仪器记录的瞬时含水率中的非值统一赋值为0.0001或者0.9999;
S4:利用泵向状态参数确定泵抽周期的开始和结束时间:若泵向状态参数为111000111000…,则将状态参数由0变为1的时刻作为一个周期的开始时刻,同时也是上一个周期的结束时刻;
S5:将瞬时含水率、泵抽时间等数据按泵向状态参数进行周期化平均处理:将相同周期内的瞬时含水率、泵抽时间进行累加求和,然后获取周期平均值;
S6:确定油气突破时间:将同一泵抽周期内的瞬时含水率最大值和最小值进行对比,若二者之差大于0.1且连续3个周期都出现这一情况,则确定相应3个周期中第一个周期的时间点为油气突破时间;
S7:将周期平均含水率数组的最后5个元素互相做差,若差值大于0.05的点小于3个,则初步判断泵抽结束时含水率近稳定,若出现3个或更多差值大于0.05,则泵抽结束时含水率还未稳定;
S8:求取近稳定泵抽点的周期平均含水率数组最后5个元素的平均值,将平均含水率数组的元素从前往后逐个与平均值对比,如果连续3个点满足二者之差小于0.02,则确定相应3个点中第一个点对应的时间为泵抽过程中含水率开始稳定时间;
S9:在安装内存为4G,处理器为酷睿I7,64位操作系统或以上配置的电脑上安装matlab7.0版本;
S10:利用matlab软件按上述方法原理编译算法实现油气突破时间和含水率稳定时间的自动拾取,同时实现多泵抽点数据的批量处理,批量输出各点拾取的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710704338.9A CN107605477B (zh) | 2017-08-17 | 2017-08-17 | 一种确定井下泵抽取样油气突破和含水率稳定时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710704338.9A CN107605477B (zh) | 2017-08-17 | 2017-08-17 | 一种确定井下泵抽取样油气突破和含水率稳定时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107605477A true CN107605477A (zh) | 2018-01-19 |
CN107605477B CN107605477B (zh) | 2020-07-24 |
Family
ID=61065424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710704338.9A Active CN107605477B (zh) | 2017-08-17 | 2017-08-17 | 一种确定井下泵抽取样油气突破和含水率稳定时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107605477B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109989743A (zh) * | 2019-05-06 | 2019-07-09 | 中海油田服务股份有限公司 | 一种确定泥浆滤液侵入深度的方法及系统 |
CN110162851A (zh) * | 2019-05-09 | 2019-08-23 | 中国海洋石油集团有限公司 | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 |
CN112858397A (zh) * | 2021-01-11 | 2021-05-28 | 中国海洋石油集团有限公司 | 一种地层流体中原油饱和度的测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2252331Y (zh) * | 1995-04-28 | 1997-04-16 | 锦州电子应用技术研究所 | 单井原油自动控制监测仪 |
US20070203681A1 (en) * | 2006-02-24 | 2007-08-30 | Saudi Arabian Oil Company | Monte carlo simulation of well logging data |
CN105587298A (zh) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | 多流管模式的水驱油井含水率反演方法 |
US20160356128A1 (en) * | 2013-03-18 | 2016-12-08 | Yevgeny B. Levitov | Power wave optimization for oil and gas extracting processes |
-
2017
- 2017-08-17 CN CN201710704338.9A patent/CN107605477B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2252331Y (zh) * | 1995-04-28 | 1997-04-16 | 锦州电子应用技术研究所 | 单井原油自动控制监测仪 |
US20070203681A1 (en) * | 2006-02-24 | 2007-08-30 | Saudi Arabian Oil Company | Monte carlo simulation of well logging data |
US20160356128A1 (en) * | 2013-03-18 | 2016-12-08 | Yevgeny B. Levitov | Power wave optimization for oil and gas extracting processes |
CN105587298A (zh) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | 多流管模式的水驱油井含水率反演方法 |
Non-Patent Citations (2)
Title |
---|
E.B.DUSSAN V.: "A robust method for calculating formation mobility with a formation tester", 《SPE》 * |
张国栋等: "东海西湖凹陷低孔低渗储层电缆取样方法创新与应用", 《中国海上油气》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109989743A (zh) * | 2019-05-06 | 2019-07-09 | 中海油田服务股份有限公司 | 一种确定泥浆滤液侵入深度的方法及系统 |
CN110162851A (zh) * | 2019-05-09 | 2019-08-23 | 中国海洋石油集团有限公司 | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 |
CN110162851B (zh) * | 2019-05-09 | 2022-11-22 | 中国海洋石油集团有限公司 | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 |
CN112858397A (zh) * | 2021-01-11 | 2021-05-28 | 中国海洋石油集团有限公司 | 一种地层流体中原油饱和度的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107605477B (zh) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8141632B2 (en) | Method for hydraulic fracture dimensions determination | |
US9176255B2 (en) | Permeability prediction systems and methods using quadratic discriminant analysis | |
BR112015009197B1 (pt) | Método e sistema para realizar uma operação de perfuração | |
US10294785B2 (en) | Data extraction for OBM contamination monitoring | |
CN110397402B (zh) | 钻井方法及装置 | |
CN107605477A (zh) | 一种确定井下泵抽取样油气突破和含水率稳定时间的方法 | |
CN108982320A (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
CN106599442B (zh) | 基于综合录井参数的随钻储层物性识别评价方法和装置 | |
US11008835B2 (en) | Horizontal reservoir description systems | |
US3899926A (en) | Method and apparatus for continual compilation of a well data log | |
WO2006120366A1 (en) | Methods for analysis of pressure response in underground formations | |
US10386215B2 (en) | Method for monitoring a flow using distributed acoustic sensing | |
CN109538199A (zh) | 一种煤系地层含气量评价方法、装置及电子设备 | |
US11249218B2 (en) | Horizontal reservoir description systems | |
CN110863814B (zh) | 巨厚型生屑灰岩油藏单井分段比采液指数确定方法及装置 | |
CN114575833B (zh) | 一种水侵气藏动态储量预测方法、装置、设备及存储介质 | |
EP3012671A1 (en) | System and method for estimating properties of geological formations drilled using underreamer | |
WO2017048258A1 (en) | Real-time determination of formation water-filled porosity using dielectric measurement data | |
Jelmert | The effect of a distributed block length function on double porosity transitions during linear flow | |
CN118091786B (zh) | 基于井下测井资料的裂缝线密度校正方法和装置 | |
RU2731013C2 (ru) | Способ интерпретации краткосрочных гидродинамических исследований горизонтальных скважин и скважин с гидроразрывом пласта на неустановившемся режиме фильтрации | |
US11879834B2 (en) | Fluid identification using optical data measurements | |
US20170053389A1 (en) | Downhole imaging techniques using lamination analysis | |
CN113125412B (zh) | 基于激光元素信息的砂岩-泥岩识别图版识别岩性方法 | |
CN111221050B (zh) | 油水层识别方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant after: CHINA NATIONAL OFFSHORE OIL Corp. Applicant after: ZHANJIANG BRANCH OF CHINA NATIONAL OFFSHORE OIL Corp. Address before: Chinese CNOOC Building No. 25 Beijing 100010 Dongcheng District Chaoyangmen North Street Applicant before: CHINA NATIONAL OFFSHORE OIL Corp. Applicant before: ZHANJIANG BRANCH OF CHINA NATIONAL OFFSHORE OIL Corp. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |