Data extraction for OBM contamination monitoring

Patent metadata

Field	Value
Publication number	US-10294785-B2
Application number	US-201514975700-A
Country	US
Kind code	B2
Filing date	Dec 18, 2015
Priority date	Dec 30, 2014
Publication date	May 21, 2019
Grant date	May 21, 2019

Abstract

Disclosed are methods and apparatus obtaining in-situ, real-time data associated with a sample stream obtained by a downhole sampling apparatus disposed in a borehole that extends into a subterranean formation. The obtained data includes multiple fluid properties of the sample stream. The sample stream includes native formation fluid from the subterranean formation and filtrate contamination resulting from formation of the borehole in the subterranean formation. The obtained data is filtered to remove outliers. The filtered data is fit to each of a plurality of models each characterizing a corresponding one of the fluid properties as a function of a pumpout volume or time of the sample stream. Based on the fitted data, a start of a developed flow regime of the native formation fluid within the subterranean formation surrounding the borehole is identified.

First claim

What is claimed is:

1. A method comprising:
obtaining in-situ, real-time data associated with a sample stream obtained by a downhole sampling apparatus disposed in a borehole that extends into a subterranean formation, wherein the downhole sampling apparatus is in electrical communication with surface equipment disposed at a wellsite surface from which the borehole extends, wherein the obtained data includes a plurality of fluid properties of the sample stream, and wherein the sample stream comprises:
native formation fluid from the subterranean formation; and
filtrate contamination resulting from formation of the borehole in the subterranean formation;

via operation of a processor at at least one of the downhole sampling apparatus and the surface equipment:
using the processor to statistically filter the obtained data to remove outliers from the obtained data;
fitting the filtered data to each of a plurality of models each characterizing a corresponding one of the fluid properties as a function of a pumpout volume (V) or time (t) of the sample stream; and
identifying a start of a developed flow regime of the native formation fluid within the subterranean formation surrounding the borehole, wherein identifying the start of the developed flow regime is based on the fitted data; and

operating the downhole sampling apparatus based on the identified start of the developed flow regime of the native formation fluid.

2. The method of claim 1 wherein identifying the start of the developed flow regime based on the fitted data comprises:
generating a flow regime identification (FRID) plot comprising:
the fitted data corresponding to each of the fluid properties, relative to V or t; and
an exponential factor of V or t, relative to V or t; and

identifying from the FRID plot the minimum V or t at which the fitted data for each of the fluid properties substantially coincide with the exponential factor of V or t, wherein the start of the developed flow regime is the identified minimum V or t.

3. The method of claim 2 wherein the exponential factor of V or t is V −y or t −y , and wherein y is an adjustable parameter obtained based on the fitted data.

4. The method of claim 2 wherein:
the fluid properties comprise apparent optical density (OD) of the sample stream, apparent mass density (ρ) of the sample stream, and apparent gas-oil ratio (GOR) of the sample stream, each determined by the downhole sampling apparatus;
obtaining the in-situ, real-time data comprises obtaining OD data, ρ data, and GOR data;
filtering the obtained data comprises removing outliers from the OD data, the ρ data, and the GOR data;
fitting the filtered data comprises fitting the OD data, the ρ data, and the GOR data to corresponding models each characterizing a corresponding one of OD, ρ, and GOR as a function of V or t;
generating the FRID plot comprises plotting each of the fitted OD data, the fitted ρ data, and the fitted GOR data relative to V or t; and
identifying the start of the developed flow regime comprises identifying, from the FRID plot, the minimum V or t at which the plotted OD, ρ, and GOR data each substantially coincide with the exponential factor of V or t, wherein the start of the developed flow regime is the identified minimum V or t.

5. The method of claim 1 wherein filtering the obtained data to remove outliers comprises truncating the obtained data based on a range of one of the fluid properties.

