Method of modelling the production of an oil reservoir

[0001] 1. Field of the Invention

[0002] The present invention relates to the study and to the optimization of oil reservoir production schemes and models of production behavior of an oil reservoir in order to compare production schemes and to define an optimum scheme considering a given production criterion such as oil recovery, water inflow, production rate, etc.

[0003] 2. Description of the Prior Art

[0004] The study of a reservoir comprises two main stages.

[0005] The reservoir characterization stage determines a numerical flow model or flow simulator that is compatible with the real data collected in the field. Engineers have access to only a small part of the reservoir under study (core analysis, logging, well tests) and have to extrapolate these limited data over the entire oilfield to construct the numerical simulation model.

[0006] The production prediction stage uses the numerical simulation model to estimate the reserves and production to be obtained in the future or to improve the production scheme in place. This stage is carried out using a numerical simulation model constructed from many data sources, obtained from only a small part of the reservoir. Consequently, an uncertainty notion has to be constantly accounted for in this stage.

[0007] In order to properly characterize the impact of each uncertainty on the oil production, the largest possible number of production scenarios has to be tested, which therefore requires a large number of reservoir simulations. Considering the long time required for a flow simulation, it is clearly not conceivable to test all possible scenarios via the numerical flow model. In this context, using the experimental design method can allow construction of a simplified model of the flow simulator as a function of a reduced number of parameters. Experimental designs allow determination of the number and the location in space of the parameters of the simulations to be carried out so as to have a maximum amount of pertinent data at the lowest possible cost. This simple model translates the behavior of a given response (for example the 10-year cumulative oil production) as a function of some parameters. Its construction requires a reduced number of simulations previously defined by experimental designs.

[0008] During the production prediction stage, the simplified model is used because it is simple and analytical and, therefore, each simulation obtained by this model is immediate. This saves considerable time. Using this model allows the reservoir engineer to test as many scenarios as are wanted, without having to care about the time required to perform a numerical flow simulation.

[0009] The methods disclosed in French patents 2,855,631 and 2,855,633 use simplified models to optimize the production of an oil reservoir or as decision support for managing an oil reservoir, in the presence of uncertainties.

[0010] The simplified model obtained using experimental designs implies that the response obtained by the model is a linear function of the parameters taken into account. However, in most cases, this is not true. When the range of evolution of a parameter(s) (permeability, porosity, . . . ) is relatively limited and its contribution is reasonable, its behavior can be assumed to be linear. But when this range of evolution of the parameter(s) becomes too large or when the contribution of the parameter is no longer linear, the linearity hypothesis biases the knowledge of the oil reservoir.

[0011] It is therefore necessary to set a criterion allowing detection of non-linearities and to establish an efficient and fast methodology allowing prediction, in an effective manner, of non-linear response behaviors.

[0012] The present invention models an oil reservoir by iterative adjustments so as to best reproduce the behavior of the oil reservoir, while controlling the number of simulations.

[0013] In general terms, the present invention relates to a method for simulating production of an oil reservoir wherein the following stages are carried out:

[0014] a) providing a flow simulator utilizing physical data obtained from measurements in the oil reservoir;

[0015] b) determining a first analytical model representing the production of the reservoir as a function of time by accounting for parameters which influence production of the reservoir, which provides adjustment to production values provided by the flow simulator;

[0016] c) selecting at least one new production value associated with a point within an area of the reservoir selected as a function of non-linearity of production in the area of the reservoir with the new value being obtained from the flow simulator; and

[0017] d) determining a second model by adjusting the first model so that a response of the second model at the point corresponding to the new production value.

[0018] According to the invention, in c), the following steps can be carried out:

[0019] determining a sub-model that optimally adjusts to production values, except for a test value selected from among the production values;

[0020] calculating a prediction residue associated with the test value by determining a difference between a response of the sub-model and the test value;

[0021] calculating the prediction residue associated with each production value by repeating determining a sub-model and calculating a prediction residue associated with the test value by assigning successively to the test value each value of the production values; and

[0022] selecting the new production value in another area of the reservoir in a vicinity of the point with the new production value having a greatest prediction residue.

[0023] The new production value can be selected accounting for a gradient of production at a point associated with the production value having the greatest prediction residue.

[0024] Furthermore, a new value can be selected in c) and d) can be carried out provided that the greatest prediction residue is greater than a previously set value.

