Patent EA-003475-B1 (EA)
Title: Process for constructing three-dimensional geologic models having adjustable geologic interfaces
Inventor: CALVERT CRAIG S
IPC: G01V1/28
Filed: 2000-10-12
Abstract: 1. A process for constructing a three-dimensional geologic model of a subsurface earth volume containing one or more geologic interfaces, said process comprising the steps of: (a) generating a tentative geologic model of said subsurface earth volume, said tentative geologic model comprising a three-dimensional array of contiguous model blocks, each model block having tentative values for one or more rock properties assigned thereto, said tentative geologic model containing tentative positions for said one or more geologic interfaces; (b) specifying training criteria which define the spatial attributes of geologic interfaces and the characteristics of rock properties in said subsurface earth volume; (c) calculating statistics that describe the spatial attributes of said tentative geologic-interface positions and the characteristics of said tentative rock-property values in said tentative geologic mode; and (d) comparing statistics calculated in step (c) with corresponding training criteria specified in step (b) and either (1) if said statistics do not match said training criteria within specified limits, perturbing the tentative position within said tentative geologic model of at least a portion of at least one of said geologic interfaces, updating said tentative geologic model, and repeating steps (c) and (d), or (2) if said statistics match said training criteria within said specified limits, accepting said tentative geologic model as the three-dimensional geologic model for said subsurface earth volume. 2. The process of claim 1, wherein each geologic interface in said subsurface earth volume is represented in said tentative geologic model by a set of interface nodes, and wherein said perturbation in step (d) (1) comprises changing the position of only one of said interface nodes. 3. The process of claim 1, wherein each geologic interface in said subsurface earth volume is represented in said tentative geologic model by a set of interface nodes, and wherein said perturbation in step (d) (1) comprises changing the positions of all interface nodes within a specified distance from a specific interface node. 4. The process of claim 1, wherein step (d) (1) further includes perturbing the tentative value of at least one rock property for at least one of said model blocks prior to updating said tentative geologic model. 5. The process of claim 4, said method further comprising the step of moving the tentative positions within said tentative geologic model of said geologic interfaces as required to be consistent with said perturbed rock property values. 6. The process of claim 1, wherein said training criteria include the match between synthetic seismic traces derived from said tentative geologic model and actual seismic data traces for said subsurface earth volume. 7. The process of claim 1, wherein one or more wells are located in said subsurface earth volume, and wherein actual geologic-interface locations obtained from said wells are inserted into said tentative geologic model at fixed positions corresponding to their locations in said wells, said fixed geologic-interface positions being accepted as correct and not subject to perturbation. 8. The process of claim 4, wherein one or more wells are located in said subsurface earth volume, and wherein actual rock-property values obtained from said one or more wells are assigned to model blocks representing corresponding portions of said subsurface earth volume, said actual rock-property values being accepted as correct and not subject to perturbation. 9. The process of claim 1, wherein the vertical dimension of said tentative geologic model is seismic time. 10. The process of claim 1, wherein the vertical dimension of said tentative geologic model is vertical depth. 11. The process of claim 1, wherein separate training criteria are specified for each zone and each region in said subsurface earth volume. 12. The process of claim 1, wherein said array of contiguous model blocks is replace by an array of discrete model points, and wherein said tentative rock-property values are assigned to said model points. 13. The process of claim 12, wherein at least a portion of said discrete model points correspond to sample points in a three-dimensional seismic data volume for said subsurface earth volume. 14. The process of claim 1, wherein said process is performed using a suitably programmed digital computer. 15. A process for constructing a three-dimensional geologic model of a subsurface earth volume containing one or more geologic interfaces, said process comprising the steps of: (a) generating an initial geologic model of said subsurface earth volume by (1) specifying initial positions within said initial geologic model for said one or more geologic interfaces, (2) defining a three-dimensional array of contiguous model blocks to represent said subsurface earth volume, and (3) assigning initial rock-property values to each of said model blocks; (b) specifying training criteria that define the characteristics of rock properties in said subsurface earth volume and the spatial attributes of geologic interfaces in said subsurface earth volume, and specifying desired values for said training criteria; (c) calculating statistics comparing the characteristics of rock properties in said initial geologic model and the spatial attributes of geologic interfaces in said initial geologic model to the desired values for corresponding training criteria; (d) combining the statistics calculated in step (c) into an objective function and calculating an initial value for said objective function; (e) randomly selecting at least one model block and tentatively perturbing the value of at least one rock property thereof; (f) if said selected model block is within a preselected distance from a geologic interface, determining whether or not to tentatively shift the position of said geologic interface to the opposite side of said selected model block; (g) calculating tentative statistics and a tentative value for said objective function based on said rock-property perturbation and said geologic-interface position shift, if any; (h) determining whether to accept or reject said rock-property perturbation and said geologic-interface position shift, if any, based on a preselected optimization procedure; (i) if either or both of said rock-property perturbation and said geologic-interface position shift are accepted, updating said geologic model and said objective-function value ; and (j) repeating steps (e) through (i) until a preselected number of perturbations have been attempted or predetermined criteria for said objective function have been attained. 16. The process of claim 15, wherein the determination in step (f) of whether or not to tentatively shift the position of said geologic interface to the opposite side of said selected model block is made randomly. 17. The process of claim 15, wherein the determination in step (f) of whether or not to tentatively shift the position of said geologic interface to the opposite side of said selected model block is made by calculating statistics to measure the effects of moving and of not moving the position of said geologic interface and selecting the option that best matches the desired values of the training criteria. 18. The process of claim 15, wherein said training criteria include the match between synthetic seismic traces derived from said tentative geologic model and actual seismic data traces for said subsurface earth volume. 19. The process of claim 15, wherein one or more wells are located in said subsurface earth volume, and wherein actual rock-property values and actual geologic-interface positions obtained from said wells are inserted into said initial geologic model at corresponding model positions, said actual rock property values and said actual geologic-interface positions being accepted as correct and not subject to perturbation. 20. The process of claim 19, wherein initial rock-property values for all model blocks that do not correspond to wells are determined at random. 21. The process of claim 15, wherein said preselected optimization procedure is simulated annealing. 22. The process of claim 15, wherein said preselected optimization procedure is the maximum a posteriori optimization procedure.