WO2017053080A1 - Évaluation de volume en subsurface - Google Patents

Évaluation de volume en subsurface Download PDF

Info

Publication number
WO2017053080A1
WO2017053080A1 PCT/US2016/050842 US2016050842W WO2017053080A1 WO 2017053080 A1 WO2017053080 A1 WO 2017053080A1 US 2016050842 W US2016050842 W US 2016050842W WO 2017053080 A1 WO2017053080 A1 WO 2017053080A1
Authority
WO
WIPO (PCT)
Prior art keywords
cells
peak
cell
location
interpolated
Prior art date
Application number
PCT/US2016/050842
Other languages
English (en)
Inventor
Stephen Freeman
Original Assignee
Schlumberger Technology Corporation
Schlumberger Canada Limited
Services Petroliers Schlumberger
Geoquest Systems B.V.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Schlumberger Technology Corporation, Schlumberger Canada Limited, Services Petroliers Schlumberger, Geoquest Systems B.V. filed Critical Schlumberger Technology Corporation
Publication of WO2017053080A1 publication Critical patent/WO2017053080A1/fr

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/641Continuity of geobodies
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling
    • G01V2210/665Subsurface modeling using geostatistical modeling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/70Other details related to processing
    • G01V2210/74Visualisation of seismic data

Definitions

  • Patent Application Serial Number 62/222,820 filed on September 24, 2015, and entitled "RESEVOIR MODEL VERIFICATION SYSTEM,” which is incorporated herein by reference in its entirety.
  • embodiments are directed to evaluating the subsurface volumes generated by data interpolation procedures.
  • An interpolated subsurface volume is obtained using control cells and derived cells.
  • the derived cells include attribute values derived from the control cells.
  • a set of peak cells is selected, according to the attribute values.
  • a correlation coefficient is determined that characterizes a spatial relationship between the location of the peak cells and the location of the control cell. The correlation coefficient may be used to evaluate the interpolated subsurface volume.
  • FIG. 1 is a schematic view, partially in cross-section, of a field in which one or more embodiments of bulls eye identification analysis may be implemented.
  • FIG. 2 shows a schematic diagram of a system in accordance with one or more embodiments of the technology.
  • FIGs. 3, 4.1, 4.2, 4.3, and 4.4 show flowcharts in accordance with one or more embodiments of the technology.
  • FIGs. 5, 6.1, 6.2, 6.3, 6.4, 7.1, 7.2 and 8 show examples in accordance with one or more embodiments.
  • FIGs. 9.1 and 9.2 show E&P a computing system in accordance with one or more embodiments.
  • ordinal numbers e.g., first, second, third, etc.
  • an element i.e., any noun in the application.
  • the use of ordinal numbers is not to imply or create any particular ordering of the elements nor to limit any element to being a single element unless expressly disclosed, such as by the use of the terms "before”, “after”, “single”, and other such terminology. Rather, the use of ordinal numbers is to distinguish between the elements.
  • a first element is distinct from a second element, and the first element may encompass more than one element and succeed (or precede) the second element in an ordering of elements.
  • embodiments are directed to a method and system for analysis of the location and likelihood of extreme values clustering in interpolated data.
  • One or more embodiments include functionality to display a map of the interpolated data along with the data resulted from field measurements.
  • One or more embodiments include functionality to quantify the magnitude of the extreme values clustering effect on the interpolated data.
  • FIG. 1 depicts a schematic view, partially in cross section, of a field (100) in which one or more embodiments may be implemented.
  • one or more of the modules and elements shown in FIG. 1 may be omitted, repeated, and/or substituted. Accordingly, embodiments should not be considered limited to the specific arrangements of modules shown in FIG. 1.
  • a geologic sedimentary basin contains a subterranean formation (104).
  • the subterranean formation (104) may include several geological structures (106-1 through 106-4).
  • the subterranean formation (104) may include a shale layer (106-1), a limestone layer (106-2), a sandstone layer (106-3), and another shale layer (106-4).
  • a fault plane (107) may extend through the subterranean formation (104).
  • various survey tools and/or data acquisition tools are adapted to measure the subterranean formation (104) and detect the characteristics of the geological structures of the subterranean formation (104).
  • survey operations and wellbore operations are referred to as field operations of the field (100). These field operations may be performed as directed by the surface unit (112).
  • the surface unit (112) is communicatively coupled to the E&P computer system (118).
  • the data received by the surface unit (112) may be sent to the E&P computer system (118) for further analysis.
  • the E&P computer system (118) is configured to analyze, model, control, optimize, or perform management tasks of the aforementioned field operations based on the data provided from the surface unit (112).
  • the E&P computer system (118) is provided with functionality for manipulating and analyzing the data, such as performing simulation, planning, and optimization of production operations of the wellsite system A (114-1), wellsite system B (114-2), and/or wellsite system C (114-3).
  • the result generated by the E&P computer system (118) may be displayed for an analyst user to view the result in a two-dimensional (2D) display, three-dimensional (3D) display, or other suitable displays.
  • 2D two-dimensional
  • 3D three-dimensional
  • the surface unit (112) is shown as separate from the E&P computer system (118) in FIG. 1, in other examples, the surface unit (112) and the E&P computer system (118) may also be combined.
  • FIG. 1 shows a field (100) on the land
  • the field (100) may be an offshore field.
  • the subterranean formation (104) may be in the sea floor.
  • field data may be gathered from the field (100) that is an offshore field using a variety of offshore techniques for gathering field data.
  • the trajectory (i.e., path) of the wellbore is vertical.
  • FIG. 1 shows a wellbore as being vertical, the wellbore may be horizontal or may follow a winding path to the reservoir.
  • oilfield equipment may gather and store data in one or more well logs.
  • the data received by the surface unit (112) represents characteristics of the subterranean formation (104) and may include seismic data and/or information related to location of the horizon and fault surfaces or characteristics of the formation rocks like porosity, saturation, permeability, natural fractures, stress magnitude and orientations, elastic properties, etc., during a drilling, fracturing, logging, or production operation of the wellbore (103) at the wellsite system.
  • the data received by the surface unit (112) may be combined and included in a subsurface volume of data.
  • FIG. 2 shows a diagram of a system in which one or more embodiments of the technology may be implemented.
  • one or more of the modules and elements shown in FIG. 2 may be omitted, repeated, and/or substituted. Accordingly, embodiments should not be considered limited to the specific arrangements of modules shown in FIG. 2.
  • the E&P computer system (1 18) includes a data repository (210) for storing input data, intermediate data, and resultant outputs of the analysis data, an analysis tool (240), and a field equipment module (250) for performing various tasks of the field operation.
  • the data repository (210) may include one or more disk drive storage devices, one or more semiconductor storage devices, other suitable computer data storage devices, or combinations thereof.
  • content stored in the data repository (210) may be stored as a data file, a linked list, a data sequence, a database, a graphical representation, any other suitable data structure, or combinations thereof.
  • the data repository (210) includes functionality to store subsurface volume data (220) and one or more property maps (230).
  • the data repository (210) is any type of storage unit and/or device (e.g., a file system, database, collection of tables, or any other storage mechanism) for storing data.
  • the data repository (210) may include multiple different storage units and/or devices. The multiple different storage units and/or devices may or may not be of the same type or located at the same physical site.
  • the subsurface volume (220) contains data related to one or more attribute of the subterranean formation (104) presented in FIG. l.
  • the subsurface volume (220) includes volume cells, each volume cell being characterized by a single value of an attribute over the entire cell.
  • the volume cells may be created as a result of a gridding operation on the subsurface volume, each grid cell being a volume cell.
  • the volume cells may be control cells (221), interpolated cells (223), or peak cells (225).
  • a control cell (221) is a subsurface volume cell that has one or more attribute values extracted from a field measurement operation.
  • the field measurement operation directly obtains a measurement for an attribute of the location that is represented by the control cell.
  • the measurement for the attribute may be the attribute value or may be used to determine the attribute value.
  • the control cells are regarded as control cells in relation to the attribute.
  • a control cell may be a control cell for a first attribute and not a control cell for a second attribute.
  • a set of volume cells may be referred to as control cells regarding an attribute and not control cells regarding a different attribute.
  • an attribute for a volume cell may be an environment parameter
  • a pore fluid property e.g., fluid density, fluid viscosity, etc.
  • a surface or laver eeometrv Te e., surface depth, layer thickness, etc.
  • An interpolated cell (223) is a subsurface volume cell that has one or more attribute values derived from the attribute values of the control cell.
  • the attribute values of the interpolated cells are obtained from an interpolation operation of the control cells.
  • the attribute values of interpolated cells are derived from the control cells.
  • the measured data points are associated with a value of an attribute of the control cells (221).
  • the interpolated data points are associated with a value of an attribute of the interpolated cells (223).
  • a peak cell (225) is a subsurface volume cell that has the value of an attribute that represents a local extrema.
  • a local extrema is defined as a value of the attribute that is either the biggest or the smallest among the attribute values of volume cells situated within a threshold distance in each direction.
  • the threshold distance is defined in terms of number of volume cells in each direction of the peak cell.
  • the property maps (230) contain spatially defined attribute values that characterizes the subsurface volume.
  • the property map may be a density distribution map (231) or a correlation coefficient map (233).
  • the density distribution map (231) is a map that displays the level of clustering of the peak cells across the subsurface volume.
  • the density distribution map shows areas of convergence of the peak cells.
  • the density distribution map is obtained by projecting vertically the peak cells on a surface parallel to the horizon and interpolating the projected points to obtain a surface property map. The projected display representation of the peak cells is therefore characterizing the entire vertical space through the subsurface volume.
  • the correlation coefficient map (233) is a surface property map that displays the spatial distribution of the likelihood of extreme values clustering affected interpolation results.
  • the correlation coefficient of a point within the subsurface volume is defined as the attribute that quantify the degree of spatial separation of the peak cells from the control cells at the location of the point.
  • the correlation coefficient map is obtained by projecting vertically the peak cells on a surface parallel to the horizon and interpolating the correlation coefficient values of the projected points to obtain a surface property map.
  • the correlation coefficient volume (235) is a subsurface property volume that displays the spatial distribution of the likelihood of extreme values clustering affected interpolation results.
  • the correlation coefficient volume is obtained by interpolating the correlation coefficient values of the peak cells within the subsurface volume to obtain a subsurface property volume.
  • the E&P computer system (118) additionally includes an analysis tool (240) in accordance with one or more embodiments.
  • the analysis tool (240) includes a user interface (241), a cell interpolator (243), a volume calculator (245), a surface calculator (247), and a visualization creator (249). Each of these components is described below.
  • the user interface (241) corresponds to a graphical user interface (GUI) that includes functionality to receive input from a user and present or display graphical data to the user.
  • GUI graphical user interface
  • the user interface (241) may include a 3D volume data viewer, a 2D section of subsurface volume viewer, and parameters value input fields in accordance with one or more embodiments.
  • the cell interpolator (243) is software component that is configured to interpolate the attribute to obtain attribute values for interpolated cells from the attribute values of the control cells.
  • the cell interpolator (243) includes functionality to generate a set of calculated data points at locations within the subsurface volume where no data measurements exist.
  • the result of the interpolation process is a continuous property map or a continuous property volume.
  • the volume calculator (245) is a software component that is configured to execute calculations on volume cell attributes and execute spatial analysis within the subsurface volume. In one or more embodiments, the volume calculator (245) may perform calculations on multiple attribute data that cover the same area to study the relationships between different attributes.
  • the surface calculator (247) is a software component that is configured to execute calculations on surface properties maps and execute spatial analysis on the surfaces associated with the subsurface volume.
  • the visualization creator (249) is a software component that is configured to render 3D perspective visualizations of the subsurface volume components.
  • the visualization creator (249) is operatively connected to the user interface (241) and the data repository (210).
  • the visualization creator (249) includes functionality to generate display representations of subsurface geologic layers for presentation using coded maps (e.g., greyscale coded, color coded, pattern coded, etc.) linked to volume cells attributes.
  • the visualization creator (249) includes functionality to generate 3D surfaces for display using a color map linked to a surface attribute.
  • the visualization creator (249) may include functionality to perform final rendering and show the components of the subsurface volume.
  • the subsurface geometry may include, for example, natural faults, subsurface layers of different density, density distribution maps, and other such information.
  • the visualization creator (249) includes functionality to display subsurface field measurements or any other data points on top of the layers.
  • the visualization creator (249) includes functionality to highlight the control cells at the control cell's location within subsurface volume.
  • the visualization creator (249) displays a rendered image of the subsurface volume in the user interface (241). In such embodiments, a portion of the volume cells are shown on the display. For example, the portion may include a left portion, a rear portion, and a slice through subsurface volume with the volume cells of the front right section removed from display.
  • the rendering is defined by an observation point angle.
  • the observation point defines the perspective from which the volume cells are viewed.
  • the observation point is a position of a virtual camera or theoretical user from which the volume cells are viewed.
  • a top portion may be shown where the view is from above the surface of the earth looking downward at the subsurface volume.
  • the E&P computer system (118) includes the field equipment module (250) that is configured to generate a field operation control signal based at least on a result generated by the E&P computer system (118), such as based on the likelihood of unre liable subsurface information due to a bull-eyes effect in some regions of the field (100) depicted in FIG. 1 above.
  • the field operation equipment depicted in FIG. 1 above may be controlled by the field operation control signal.
  • the field operation control signal may be used to control drilling equipment, an actuator, a fluid valve, or other electrical and/or mechanical devices disposed about the field (100).
  • FIGs. 1 and 2 show a configuration of components, other configurations may be used without departing from the scope. For example, various components may be combined to create a single component. As another example, the functionality performed by a single component may be performed by two or more components.
  • FIGs. 3, 4.1, 4.2, 4.3, and 4.4 show flowcharts in accordance with one or more embodiments. While the various blocks in these flowcharts are presented and described sequentially, one of ordinary skill will appreciate that at least some of the blocks may be executed in different orders, or may be combined or omitted, and at least some of the blocks may be executed in parallel. Furthermore, the actions in the blocks may be performed actively or passively. For example, some actions may be performed using polling or be interrupt driven in accordance with one or more embodiments.
  • determination blocks may not use a processor to process an instruction unless an interrupt is received to signify that condition exists in accordance with one or more embodiments. As another example, determination blocks may be performed by performing a test, such as checking a data value to test whether the value is consistent with the tested condition in accordance with one or more embodiments.
  • FIG. 3 shows a flowchart for performing analysis of the location and likelihood of extreme values clustering in interpolated data in accordance with one or more embodiments.
  • an interpolated subsurface volume including volume cells is obtained.
  • the interpolated subsurface volume may be obtained from the data repository in accordance with one or more embodiments.
  • the volume cells include interpolated cells or control cells.
  • the obtained volume is spatially defined based on a set of location boundaries specified by the user or in a file. For example, the user may draw the boundaries in the user interface. As another example, the boundaries may be selected to match the boundaries of a seismic data cube.
  • a set of peak cells is selected from the volume cells in accordance with one or more embodiments.
  • the volume cells with an attribute value equal to a local extrema are selected as peak cells.
  • the selection is performed by applying a local moving data volume through the subsurface volume, wherein the moving volume is selected as to accommodate one or more volume cells at any time.
  • the moving data volume returns a flag when the centroid of the volume is either the maxima or the minima within the local volume.
  • the moving volume can be vertically thin (e.g. 1 volume cell thick within the layer of interest).
  • the lateral size of the volume is dependent on the scale of maxima or minima chosen but may be as small as a 3*3 volume cell window.
  • the identification of maxima and minima is performed through the volume of interest and results in a property that is the maxima and minima flags that identify the local highs and lows in the data. Further, once a local maxima or a local minima is detected, the volume cell that contains the centroid of the flagged moving volume is selected as a peak cell. In one or more embodiments, the presence and location of the peak cells may indicate where an extreme values clustering condition could be occurring.
  • a correlation coefficient that characterizes a spatial relationship between the location of the peak cells and the location of the control cells is determined in accordance with one or more embodiments.
  • the correlation coefficient is a measurement of the degree of clustering of the peak cells around the location of the control cells.
  • the degree of clustering is calculated based on the distance in space between the r>eak cell locations and the nearest control cell location.
  • an extreme values clustering condition present as a result of interpolation is characterized by a very short distance (i.e., less than a threshold distance) between peak cells and the control cells in comparison to the average distance of entire volume cells to the nearest control cell.
  • a condition where no or minimum extreme values clustering is present is characterized by an even spread in distances between peak cells and the control cells in comparison to the average distance of entire volume cells to the nearest control cell.
  • the correlation coefficient at each peak cell location can be interpolated to generate a subsurface property volume or a surface property map for further spatial analysis.
  • the correlation coefficient is presented for visualization.
  • the correlation coefficient subsurface property volume or the correlation coefficient surface property map is sent to the visualization creator to generate 3D visualizations of correlation coefficient property for display using a color map linked to the correlation coefficient values.
  • the correlation coefficient may be scaled to a range of extreme values clustering occurrence probability risk intervals and presented as a risk intervals volume or risk intervals map. For example, a user can analyze the risk interval volume and identify zones with unreliable data.
  • the displayed correlation coefficient may be used to perform oilfield operations. For example, based on the position of zones with umeliable data along with other field simulations and modeling, oilfield equipment may be physically adjusted to change the monitoring of the wellbore, inject fluids, and/or extract hydrocarbons from the wellbore.
  • FIG. 4.1 shows a flowchart for interpolating a set of control cells in accordance with one or more embodiments.
  • the control cells information is read by the cell interpolator.
  • the information about control cells includes location information and the value of the attribute being interpolated.
  • the interpolation function is selected in one or more embodiments.
  • the interpolation function is defined as the algorithm that predicts values for cells in a raster from a limited number of sample data points.
  • the interpolation function is used to predict values for any volume cell attribute at any location between the control cells.
  • the interpolation function may be used to predict values for subsurface depth, geothermal gradient, porosity, or seismic velocity.
  • values for an attribute at interpolated cells location are obtained based on the interpolation of the control cells by the interpolation function.
  • the interpolation operation results in a spatial correlation of the interpolated cells. In other words, cells that are in close proximity tend to have similar attribute value using the interpolation operation.
  • control cells and interpolated cells attribute values are combined to create an interpolated subsurface volume.
  • the subsurface volume contains entire control cells and entire interpolated cells.
  • the subsurface volume contains volume cells, wherein a volume cell is either a control cell or an interpolated cell. The process then proceeds to end.
  • FIG. 4.2 shows a flowchart for selecting the peak cells in accordance with one or more embodiments.
  • a volume cell is selected from the interpolated subsurface volume, and the volume's cell attribute value is read.
  • the predefined range may be defined in terms of the number of volume cells in the proximity of the selected volume cell.
  • the predefined range may be extended by different lengths in different directions.
  • the predefined range can be defined like 5 cells in lateral direction and 3 cells in vertical direction.
  • the size of the predefined range is a constant through the selection process and is considered therefore a moving volume as discussed in FIG 3.
  • the predefined range is centered on the selected cell location. In other words, the centroid of the moving volume falls within the selected volume cell.
  • the set of extrema values is selected from each of the minimum value of entire attribute values of cells within the selected range, and the maximum value of entire attribute values of cells within the selected range.
  • the selected volume cell attribute value is tested against each of extrema values of the set of extrema values within the predefined range. In one or more embodiments, a determination is made as to whether the selected volume cell attribute value is lower than or equal to the minimum value or the selected volume cell attribute value is higher than or equal to the maximum value. If the selected volume cell attribute value is either less than or equal to the minimum value, or greater than or equal to the maximum value, then the selected volume cell attribute value is defined as a local extrema.
  • Block 424 if the selected volume cell attribute value is identified as a local extrema, then the selected volume cell is assigned to the peak cells set. If the selected volume cell attribute value is not identified as a local extrema, then the process proceeds to Block 425.
  • FIG. 4.3 shows a flowchart for determining a correlation coefficient in accordance with one or more embodiments.
  • an interpolated cell location is read from the subsurface volume.
  • the cell location is defined by a set of X, Y, Z coordinates.
  • Block 432 the distance between the location of the interpolated cell and the location of the nearest control cell is calculated.
  • the interpolated cell is compared to the list of cells in the set of peak cells, and a determination is made as to whether the interpolated cell status is identified as being part of the peak cell set.
  • Block 434 if the interpolated cell is identified as a peak cell, the calculated distance to the peak cell is added to a peak cell total distance variable.
  • the peak cell total distance variable cumulates the distance of each peak cells to the peak cells nearest control cell.
  • the interpolated cell total distance variable cumulates the distance of each interpolated cells to the interpolated cells nearest control cell.
  • the interpolated cell total distance summarize entire interpolated cells distances to the interpolated cell's nearest control cell including the peak cells' distances.
  • Block 436 the subsurface volume is scanned, and a determination is made as to whether there is any cell whose location coordinate was not read. If an interpolated cell is found that was not previously selected, the process proceeds to Block 431.
  • an average peak cell distance is calculated based on peak cell total distance.
  • the average peak cell distance is obtained by dividing the peak cell total distance y the total number of peak cells.
  • an average interpolated cell distance is calculated based on interpolated cell total distance.
  • the average interpolated cell distance is obtained by dividing the interpolated cell total distance by the total number of interpolated cells.
  • a correlation coefficient is determined based on the ratio between the average peak cell distance and the average interpolated cell distance.
  • the correlation coefficient may be expressed as a percentage ratio.
  • the correlation coefficient is characterizing the degree of interpolated value clustering around the location of cells with attribute values used to generate the interpolated values.
  • a low value for a correlation coefficient is characterized by a high degree of interpolated value clustering.
  • a correlation coefficient at each interpolated cell location can be calculated and a subsurface property volume derived, wherein the property is the correlation coefficient. Further, the volume can be rendered by the visualization creator and displayed in the user interface for further spatial analysis.
  • results of the interpolation operation can be assessed in terms of the extreme values clustering probability of occurrence based on an established threshold for the correlation coefficient. For example, a correlation coefficient value of 5% or lower may suggest an unacceptable level of extreme values clustering probability.
  • a field operation may be initiated. The process then proceeds to end.
  • FIG. 4.4 shows a flowchart for determining a surface correlation coefficient in accordance with one or more embodiments.
  • Block 441 the location coordinates of the peak cells within the subsurface volume are read by the volume calculator.
  • the subsurface volume is divided into a set of vertical columns based on a predefined reference surface lattice.
  • the lattice is defined as a repetitive arrangement of polygons that cover a reference surface.
  • the polygons may be either regular or irregular.
  • the reference surface may be represented by the top of the subsurface volume, top of a geologic layer, or may be a horizontal plane overlapping the subsurface volume.
  • Block 445 the number of peak cells within each column's volume is counted for each column.
  • each column has a corresponding value specifying the number of peak cells in the column.
  • a density distribution map is computed for the entire subsurface volume, based on the number of peak cells within each column's volume.
  • the density map is derived by plotting the results of the peak cell count for each volume in the center of the corresponding draped polygon. Further, by interpolating the center points across the reference surface, a property raster map is created, wherein the property is the density of the peak cells.
  • the entire draped polygon may be assigned the value of the peak count, producing a polygon vector map.
  • a correlation coefficient that characterizes a spatial relationship between the control cells and the density distribution map is determined in the surface calculator.
  • the correlation coefficient may be defined as the percentage of the number of the peak cells to the total number of cells within the subsurface volume.
  • a correlation coefficient map may be calculated based on the density distribution map.
  • the correlation coefficient map is derived by dividing the density distribution map by the horizontal distance from the center point to the nearest control cell, for each draped polygon. Further, the surface property map can be rendered by the visualization creator and displayed in the user interface for further spatial analysis.
  • the blocks may be performed actively or passively.
  • some blocks may be performed using polling or be interrupt driven in accordance with one or more embodiments.
  • determination blocks may not use a processor to process an instruction unless an interrupt is received to signify that condition exists in accordance with one or more embodiments.
  • determination blocks may be performed by performing a test, such as checking a data value to test whether the value is consistent with the tested condition in accordance with one or more embodiments.
  • FIG. 5 shows a visualization of the spatial distribution of an interpolated property derived from the input well control data.
  • the interpolation resulted in a low extreme values clustering density.
  • the labeled pipes (501) depicted in FIG. 5 represents the locations of the input data associated with the control cells.
  • the input data is obtained from the wells (101) of an oilfield (100).
  • FIGs. 6.1 - 6.4 show examples of interpolation results.
  • FIG. 6.1 shows a sample map of an attribute values distribution. Each square is a representation of a volume cell and the grey shades represent the value of the attribute. In one or more embodiments, the attribute value is mapped to a grey scale.
  • the interpolated cell values are based on the control cells values.
  • the labeled pipes (610) depicted in FIG. 6.1, FIG. 6.2, FIG. 6.3, FIG. 6.4 represents the locations of the control cells. In one or more embodiments, the control cell values represent measurements from wells (101) of an oilfield (100).
  • FIG. 6.2 shows a sample map over the same area with the result of local minima detection.
  • FIG. 7.1 and 7.2 show examples of density distribution maps.
  • FIG. 7.1 shows an example of a non-clustered extreme values data set resulted from the interpolation of the control cells (710) shown as labeled pipes in FIG. 7.1.
  • the control cells have attribute values derived well measurements. For example, based on well logging data, a set of porosity measurements are derived along a set of vertical wellbore trajectories.
  • the control cells (710) are defined as the volume cells crossed by the wellbore trajectories. Based on the porosity values at the control cells location a porosity volume is interpolated for a reservoir layer similar with the sandstone layer (106- 3) of the oilfield (100).
  • the volume porosity is used to determine the volume of hydrocarbons that may be hosted by the sandstone layer in an exploration project and based on the degree of clustering a risk factor is assumed for the estimation of the volume of hydrocarbons. In this case, a low risk value is assumed, as the random distribution of the peak cells location relative to the control cells locations implies a high level of confidence in the interpolation results.
  • FIG. 7.2 shows an example of a strong clustered extreme values data set resulted from the interpolation of the control cells (720) shown as labeled pipes in FIG. 7.2.
  • the light areas that represent high clustering (density distribution) level exhibit a pronounced localization of the density distribution around the control cells locations.
  • a high risk value is assumed, as the localized distribution of the peak cells location relative to the control cells locations implies a low level of confidence in the interpolation results.
  • FIG. 8 shows a 3D view of an example of the peak cells' distribution in a highly clustered extreme values interpolated data.
  • the locations of the local extrema values marking the peak cells are displayed as dark areas.
  • the input data associated with the control cells come from the wells shown as pipes (810).
  • the peak cells are clustered around the control cells.
  • a correlation coefficient derived based on the dataset presented in FIG. 8 yield a low value (peak cells situated close to control cells relatively to the average distance of volume cells). In this case the interpolation results present a low level of confidence.
  • Embodiments may be implemented on an E&P computing system. Any combination of mobile, desktop, server, router, switch, embedded device, or other types of hardware may be used.
  • the E&P computing system (900) may include one or more computer processors (902), non- persistent storage (904) (e.g., volatile memory, such as random access memory (RAM), cache memory), persistent storage (906) (e.g., a hard disk, an optical drive such as a compact disk (CD) drive or digital versatile disk (DVD) drive, a flash memory, etc.), a communication interface (912) (e.g., Bluetooth interface, infrared interface, network interface, optical interface, etc.), and numerous other elements and functionalities.
  • non- persistent storage e.g., volatile memory, such as random access memory (RAM), cache memory
  • persistent storage e.g., a hard disk, an optical drive such as a compact disk (CD) drive or digital versatile disk (DVD) drive, a flash memory, etc.
  • a communication interface (912)
  • the computer processor(s) (902) may be an integrated circuit for processing instructions.
  • the computer processor(s) may be one or more cores or micro-cores of a processor.
  • the E&P computing system (900) may also include one or more input devices (910), such as a touchscreen, keyboard, mouse, microphone, touchpad, electronic pen, or any other type of input device.
  • the communication interface (912) may include an integrated circuit for connecting the E&P computing system (900) to a network (not shown) (e.g., a local area network (LAN), a wide area network (WAN) such as the Internet, mobile network, or any other type of network) and/or to another device, such as another computing device.
  • a network not shown
  • LAN local area network
  • WAN wide area network
  • the Internet such as the Internet
  • mobile network such as another computing device.
  • the E&P computing system (900) may include one or more output devices (908), such as a screen (e.g., a liquid crystal display (LCD), a plasma display, touchscreen, cathode ray tube (CRT) monitor, projector, or other display device), a printer, external storage, or any other output device.
  • a screen e.g., a liquid crystal display (LCD), a plasma display, touchscreen, cathode ray tube (CRT) monitor, projector, or other display device
  • One or more of the output devices may be the same or different from the input device(s).
  • the input and output device(s) may be locally or remotely connected to the computer processor(s) (902), non-persistent storage (904), and persistent storage (906).
  • the computer processor(s) (902), non-persistent storage (904), and persistent storage (906 may be locally or remotely connected to the computer processor(s) (902), non-persistent storage (904), and persistent storage (906).
  • Software instructions in the form of computer readable program code to perform embodiments may be stored, in whole or in part, temporarily or permanently, on a non-transitory computer readable medium such as a CD, DVD, storage device, a diskette, a tape, flash memory, physical memory, or any other computer readable storage medium.
  • the software instructions may correspond to computer readable program code that, when executed by a processor(s), is configured to perform one or more embodiments.
  • the E&P computing system (900) in FIG. 9.1 may be connected to or be a part of a network.
  • the network (920) may include multiple nodes (e.g., node X (922), node Y (924)).
  • Each node may correspond to an E&P computing system, such as the E&P computing system shown in FIG. 9.1, or a group of nodes combined may correspond to the E&P computing system shown in FIG. 9.1.
  • embodiments may be implemented on a node of a distributed system that is connected to other nodes.
  • embodiments may be implemented on a distributed E&P computing system having multiple nodes, where each portion may be located on a different node within the distributed E&P computing system.
  • one or more elements of the aforementioned E&P computing system (900) may be located at a remote location and connected to the other elements over a network.
  • the node may correspond to a blade in a server chassis that is connected to other nodes via a backplane.
  • the node may correspond to a server in a data center.
  • the node may correspond to a computer processor or micro-core of a computer processor with shared memory and/or resources.
  • the nodes e.g., node X (922), node Y (924) in the network (920) may be configured to provide services for a client device (926).
  • the nodes may be part of a cloud E&P computing system.
  • the nodes may include functionality to receive requests from the client device (926) and transmit responses to the client device (926).
  • the client device (926) may be an E&P computing system, such as the E&P computing system shown in FIG. 9.1. Further, the client device (926) may include and/or perform at least a portion of one or more embodiments.
  • the E&P computing system or group of E&P computing systems described in FIG. 9.1 and 9.2 may include functionality to perform a variety of operations disclosed herein.
  • the E&P computing system(s) may perform communication between processes on the same or different system.
  • a variety of mechanisms, employing some form of active or passive communication, may facilitate the exchange of data between processes on the same device.
  • Examples representative of these inter-process communications include, but are not limited to, the implementation of a file, a signal, a socket, a message queue, a pipeline, a semaphore, shared memory, message passing, and a memory-mapped file. Further details pertaining to a couple of these non-limiting examples are provided below.
  • sockets may serve as interfaces or communication channel end-points enabling bidirectional data transfer between processes on the same device.
  • a server process e.g., a process that provides data
  • the server process may create a first socket object.
  • the server process binds the first socket object, thereby associating the first socket object with a unique name and/or address.
  • the server process then waits and listens for incoming connection requests from one or more client processes (e.g., processes that seek data).
  • client processes e.g., processes that seek data.
  • the client process then proceeds to generate a connection request that includes at least the second socket object and the unique name and/or address associated with the first socket object.
  • the client process then transmits the connection request to the server process.
  • the server process may accept the connection request, establishing a communication channel with the client process, or the server process, busy in handling other operations, may queue the connection request in a buffer until server process is ready.
  • An established connection informs the client process that communications may commence.
  • the client process may generate a data request specifying the data that the client process wishes to obtain.
  • the data request is subsequently transmitted to the server process.
  • the server process analyzes the request and gathers the requested data.
  • the server process then generates a reply including at least the requested data and transmits the reply to the client process.
  • the data may be transferred, more commonly, as datagrams or a stream of characters (e.g., bytes).
  • Shared memory refers to the allocation of virtual memory space in order to substantiate a mechanism for which data may be communicated and/or accessed by multiple processes.
  • an initializing process first creates a shareable segment in persistent or non-persistent storage. Post creation, the initializing process then mounts the shareable segment, subsequently mapping the shareable segment into the address space associated with the initializing process. Following the mounting, the initializing process proceeds to identify and grant access permission to one or more authorized processes that may also write and read data to and from the shareable segment. Changes made to the data in the shareable segment by one process may immediately affect other processes, which are also linked to the shareable segment. Further, when one of the authorized processes accesses the shareable segment, the shareable segment maps to the address space of that authorized process. Often, one authorized process may mount the shareable segment, other than the initializing process, at any given time.
  • the E&P computing system performing one or more embodiments may include functionality to receive data from a user.
  • a user may submit data via a GUI on the user device.
  • Data may be submitted via the GUI by a user selecting one or more GUI widgets or inserting text and other data into GUI widgets using a touchpad, a keyboard, a mouse, or any other input device.
  • information regarding the particular item may be obtained from persistent or non-persistent storage by the computer processor.
  • the contents of the obtained data regarding the particular item may be displayed on the user device in response to the user's selection.
  • a request to obtain data regarding the particular item may be sent to a server operatively connected to the user device through a network.
  • the user may select a uniform resource locator (URL) link within a web client of the user device, thereby initiating a Hypertext Transfer Protocol (HTTP) or other protocol request being sent to the network host associated with the URL.
  • HTTP Hypertext Transfer Protocol
  • the server may extract the data regarding the particular selected item and send the data to the device that initiated the request.
  • the contents of the received data regarding the particular item may be displayed on the user device in response to the user's selection.
  • the data received from the server after selecting the URL link may provide a web page in Hyper Text Markup Language (HTML) that may be rendered by the web client and displayed on the user device.
  • HTML Hyper Text Markup Language
  • the E&P computing system may extract one or more data items from the obtained data.
  • the extraction may be nerformed as follows by the E&P computing system in FIG. 9.1.
  • the organizing pattern e.g., grammar, schema, layout
  • the organizing pattern is determined, which may be based on one or more of the following: position (e.g., bit or column position, Nth token in a data stream, etc.), attribute (where the attribute is associated with one or more values), or a hierarchical/tree structure (consisting of layers of nodes at different levels of detail-such as in nested packet headers or nested document sections).
  • the raw, unprocessed stream of data symbols is parsed, in the context of the organizing pattern, into a stream (or layered structure) of tokens (where each token may have an associated token "type").
  • extraction criteria are used to extract one or more data items from the token stream or structure, where the extraction criteria are processed according to the organizing pattern to extract one or more tokens (or nodes from a layered structure).
  • the token(s) at the position(s) identified by the extraction criteria are extracted.
  • the token(s) and/or node(s) associated with the attribute(s) satisfying the extraction criteria are extracted.
  • the token(s) associated with the node(s) matching the extraction criteria are extracted.
  • the extraction criteria may be as simple as an identifier string or may be a query presented to a structured data repository (where the data repository may be organized according to a database schema or data format, such as XML).
  • the extracted data may be used for further processing by the E&P computing system.
  • the E&P computing system of FIG. 9.1 while performing one or more embodiments, may perform data comparison.
  • the comparison may be performed by submitting A, B, and an opcode specifying an operation related to the comparison into an arithmetic logic unit (ALU) (i.e.. circuitry that performs arithmetic and/or bitwise logical operations on the two data values).
  • ALU arithmetic logic unit
  • the ALU outputs the numerical result of the operation and/or one or more status flags related to the numerical result.
  • the status flags may indicate whether the numerical result is a positive number, a negative number, zero, etc.
  • B may be subtracted from A (i.e., A - B), and the status flags may be read to determine if the result is positive (i.e., if A > B, then A - B > 0).
  • a and B may be vectors, and comparing A with B includes comparing the first element of vector A with the first element of vector B, the second element of vector A with the second element of vector B, etc.
  • if A and B are strings, the binary values of the strings may be compared.
  • the E&P computing system in FIG. 9.1 may implement and/or be connected to a data repository.
  • a data repository is a database.
  • a database is a collection of information configured for ease of data retrieval, modification, re-organization, and deletion.
  • Database Management System is a software application that provides an interface for users to define, create, query, update, or administer databases.
  • the user, or software application may submit a statement or query into the
  • the DBMS interprets the statement.
  • the statement may be a select statement to request information, update statement, create statement, delete statement, etc.
  • the statement may include parameters that specify data, or data container (database, table, record, column, view, etc.), identifier(s), conditions (comparison operators), functions (e.g. join, full join, count, average, etc.), sort (e.g. ascending, descending), or others.
  • the DBMS may execute the statement. For exarrmle. the DBMS may access a memory buffer, a reference or index a file for read, write, deletion, or any combination thereof, for responding to the statement.
  • the DBMS may load the data from persistent or non-persistent storage and perform computations to respond to the query.
  • the DBMS may return the result(s) to the user or software application.
  • the E&P computing system of FIG. 9.1 may include functionality to present raw and/or processed data, such as results of comparisons and other processing.
  • presenting data may be accomplished through various presenting methods.
  • data may be presented through a user interface provided by a computing device.
  • the user interface may include a GUI that displays information on a display device, such as a computer monitor or a touchscreen on a handheld computer device.
  • the GUI may include various GUI widgets that organize what data is shown as well as how data is presented to a user.
  • the GUI may present data directly to the user, e.g., data presented as actual data values through text, or rendered by the computing device into a visual representation of the data, such as through visualizing a data model.
  • a GUI may first obtain a notification from a software application requesting that a particular data object be presented within the GUI.
  • the GUI may determine a data object type associated with the particular data object, e.g., by obtaining data from a data attribute within the data object that identifies the data object type.
  • the GUI may determine any rules designated for displaying that data object type, e.g., rules specified by a software framework for a data object class or according to any local parameters defined by the GUI for presenting that data object type.
  • the GUI may obtain data values from the particular data object and render a visual representation of the data values within a display device according to the designated rules for that data object type.
  • Data may also be presented through various audio methods. In particular, data may be rendered into an audio format and presented as sound through one or more speakers operably connected to a computing device.
  • Data may also be presented to a user through haptic methods.
  • haptic methods may include vibrations or other physical signals generated by the E&P computing system.
  • data may be presented to a user using a vibration generated by a handheld computer device with a predefined duration and intensity of the vibration to communicate the data.
  • embodiments may be utilized in conjunction with a handheld system (i.e., a phone, wrist or forearm mounted computer, tablet, or other handheld device), portable system (i.e., a laptop or portable computing system), a fixed computing system (i.e., a desktop, server, cluster, or high performance computing system), or across a network (i.e., a cloud-based system).
  • a handheld system i.e., a phone, wrist or forearm mounted computer, tablet, or other handheld device
  • portable system i.e., a laptop or portable computing system
  • a fixed computing system i.e., a desktop, server, cluster, or high performance computing system
  • a network i.e., a cloud-based system

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

Selon l'invention, un volume en subsurface interpolé comprenant des mailles témoins et des mailles déduites est obtenu, les mailles déduites comprenant des valeurs d'attribut déduites des mailles témoins. Un ensemble de mailles de crête est choisi en fonction de valeurs d'attribut. La valeur d'attribut provenant des mailles déduites ou des mailles témoins dans la maille de crête est un extremum local. Un coefficient de corrélation est déterminé, lequel coefficient caractérise une relation spatiale entre l'emplacement des mailles de crête et l'emplacement des mailles témoins. Le coefficient de corrélation est présenté.
PCT/US2016/050842 2015-09-24 2016-09-09 Évaluation de volume en subsurface WO2017053080A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201562222820P 2015-09-24 2015-09-24
US62/222,820 2015-09-24

Publications (1)

Publication Number Publication Date
WO2017053080A1 true WO2017053080A1 (fr) 2017-03-30

Family

ID=58387005

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2016/050842 WO2017053080A1 (fr) 2015-09-24 2016-09-09 Évaluation de volume en subsurface

Country Status (1)

Country Link
WO (1) WO2017053080A1 (fr)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6016462A (en) * 1997-08-29 2000-01-18 Exxon Production Research Company Analysis of statistical attributes for parameter estimation
US6757614B2 (en) * 2000-12-18 2004-06-29 Schlumberger Technology Corporation Seismic signal processing method and apparatus for generating correlation spectral volumes to determine geologic features
US20050171700A1 (en) * 2004-01-30 2005-08-04 Chroma Energy, Inc. Device and system for calculating 3D seismic classification features and process for geoprospecting material seams
US20120090834A1 (en) * 2009-07-06 2012-04-19 Matthias Imhof Method For Seismic Interpretation Using Seismic Texture Attributes
US20130020131A1 (en) * 2011-06-09 2013-01-24 Mickaele Le Ravalec Method of developing a petroleum reservoir from a technique for selecting the positions of the wells to be drilled

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6016462A (en) * 1997-08-29 2000-01-18 Exxon Production Research Company Analysis of statistical attributes for parameter estimation
US6757614B2 (en) * 2000-12-18 2004-06-29 Schlumberger Technology Corporation Seismic signal processing method and apparatus for generating correlation spectral volumes to determine geologic features
US20050171700A1 (en) * 2004-01-30 2005-08-04 Chroma Energy, Inc. Device and system for calculating 3D seismic classification features and process for geoprospecting material seams
US20120090834A1 (en) * 2009-07-06 2012-04-19 Matthias Imhof Method For Seismic Interpretation Using Seismic Texture Attributes
US20130020131A1 (en) * 2011-06-09 2013-01-24 Mickaele Le Ravalec Method of developing a petroleum reservoir from a technique for selecting the positions of the wells to be drilled

Similar Documents

Publication Publication Date Title
US11775858B2 (en) Runtime parameter selection in simulations
US11269110B2 (en) Computing system assessment of geological similarity of wells employing well-log data
US20210125312A1 (en) Artificial intelligence technique to fill missing well data
US11782184B2 (en) Enhanced two point flux approximation scheme for reservoir simulation
EP3469498A1 (fr) Étalonnage automatique de la modélisation d'un champ
WO2017024113A1 (fr) Procédé de l'évaluation des propriétés de transport de fluide dans une formation géologique hétérogène
NO347918B1 (en) Efficient algorithms for volume visualization on irregular grids
WO2017058738A1 (fr) Déroulement des opérations de simulation basée sur un réseau
EP3398090B1 (fr) Imagerie géologique et inversion à l'aide d'un dispositif de stockage d'objet
EP3469404B1 (fr) Segmentation d'un volume structural
EP4104077A1 (fr) Reconnaissance automatique d'activités de forage sur la base de codes opérationnels rapportés quotidiennement
WO2017053080A1 (fr) Évaluation de volume en subsurface
EP3475734B1 (fr) Mise en oeuvre d'advection libre en modélisation de bassin
US11422874B2 (en) Visualization infrastructure for web applications
EP3510425B1 (fr) Calcul de zone d'infiltration de puits à l'aide de données de diagraphie en cours de forage
WO2016137888A1 (fr) Visualisation d'ensembles de données

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16849338

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16849338

Country of ref document: EP

Kind code of ref document: A1