[0025] According to a variant of the invention, in step c), the following steps can be carried out:

[0026] determining a first kriging variance of the first model for the production values obtained from the flow simulator;

[0027] selecting a first pilot point in the reservoir at a location where the first kriging variance is maximum;

[0028] determining a second kriging variance of the first model for the production values obtained from the flow simulator and the first pilot point;

[0029] selecting a second pilot point in the reservoir at a location where the second kriging variance is maximum; and

[0030] assigning a value to each pilot point by carrying out the following steps for each pilot point:

[0031] determining a sub-model that adjusts the production values and to a value associated with one of the pilot points, except for a test value selected from among production values and the value associated with the one pilot point;

[0032] calculating a prediction residue associated with the test value determining the difference between a response of the sub-model and the test value;

[0033] calculating the prediction residue associated with each sub-model response by repeating determining the sub-model and calculating the prediction residue by assigning successively to the test value each value contained in a set of production values and a value associated with the one pilot point;

[0034] calculating a sum of absolute values of the prediction residues calculated for each test value;

[0035] assigning to the one pilot point a value that minimizes this sum;

[0036] determining a second sub-model that adjusts to the production values and values of the pilot points;

[0037] for each pilot point, determining a difference between a response of the second sub-model and a response of the first model; and

[0038] associating the new production value of step c) with the pilot point for which the difference is the greatest.

[0039] Furthermore, in d), the second model can be determined by adjusting the first model so that the response of the second model at the pilot point selected corresponds to the new production value and, furthermore, to the values assigned to other pilot points.

[0040] According to another variant of the invention, in c), the following steps can be carried out:

[0041] determining an analytical model expressing a derivative of the reservoir production as a function of time, the model adjusting to derivatives at points associated with the production values used in b); and

[0042] from the model expressing the derivative, selecting at least one new production value associated with a point whose response of the model expressing the derivative is zero.

[0043] It is possible to select a new value in c) and d) can be carried out, provided that the prediction residue of the new selected value is greater than a previously set value.

[0044] According to the invention, after d), the following steps are carried out:

[0045] determining a third analytical model expressing the derivative of the reservoir production as a function of time, the third model adjusting to the derivatives at the points associated with the production values and the production values selected in c);

[0046] if the response of the third analytical model at the point selected in c) is greater than zero, determining a point associated with the maximum value of the response of the second model in a vicinity of the point selected in c);

[0047] if the response of the third analytical model at the point selected in c) is less than zero, determining a point associated with a minimum value of the response of the second model in a vicinity of the point selected in c),

[0048] determining a new production value using the flow simulator at a point associated with a previously determined minimum or maximum value;

[0049] determining a fourth model by adjusting the second model so that the response of the fourth model corresponds to a new value determined in the determining a new production value using the flow simulator.

[0050] According to the invention, c) and d) can be repeated.

[0051] In b), the production values can be selected using an experimental design.

[0052] In b), the first model can be adjusted using one of the following approximation methods: polynomial approximation, neural networks or support vector machines.

[0053] In d), one of the following interpolation methods can be used: kriging or spline methods.

[0054] Thus, the method according to the invention provides the reservoir engineer with a simple and inexpensive form of numerical simulation for scenario management and production scheme optimization, as a support for decision-making for minimizing risks.

[0059] The method according to the invention is illustrated by the diagram of FIG. 1.

[0060] Step 1: Providing a Reservoir Flow Simulator

[0061] The oil reservoir is modelled using a numerical reservoir simulator. The reservoir simulator or flow simulator notably allows calculation of production of hydrocarbons or of water in time as a function of technical parameters such as a number of layers in the reservoir, permeability of the layers, aquifer force, position of oil wells, etc. Furthermore, the flow simulator calculates the derivative of the production value at the point which is considered.

[0062] The numerical simulator is provided using characteristic data of the oil reservoir. For example, the data are obtained by measurements performed in the laboratory on cores and fluids taken from the oil reservoir, by logging, well tests, etc.

[0063] Step 2: Approximation of the Flow Simulator

[0064] Since the flow simulator is complex and calculation time is consuming, a simplified model of the behavior of the oil reservoir is constructed.

[0065] Parameters having an influence on the hydrocarbon or water production profiles of the reservoir are selected. Selection of the parameters can be done either through physical knowledge of the oil reservoir, or by means of a sensitivity analysis. For example, it is possible to use a statistical Student or Fischer test.