6. The method of claim 5 wherein one of the fluid properties is apparent optical density (OD) of the sample stream, and wherein truncating the obtained data comprises truncating the obtained data to data points in which the OD is between about −1.0 and about 5.0.

7. The method of claim 5 wherein one of the fluid properties is apparent mass density (ρ) of the sample stream, and wherein truncating the obtained data comprises truncating the obtained data to data points in which the ρ is less than about 1.5 g/cm 3 .

8. The method of claim 5 wherein one of the fluid properties is apparent gas-oil-ratio (GOR) of the sample stream, and wherein truncating the obtained data comprises truncating the obtained data to data points in which the GOR is less than about 1,000,000 scf/bbl.

9. The method of claim 5 wherein filtering the obtained data further comprises truncating the obtained data based on a range of V or t.

10. The method of claim 9 wherein truncating the obtained data based on the range of V or t comprises truncating the obtained data to those data points in which the V or t is greater than the V or t at which native formation fluid is first detected in the sample stream.

11. The method of claim 5 wherein filtering the obtained data further comprises downsampling the obtained data.

12. The method of claim 5 wherein filtering the obtained data further comprises filtering the obtained data utilizing a median filter, a Winsorized mean filter, or a Hampel filter.

13. The method of claim 1 wherein the fluid properties comprise apparent optical density (OD) of the sample stream, apparent mass density (ρ) of the sample stream, and apparent gas-oil-ratio (GOR) of the sample stream, and wherein filtering the obtained data comprises:
truncating the obtained data to data points in which:
the OD is between zero and about 3.0;
the ρ is between about 0.1 g/cm 3 and about 1.2 g/cm 3 ;
the GOR is less than about 50,000 scf/bbl; and
the V or t is greater than the V or t at which native formation fluid is first detected in the sample stream;

downsampling the truncated data; and
filtering the downsampled data utilizing a Hampel filter.

14. The method of claim 1 wherein the fluid properties comprise apparent optical density (OD) of the sample stream, apparent mass density (ρ) of the sample stream, and apparent gas-oil-ratio (GOR) of the sample stream, and wherein filtering the obtained data comprises:
truncating the obtained data to data points in which:
the OD is between about zero and about 1.5;
the ρ is between about 0.6 g/cm 3 and about 0.9 g/cm 3 ;
the GOR is less than about 1,000 scf/bbl; and
the V or t is greater than the V or t at which native formation fluid is first detected in the sample stream; and

filtering the truncated data utilizing a median filter or a Winsorized mean filter.

15. The method of claim 1 wherein the filtered data includes a number n of data points corresponding to each model, and wherein fitting the filtered data to each model comprises:
with respect to each model, performing a plurality of iterations that each comprise:
adjusting a threshold fraction and/or a fit start/end of the model;
randomly selecting a sample of m data points from the n data points corresponding to the model;
fitting the model to the m data points utilizing the adjusted threshold fraction and/or fit start/end;
determining an error function for each of the m data points based on the fitting; and
selecting ones of the m data points that are inliers supporting the model for the current iteration based on the error function determined for each of the m data points;

determining an optimal threshold fraction and/or fit start/end based on which of the iterations has the highest percentage of inliers among the m data points of that iteration; and
linearly fitting the inliers selected during the iteration corresponding to the optimal threshold fraction and/or fit start/end.

16. A method comprising:
obtaining in-situ, real-time data associated with a sample stream obtained by a downhole sampling apparatus disposed in a borehole that extends into a subterranean formation, where

Assignees

• Schlumberger Technology Corp

Inventors

• Wang Kang
• Lee Ryan Sangjun
• Zuo Youxiang
• Venkatakrishnan Sai
• Sheng Xiaowei
• Chokthanyawat Suchart
• Jiang Lei
• Sun Lixiang

Classifications

• E21B49/10 Primary

  using side-wall fluid samplers or testers · CPC title

• E21B49/0875

  determining specific fluid parameters · CPC title

• E21B2049/085

  Fixed Constructions · mapped topic