[0066] Some parameters can be intrinsic to the oil reservoir. For example, the following parameters can be considered: a permeability multiplier for particular reservoir layers, aquifer force of residual oil saturation after waterflooding.

[0067] Some parameters can correspond to reservoir development options. These parameters can be well position, completion level or drilling technique.

[0068] Points for which the numerical flow simulations will be performed are selected from the experimental domain. These points are used to provide a simplified model that best reproduces the reservoir flow simulator. These points are selected using an experimental design method, which allows determination of the number and the location of the simulations to be carried out so as to have a maximum amount of information at the lowest possible cost, and thus to determine a reliable model best representing the production profile. It can be noted that selection of the experimental design method is very important: the initial experimental design method plays an essential part in determination of the modelling of the first model, and the results greatly depend on a pattern of experimentation.

[0069] Selection of the simulation points is determined using experimental designs, for example factorial designs, composite designs, Latin hypercubes, maximum distance designs, etc. It is possible to use the experimental designs described in the following documents:

1. Dejean, J. P. and Blanc, G., “Managing Uncertainties on Production Predictions Using Integrated Statistical Methods”, SPE 56696, SPE Annual Technical Conference and Exhibition, Houston, USA, Oct. 3-6, 1999.

2. Box, G. E. P. and Hunter, J. S., “The 2k-p Fractional Factorial Designs”, Part I, Technometrics, 2, 311-352, 1961a

3. Box, G. E. P. and Hunter, J. S., “The 2k-p Fractional Factorial Designs”, Part II, Technometrics, 3, 449-458, 1961b

4. Box, G. E. P and Wilson, K. B., “On the Experimental Attainment of Optimum Conditions”, Journal of the Royal Statistical Society, Series B, 13, 1-45

5. Draper, N. R., “Small Composite Designs”, Technometrics, 27, 173-180, 1985

6. Atkinson, A. C. and Donev, A. N., “Optimum Experimental Designs”, Oxford University press, 1992.

[0070] After determination of a first experimental design and when the numerical simulations are performed, an approximation method is used to provide a first model representing a trend of behaviour of the response function, that is which approximates a flow simulator.

[0071] The first model expresses a production criterion studied over the course of time. The production criterion is expressed as a function of the selected parameters. The production criterion can be for oil recovery, water inflow or a rate of production. The first analytical model is determined using previously selected values of the production criterion obtained from the flow simulator.

[0072] When referring to approximation methods, consideration is given to polynomials of the first or second order, neural networks, support vector machines or possibly polynomials of an order greater than two. Selection of this model depends on one hand on a maximum number of simulations that can be utilized by the user and, on the other hand, on the initial experimental design that is used.

[0073] Step 3: Adjustment of the First Model

[0074] There may be a difference between the production value given by the first analytical model obtained in step 2 and the simulated production values used to provide the first model.

[0075] In this case, the residues are determined at the various simulation points. The residues correspond to the difference between a response of the first model and the value obtained by the reservoir flow simulator. Then, the residues are interpolated. Any n-dimensional interpolation method is suitable. The kriging or the spline method can be used in particular. These methods are explained in the book entitled “Statistics for Spatial Data” by Cressie, N., Wiley, New York 1991.

[0076] The residue interpolation structure lends itself well to this sequential approach because it is divided up into two parts: a linear model, which corresponds to the first model determined in step 2, and a “correcting” term allowing making up the difference between the prediction of the first model and the simulation point. In cases where the analytical model should be satisfactory, it is not necessary to add this “correcting” term. In the opposite case, it allows interpolation of the responses and, thus, accounting for non-linearities detected at the surface.

[0077] An adjusted second model is thus provided by adding the results of the interpolations of the residues to the first model determined in step 2.

[0078] Step 4: Model Predictivity Test and Selection of Additional Simulation Points

[0079] At this step of the modelling procedure, the second model interpolates exactly the simulations, therefore adjustment of the response function is optimum. Considering that the interpolation method is exact, the “conventional” residues are zero. Therefore, according to the invention, an interest is taken in the prediction residues. Therefore, the predictivity of the model is examined for the points outside the experimental design. The predictions have to be as accurate as possible. Consequently, a model predictivity test is carried out to evaluate the approximation quality so as to judge whether an improvement is necessary by addition of new points to the initial design.

[0080] Two criteria are involved in the predictivity test:

[0081] (1) a priori predictivity calculation with prediction residues calculation

[0082] (2) a posteriori predictivity calculation with use of confirmation points.

[0083] A Priori Predictivity

[0084] The prediction residues are the residues obtained at a point of the design by carrying out adjustment of the first model without this point. Removing a point and re-estimating the model will allow determination of whether this point (or the zone of the design close to this point) provides decisive information or not. Calculation of these prediction residues is carried out for each point of the initial experimental design. In the vicinity of the points considered the least predictive of the current design, that is the points having the greatest prediction residue, new points are simulated. A sub-sampling zone is therefore defined in the vicinity of the points. Addition of these points can be conditioned by the fact that the residues are greater than a value set by the user.

[0085] The size of this sub-sampling zone can be defined using the information on the gradients of the production at the points and/or the value of the prediction residues. In fact, a high gradient value expresses a high variation of the response. It can therefore be informative to add a new point close to the existing one. On the other hand, a low gradient value in a given direction shows that there are no irregularities in this direction. It is therefore not necessary to investigate a wide variation range in this direction. To the contrary, the variation range for one of the parameters is all the wider as the value of the gradient is high in this direction. This approach allows elimination of certain directions (where the value of the gradient is not significant) and thus to reduce the number of simulations to be performed. This sub-sampling can for example result from the construction of a new experimental design defined in this zone. Selection of this experimental design (factorial design, composite design or Latin hypercube) results from the necessary compromise between the modelling cost and quality.

[0086] Alternatively, the pilot point method can be used to improve the second model.

[0087] For a given number of experimentations, there is a large number of estimators (exact interpolators) going through all the experimentations and respecting the spatial structure (expectation and covariance) of the process. In this class of estimators regarding the data, the estimation is sought that maximizes the a priori predictivity. In order to go through this class of estimators, fictitious information is added, that is, pilot points are added to the simulated experimentations. These pilot points are then considered to be data although no simulation has been carried out and allow going through all the estimators passing through all the experimentations. The goal is to select the interpolator that maximizes the a priori predictivity coefficient of the model, that is, the pilot points are positioned so as to obtain the maximum predictivity realization.

[0088] The location of a pilot point is determined by accounting for the following two criteria:

[0089] (1) the capacity of the pilot point to reduce the difference between the observations and the results of numerical flow simulations; and

[0090] (2) the contribution of the pilot point to the reduction of the uncertainties on the current approximation model.

[0091] For this selection to be made in an optimum way, the impact of a possible pilot point on each one of these two criteria has to be quantified.

[0092] In order to remove the prediction uncertainty regarding places which have small representation, it is interesting to apply local perturbations to the zones with a high kriging variance (absence of observations). A pilot point is thus placed where the kriging variance is maximum. Methods for determining the kriging variance are described in the book entitled “Statistics for Spatial Data” by Cressie, N., Wiley, New York 1991.

[0093] The following operations are carried out to determine the location of a pilot point:

[0094] determining the kriging variance in the uncertain domain of the second model determined in step 3 for the finite number of production values obtained by the flow simulator,

[0095] placing a first pilot point where the kriging variance is maximum.

[0096] It is assumed that, besides the production values obtained by the flow simulator, a certain number of pilot points has already been positioned in the uncertain domain and new pilot points are to be positioned to improve the model predictivity. The existing pilot points are then considered as local data of zero variance. By taking account of the location of already existing points, optimizing of the location of the pilot points sequentially occurs.

[0097] Thus, to determine the location of a second pilot point, the following operations are carried out:

[0098] (1) determining the kriging variance of the first model for the finite number of production values obtained by the flow simulator and the first pilot point;

[0099] (2) determining the location of a second pilot point where the kriging variance is maximum.

[0100] Several pilot points can be added by repeating the previous two operations.

[0101] It is preferable to add a number of pilot points that is less than or equal to the number of real experiments so as not to perturb the model. Once the optimum location of the pilot points is determined, a “fictitious” response value has to be assigned at these points.

[0102] Since the goal of the addition of pilot points is to improve the a priori predictivity of the model, the value of the pilot points have to be defined from an objective function that measures this predictivity. Since kriging is an exact interpolation method, the “conventional” residues are zero. These residues therefore provide no information on the predictivity and consequently the prediction residues are considered. What is referred to as a priori predictivity is the calculation of the prediction residues at each point of the initial experimental design. The prediction residues are the residues obtained at a point of the initial experimental design by adjusting the first model without this point.

[0103] The following steps can be carried out to determine the production value associated with one of the pilot points whose location has been previously determined:

[0104] (a) determining a sub-model that adjusts to a finite number of production values and to a value associated with the pilot point, except for a test value selected from among the finite number of production values and the value associated with the pilot point;

[0105] (b) calculating a prediction residue associated with the test value by determining a difference between the sub-model response and the test value;

[0106] (c) calculating the prediction residue associated with each response of the prediction sub-model by repeating (a) and (b) by assigning successively to the test value each one of the values contained in the finite number of production values and the value associated with the pilot point;

[0107] (d) calculating a sum of the absolute values or of squares of the prediction residues determined for each test value; and

[0108] (e) assigning to the pilot point a value that minimizes this sum.

[0109] Removing a point and re-estimating the model allows determination of whether this point or the zone of the experimental domain close to this point provides decisive information or does not. Calculation of the prediction residues is carried out in a vicinity of the pilot point to be optimized. Initial values for the pilot points are set, then these data are considered as real and the value of the pilot point is varied to obtain a model that is as predictive as possible, that is, it is desired to minimize the mean prediction error of the model.

[0110] Determination of the optimum value of the pilot point is thus performed to minimize a mean prediction error of the model throughout the uncertain domain. Similarly, this determination of the optimum value of the pilot point can be carried out so as to minimize local prediction error of the model (that is in a vicinity of the pilot point, regardless of the other prediction errors).

[0111] Once the value and the position of the pilot points are determined, testing of the sensitivity of the model occurs at the new points which have been added and then simulations are carried out at the points that are very sensitive to change in the approximation. The estimator obtained without pilot points is compared with the estimator obtained by kriging with pilot points (that is the maximum predictivity realization).

[0112] The points exhibiting the greatest disagreement, that is with the greatest difference, translate to a high approximation instability. Consequently, it is essential to improve the approximation quality in these places. Thus, the simulations corresponding to the points with the greatest disagreement are carried out in order to stabilize the approximation.

[0113] In order to select the pilot points for which a simulation will be carried out, the following stages can be carried out:

[0114] determining a sub-model from the pilot points and the finite number of production values; and

[0115] for each pilot point, calculating a difference between the response of this sub-model and the response of the second model determined in step 3,

According to a first variant:

[0116] Selecting the pilot point for which the difference between the response of the sub-model and the response of the second model is the greatest. It is the point selected for improving the first model, the other pilot points are then ignored in the rest of the procedure.

According to a second variant:

[0117] Selecting one or more pilot points for which the predictivity is the poorest (less than a threshold below 1) since this low predictivity expresses a high sensitivity of the point. In the rest of the procedure, it is taken into account, on the one hand, the production values associated with the selected pilot points with production values being obtained by the flow simulator, and, on the other hand, the production values associated with the other pilot points whose predictivity is better with production values corresponding to the values estimated according to the aforementioned a priori predictivity.

[0118] According to the second variant, if the procedure is repeated, the local predictivity at non-simulated pilot points then has to be evaluated again to ensure that this value still corresponds to a satisfactory stabilization. If this is not the case, the non-simulated pilot point is no longer considered in the new estimation.

[0119] Addition of these new simulations then allows the residues to be studied. What is referred to as residues here is, for each pilot point, the difference between the simulated value and the value obtained upon optimization of the pilot points.

[0120] As before, if the residues are too great, there is a disagreement between the current approximation with the pilot points and the simulations; this expresses a predictivity defect of the model. In this case, the current model has to be improved, which again requires new simulations. One or more new iterations therefore have to be carried out.

[0121] On the other hand, if the residues are small, the prediction at these points is good and therefore the model seems to be predictive in the considered domains. The global predictivity of the model however needs to be confirmed by adding confirmation points. These new simulations allow determination of whether the iteration procedure has to be continued or not.

[0122] A Posteriori Predictivity

[0123] It is possible to add confirmation points, which are production values obtained by the flow simulator determined in step 1, to the experimental design by examining the derivative of the production values. In fact, a simulation addition criterion can be based on: (1) the value of the derivative of the production values obtained by the flow simulator, (2) direct identification of points whose production value is maximum or (3) direct identification of points whose production value is minimum.

[0124] A model is determined that approaches the values of the derivatives at the points selected by the experimental design in step 2. Then, a new simulation point is added in the place where the response of the derivative model is zero, provided that this point is sufficiently distant from the simulations already performed. These confirmation points allow testing of the predictivity of the second model, in this new investigated zone. If the prediction residues calculated at the new selected points exceed a value set by the user, these new points are used to carry out a new interpolation step.

[0125] Adding simulations to the current model, whether it is the consequence of a lack of a priori or a posteriori predictivity, allows increasing the quality and the quantity of information on the response function so as to obtain a more representative sampling.

[0126] Step 5: Construction and Adjustment of a Third Model

[0127] From the second model determined in step 2, the residues are determined at the new simulation points selected in step 4. The residues correspond to the difference between the response of the first model and the simulation value obtained by the reservoir flow simulator. The residues are then interpolated. Any n-dimensional interpolation method is suitable. For example, kriging or the spline method can be used.

[0128] The residue interpolation structure is divided up into two parts: the first model determined in step 2, and a “correcting” term allowing making up the difference between the prediction of the first model and the new simulation(s) selected in step 4. The new simulation allows interpolation of the responses and, thus, to account for the non-linearities detected at the surface.

[0129] An adjusted second model is determined by adding the results of the interpolation of the residues to the first model determined in step 2.

[0130] Iteration

[0131] It is furthermore possible, according to the invention, to improve the model iteratively by repeating steps 4 and 5.

[0132] In this case, during the new step 4, simulations points are added in relation to the model determined during the previous step 5. During the new step 5, a new model is constructed and adjusted starting from the simulation points selected in the new step 4 and by adjusting the first model determined in step 2.

[0133] Step 6: Seeking Inflection Points

[0134] If the a posteriori method has been used in step 4, the model determined in step 5 can be improved by adding simulation points by carrying out the following steps:

[0135] (a) determining an analytical model expressing the derivative of reservoir production as a function of time, by adjusting to the derivatives at the points associated with the production values selected in steps 2 and 4;

[0136] (b) checking that, at the point added in step 4, the response of the analytical model expressing the reservoir production derivative is zero;

if this response is greater than 0, determining the maximum of the third model determined in step 5 in the vicinity of the point added in step 4;

if this response is less than 0, determining the minimum of the third model determined in step 5 in the vicinity of the point added in step 4,

[0137] (c) determining the value of the minimum or of the maximum by the flow simulator; and

[0138] (d) determining a new model by adjusting the third model so that the response of the new model corresponds to the new minimum or maximum value obtained by the flow simulator.

[0139] The advantage of the method according to the invention is illustrated hereafter in connection with FIGS. 2 and 3.

[0140] The highly non-linear analytical function which was studied comprises two parameters x and y in order to better visualize the results. It is the “camel” function, which is characterized by its high non-linearity. The expression of this function is as follows:

[0141] $F\left(x,y\right)=4x^4-\frac{21}{10}{x}^4+\frac{1}{3}{x}^6+\mathrm{xy}-4y^2+4y^4$

[0142] It is graphically represented in the unit cube [−1,1]² bearing reference A in FIG. 2.

[0143] Reference B in FIG. 2 is the graph of the estimation of the “camel” function by a linear model obtained from a 4-simulation factorial design. Reference C in FIG. 2 is the graph of the estimation of the “camel” function by a polynomial of the second order obtained from a 9-simulation centered composite design.

[0144] The disparity of the results between, on the one hand, the function to be modelled (cube A) and, on the other hand, the models (cubes B and C) confirms the limits of the theory of conventional experimental designs for modelling non-linear functions.

[0145] FIG. 3 illustrates the optimization, according to our invention, of the model approaching the “camel” function. The function represented in the unit cube [−1,1]² bearing reference D is obtained by carrying out steps 2 and 3 from a Latin hypercube of initial maximum distance containing nine tests. Then, the functions represented in the unit cube [−1,1]² bearing references E, F and G are obtained by adjusting this function obtained from a Latin hypercube and by adding seven simulation points. Stages 4 and 5 are repeated three times.

[0146] By comparing function G in FIG. 3 with the “camel” function A of FIG. 2, the curves are seen to be relatively close to one another with the non-linearities being clearly detected. The method according to the invention is suitable and the results are very satisfactory.