A multivariate Bernstein copula model for permeability stochastic simulation

array:23 ["pii" => "S0016716914714989" "issn" => "00167169" "doi" => "10.1016/S0016-7169(14)71498-9" "estado" => "S300" "fechaPublicacion" => "2014-04-01" "aid" => "71498" "copyright" => "Universidad Nacional Autónoma de México" "copyrightAnyo" => "2014" "documento" => "article" "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/" "subdocumento" => "fla" "cita" => "Geofisica Internacional. 2014;53:163-81" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 1298 "formatos" => array:3 [ "EPUB" => 25 "HTML" => 843 "PDF" => 430 ] ] "itemSiguiente" => array:18 [ "pii" => "S0016716914714990" "issn" => "00167169" "doi" => "10.1016/S0016-7169(14)71499-0" "estado" => "S300" "fechaPublicacion" => "2014-04-01" "aid" => "71499" "copyright" => "Universidad Nacional Autónoma de México" "documento" => "article" "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/" "subdocumento" => "fla" "cita" => "Geofisica Internacional. 2014;53:183-98" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 1383 "formatos" => array:3 [ "EPUB" => 39 "HTML" => 936 "PDF" => 408 ] ] "es" => array:11 [ "idiomaDefecto" => true "titulo" => "Joint interpretation of geoelectrical and volatile organic compounds data: a case study in a hydrocarbons contaminated urban site" "tienePdf" => "es" "tieneTextoCompleto" => "es" "tieneResumen" => array:2 [ 0 => "es" 1 => "en" ] "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "183" "paginaFinal" => "198" ] ] "contieneResumen" => array:2 [ "es" => true "en" => true ] "contieneTextoCompleto" => array:1 [ "es" => true ] "contienePdf" => array:1 [ "es" => true ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:7 [ "identificador" => "f0045" "etiqueta" => "Figure 9" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr9.jpeg" "Alto" => 854 "Ancho" => 1221 "Tamanyo" => 74157 ] ] "descripcion" => array:1 [ "es" => "BTEX ratio of free-phase gasoline samples." ] ] ] "autores" => array:2 [ 0 => array:2 [ "autoresLista" => "Omar Delgado-Rodríguez, David Flores-Hernández, Myriam A. Amezcua-Allieri, Andrés Rosas-Molina, Salvador Marín-Córdova" "autores" => array:5 [ 0 => array:2 [ "nombre" => "Omar" "apellidos" => "Delgado-Rodríguez" ] 1 => array:2 [ "nombre" => "David" "apellidos" => "Flores-Hernández" ] 2 => array:2 [ "nombre" => "Myriam A." "apellidos" => "Amezcua-Allieri" ] 3 => array:2 [ "nombre" => "Andrés" "apellidos" => "Rosas-Molina" ] 4 => array:2 [ "nombre" => "Salvador" "apellidos" => "Marín-Córdova" ] ] ] 1 => array:2 [ "autoresLista" => "Vladimir Shevnin" "autores" => array:1 [ 0 => array:2 [ "nombre" => "Vladimir" "apellidos" => "Shevnin" ] ] ] ] ] "idiomaDefecto" => "es" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714990?idApp=UINPBA00004N" "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714990/v2_201505081358/es/main.assets" ] "itemAnterior" => array:18 [ "pii" => "S0016716914714977" "issn" => "00167169" "doi" => "10.1016/S0016-7169(14)71497-7" "estado" => "S300" "fechaPublicacion" => "2014-04-01" "aid" => "71497" "copyright" => "Universidad Nacional Autónoma de México" "documento" => "article" "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/" "subdocumento" => "fla" "cita" => "Geofisica Internacional. 2014;53:153-62" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 825 "formatos" => array:3 [ "EPUB" => 24 "HTML" => 507 "PDF" => 294 ] ] "en" => array:11 [ "idiomaDefecto" => true "titulo" => "Sensitivity of the surface temperature to changes in total solar irradiance calculated with the WRF model" "tienePdf" => "en" "tieneTextoCompleto" => "en" "tieneResumen" => array:2 [ 0 => "es" 1 => "en" ] "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "153" "paginaFinal" => "162" ] ] "contieneResumen" => array:2 [ "es" => true "en" => true ] "contieneTextoCompleto" => array:1 [ "en" => true ] "contienePdf" => array:1 [ "en" => true ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:7 [ "identificador" => "f0035" "etiqueta" => "Figure 7" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr7.jpeg" "Alto" => 1022 "Ancho" => 1879 "Tamanyo" => 195044 ] ] "descripcion" => array:1 [ "en" => "Daily temperature sensitivity in: a) The whole area of study. b) Zone A. c) Zone B. d) Zone C. e) Zone D. f) Zone E." ] ] ] "autores" => array:3 [ 0 => array:2 [ "autoresLista" => "Carolina Cipagauta" "autores" => array:1 [ 0 => array:2 [ "nombre" => "Carolina" "apellidos" => "Cipagauta" ] ] ] 1 => array:2 [ "autoresLista" => "Blanca Mendoza" "autores" => array:1 [ 0 => array:2 [ "nombre" => "Blanca" "apellidos" => "Mendoza" ] ] ] 2 => array:2 [ "autoresLista" => "Jorge Zavala-Hidalgo" "autores" => array:1 [ 0 => array:2 [ "nombre" => "Jorge" "apellidos" => "Zavala-Hidalgo" ] ] ] ] ] "idiomaDefecto" => "en" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714977?idApp=UINPBA00004N" "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714977/v2_201505081358/en/main.assets" ] "en" => array:18 [ "idiomaDefecto" => true "titulo" => "A multivariate Bernstein copula model for permeability stochastic simulation" "tieneTextoCompleto" => true "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "163" "paginaFinal" => "181" ] ] "autores" => array:2 [ 0 => array:4 [ "autoresLista" => "Victor Hernández-Maldonado, Martín Díaz-Viera" "autores" => array:2 [ 0 => array:4 [ "nombre" => "Victor" "apellidos" => "Hernández-Maldonado" "email" => array:1 [ 0 => "vmhernann@yahoo.com.mx" ] "referencia" => array:1 [ 0 => array:2 [ "etiqueta" => "*" "identificador" => "cor0005" ] ] ] 1 => array:2 [ "nombre" => "Martín" "apellidos" => "Díaz-Viera" ] ] "afiliaciones" => array:1 [ 0 => array:2 [ "entidad" => "Instituto Mexicano del Petróleo, Eje Central Lázaro Cárdenas 152, 07730, México, D.F." "identificador" => "aff0005" ] ] "correspondencia" => array:1 [ 0 => array:3 [ "identificador" => "cor0005" "etiqueta" => "*" "correspondencia" => "Corresponding author:" ] ] ] 1 => array:3 [ "autoresLista" => "Arturo Erdely" "autores" => array:1 [ 0 => array:2 [ "nombre" => "Arturo" "apellidos" => "Erdely" ] ] "afiliaciones" => array:1 [ 0 => array:2 [ "entidad" => "Programa de Actuaria, Facultad de Estudios Superiores Acatlán, Universidad Nacional Autónoma de México, México D.F." "identificador" => "aff0010" ] ] ] ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:7 [ "identificador" => "f0040" "etiqueta" => "Figure 8" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr8.jpeg" "Alto" => 1081 "Ancho" => 1087 "Tamanyo" => 130538 ] ] "descripcion" => array:1 [ "en" => "Scatter plot and histograms of VS Meas and K (Simulated values)." ] ] ] "textoCompleto" => "<a name="p164"></a>IntroductionIntegrated Reservoir Modeling (IRM) is the most accepted way to obtain the spatial distribution of petrophysical properties in oilfields (<a class="elsevierStyleCrossRef" href="#bib0030">Cosentino, 2001</a>). An important and common task performed in this method is the estimation of permeability, because it is well known that this petrophysical property is quite informative about the oil flux patterns in a reservoir. However, it is difficult to obtain direct information about permeability, and therefore it is necessary to find models of dependence with another petrophysical property (such as porosity, water saturation, etc.) in order to have an estimation of its profile (<a class="elsevierStyleCrossRef" href="#bib0145">Landa et al., 1996</a>).The linear regression approach is the most common way to model permeability values using other petrophysical properties as covariates (<a class="elsevierStyleCrossRef" href="#bib0015">Balan et al., 1995</a>). To meet the requirement of linearity it is common to perform transformations which imply that the final result could be biased when it is back transformed; there are approaches that, in order to fit linear models, apply logarithmic transformations to induce this behavior, for example, to use the Cokriging method it is necessary to have a linear relationship since it requires the linear corregionalization model (Sanjay and Journel 1994).The main disadvantage of linear dependency models is their lack of ability to capture and model the dependence structure or pattern (<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy et al. 2005</a>). In other words, traditional methods cannot capture the complex variability of data, in terms of variance or standard deviation; hence, the predicted permeability will not reproduce extreme values of the real data. In other words, these approaches will not be able to represent impermeable barriers or high permeability zones, and from a fluid flow point of view this aspects are the most important characteristics that determine the patterns of fluid movement. In this context, the predicted permeability profile, using linear estimators, will not be an effective approximation due to its oversmoothing nature.On the other hand, model-free function estimators like artificial neural networks are very flexible tools that have been used to model permeability. However, neural networks have some disadvantages too. First, the training process has to be done with caution and can be a lengthy process. The good results obtained by this technique are reached using a comprehensive training data set, which is not always available. On the other hand, failing in correctly calibrating the network may result in aberrant results. Another point to take into consideration is that the methodology is not yet an “off the shelf” application and requires expertise by the geoscientist (<a class="elsevierStyleCrossRef" href="#bib0030">Cosentino, 2001</a>).Another alternative are Bayesian methods; however, the traditional framework of the Bayesian analysis is based on the multivariate normal distribution where the lower and upper tails are symmetrical. <a class="elsevierStyleCrossRef" href="#bib0010">Armstrong et al. (2004)</a> proposed an alternative Bayesian analysis that is based on Archimedean copulas where the joint distribution does not have to be normal and there is flexibility to have a lower tail or upper tail dependence based on the specific type of copula.Constructing numerical models of the reservoir that honor all available data (core measurements, well logs, seismic and geological interpretations, etc.) having sparse knowledge of rock properties, leads us to consider the stochastic simulation approach (<a class="elsevierStyleCrossRef" href="#bib0060">Deutsch, 1992</a>). This is not a new concept (<a class="elsevierStyleCrossRef" href="#bib0105">Haldorsen and Damsleth, 1990</a>; <a class="elsevierStyleCrossRef" href="#bib0125">Journel and Alabert, 1990</a>), stochastic models of physical systems are used extensively in many disciplines.Stochastic simulation is the process of building alternate, equally probable models of the spatial distribution of a random function. It is said that a simulation is conditional if the resulting realizations honor the raw data values at their locations. The most straightforward algorithms for generating realizations of a multivariate Gaussian field is provided by Sequential Gaussian Simulation (SGS) and Sequential Indicator Simulation (SIS), which are extensively used to perform permeability simulations (<a class="elsevierStyleCrossRef" href="#bib0120">Holden et al. 1995</a>). However despite of their improvements (<a class="elsevierStyleCrossRef" href="#bib0130">Journel and Zhu, 1990</a>; Suro-Perez, 1990) these methods are limited to cases when the spatial continuity is characterized by stationary two-point statistics and to data that is defined on the same support (<a class="elsevierStyleCrossRef" href="#bib0060">Deutsch, 1992</a>).A competitive and more systematic method for predicting permeability may be achieved by applying stochastic joint simulations, in which the correct specification of dependence pattern between petrophysical properties is crucial (<a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham, 1994</a>). According to Deutsch this approach basically consists of an annealing geostatistical cosimulation of porosity-permeability using their empirical joint distribution.<a name="p165"></a>A modification of Deutsch’s methodology was presented by <a class="elsevierStyleCrossRef" href="#bib0070">Díaz-Viera and Casar-González, (2005)</a>. Here, it was proposed the use of a bivariate t-copula to construct the joint distribution function between porosity and permeability rather than use the joint distribution function of the sample data, the dependence structure was specified by matching concordance measures such as Kendall’s τ and Spearman’s ρ. The above methodology was applied to simulate permeability from the porosity profile in a double porosity carbonate systems, restricted to one-dimensional case at well-log scale (<a class="elsevierStyleCrossRef" href="#bib0075">Díaz-Viera et al. 2006</a>).While Díaz-Viera and Casar-González proposal can reproduce adequately the observed data, and also their extreme values using a t-copula, the critical problem is that the applied copula is parametric, and consequently, it is based on a given distribution function, the student t. In that sense, expecting a single copula family to be able to model any kind of dependency relationship seems to be too restrictive, at least for the petrophysical properties under consideration. It does not mean that parametric copulas have not practical use. It means that petrophysical properties modeling, where dependence structure can be complex or non linear, using parametric copulas could be available but in a very complicated form, <a class="elsevierStyleCrossRef" href="#bib0170">Sancetta and Satchell (2004)</a>.The copula approach has been successfully used to model dependence patterns in few areas of oil industry, for example, for field development decision process (<a class="elsevierStyleCrossRef" href="#bib0050">Accioly and Chiyoshi 2004</a>) or to model dependence in petroleum decision making (<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy et al. 2005</a>).In a recent work we use the copula approach to model the bivariate dependence between petrophysical properties in a complete nonparametric fashion (<a class="elsevierStyleCrossRef" href="#bib0110">Hernández-Maldonado et al. 2012</a>), it means we did not impose a priori a parametric joint distribution function to characterize the dependence structure of the sample; instead, we used a bivariate Bernstein copula function to model the intrinsic distribution of data values. In this work we introduce a multivariate method for estimating empirical dependence among several petrophysical properties using a nonparametric copula. The dependency model obtained is then used to stochastically simulate one property (primary variable) given other ones (secondary variables). Here, we apply a method widely used in geostatistics, the Simulated Annealing method, which is a global optimization framework where we can add restrictions to simulation. In this case, in addition of using multivariate dependence by sampling the empirical copula model, a spatial correlation function (a variogram) for the secondary variable is imposed. Additionally, if any raw data of the primary variable is available, it can be exactly honored, that is, the method can be conditional.In this paper this method is applied to simulate rock permeability using a trivariate copula model where permeability is described by Vugular Porosity and Shear Wave Velocity. It is introduced a two stage simulation method in a multivariate fashion to model stochastically the spatial distribution of permeability at well log scale. Conditional and nonconditional simulations, and the median of simulations are performed in order to show the results that this method provides; also it is established that, as far as more descriptive variables integrate the Bernstein copula it will reproduce permeability values in a very precise form, consequently it will not be necessary to perform many simulations to reduce small scale variability.Brief introduction to multivariate copulasThe main disadvantage of dependency models based on a linear regression approach is their lack of ability to capture and model more general dependence structure. As an alternative a copula can model the joint distribution of petrophysical properties in better ways. The essence of the copula approach is that a joint distribution of random variables can be expressed as a function of the marginal distributions (<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy et al., 2005</a>).According to Sklar’s theorem (<a class="elsevierStyleCrossRef" href="#bib0190">Sklar, 1959</a>), the underlying copula associated to a multivariate random vector (X1, X2, ..., Xm) represents a functional link between the joint probability distribution and the univariate marginal distributions F1, F2, ..., FM respectively:<elsevierMultimedia ident="eq0005"></elsevierMultimedia>For all x1, x2, ..., xm in the extended real numbers system, where C:[0, 1]m →[0, 1] the underlying copula is unique whenever X1, X2, and Xm are continuous random variables. Therefore, all the information about the dependence between continuous random variables is contained in their corresponding copula. Several properties may be derived for copulas (<a class="elsevierStyleCrossRef" href="#bib0155">Nelsen, 2006</a>), and among them <a name="p166"></a>we have an immediate corollary from Sklar’s theorem: X1, X2 and Xm are independent continuous random variables if and only if their underlying copula is C(u1,...,un)=u,...,un.Let S={(x11, x21, ... xm1), ..., (x1n, x2n, ...xmn)} be n observations of a random vector (X1, X2, ..., Xm). We may obtain empirical estimates for the marginal distributions X1, X2, ..., Xm by means of:<elsevierMultimedia ident="eq0010"></elsevierMultimedia>where stands for an indicator function which takes value whenever its argument is true, and otherwise. It is well-known (<a class="elsevierStyleCrossRef" href="#bib0025">Billingsley, 1995</a>) that the empirical distribution Fˆj is a consistent estimator of Fj that is, Fj that is Fˆj(t)is a converges almost surely to Fj(t) as n →∞ all t.Similarly, we have the empirical copula (<a class="elsevierStyleCrossRef" href="#bib0055">Deheuvels, 1979</a>), a function Cn with domain {1n;i=0,1,…n}n (where n is the sample size) defined as<elsevierMultimedia ident="eq0015"></elsevierMultimedia>and its convergence to the true copula has also been proved (<a class="elsevierStyleCrossRef" href="#bib0100">Fermanian et al., 2004</a>). The empirical copula is not a copula, since it is only defined on a finite grid, not in the whole unit hypercube [0, 1]n, but by Sklar’s Theorem (<a class="elsevierStyleCrossRef" href="#bib0190">Sklar, 1959</a>) may be extended to a copula.Cn may be extended to a copula.Sklar’s theorem is completely general and a joint distribution function can be constructed using the copula function. The copula separates the marginal distributions from correlation and the copula itself can capture the dependence structure. This is an essential property of copulas.A general example of a multivariate copula model that describes the relationship of three random variables (a trivariate copula) according to Sklar’s theorem is shown:<elsevierMultimedia ident="eq0020"></elsevierMultimedia>There are different classes and families of copulas but, for this study, we will work with the Bernstein copula, which is a nonparametric copula since it is a nonparametric smoothing based on the empirical copula.Joint stochastic simulation method using a multivariate Bernstein copulaThe method presented here basically consists in generating a stochastic simulation of a primary variable, where the joint distribution function of its covariates is modeled using a nonparametric copula. In other words, the spatial dependence and dependence pattern between variables are decoupled; in this context, the Bernstein copula only models the dependence between variables and the spatial dependence is modeled by a variogram applying the Simulated Annealing method.Basically this proposal is a two-stage modeling algorithm method; first, the pattern relationship between petrophysical properties under study; and second, the spatial structure. In the first stage it is modeled the dependence structure of petrophysical properties using a multivariate nonparametric copula (the Bernstein copula), then a geostatistical simulation of primary variable is performed using simulated annealing technique, whose objective function is the variogram model (<a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham, 1994</a>). A detailed description of each step of the algorithm will be described.Multivariate Bernstein Copula ModelingEach petrophysical property is modeled as an absolutely continuous random variable X with unknown marginal distribution function F. For simulation of continuous random variables, the use of the empirical distribution function <a class="elsevierStyleCrossRef" href="#eq0010">(2)</a> is not appropriate since Fˆj is a step function, and therefore discontinuous, so a smoothing technique is needed. Since our main goal is to simulate a primary variable using more than one descriptive variable, it will be better to have a smooth estimation of the marginal quantile function Q(u)=F-1(u)=inf{x:F(x)≥u},0≤u≤1 which is possible by means of Bernstein polynomials as in <a class="elsevierStyleCrossRef" href="#bib0150">Munoz-Perez and Fernandez-Palacín (1987)</a>.<elsevierMultimedia ident="eq0025"></elsevierMultimedia>For a smooth estimation of the underlying copula we make use of the Bernstein copula (<a class="elsevierStyleCrossRef" href="#bib0170">Sancetta and Satchell, 2004</a>; <a class="elsevierStyleCrossRef" href="#bib0175">Sancetta, 2007</a>):<elsevierMultimedia ident="eq0030"></elsevierMultimedia><a name="p167"></a>where;<elsevierMultimedia ident="eq0035"></elsevierMultimedia>For every (u1, …, um) in the unit hypercube [0, 1]m, and αv1n1,…,vmnm is the empirical copula, defined in <a class="elsevierStyleCrossRef" href="#eq0015">(3)</a>.For example, the trivariate Bernstein copula model derived from <a class="elsevierStyleCrossRef" href="#eq0030">(6)</a> is:<elsevierMultimedia ident="eq0040"></elsevierMultimedia>Sampling algorithm of a Bernstein copula modelFor a pair of random variables (X1, X2) with joint distribution function H and underlying copula C we need to generate an observation of uniform (0, 1) random variables (U, V) whose joint distribution function is C and then transform those uniform variables as in step 3 of the sampling bivariate algorithm. For generating such pair (u, v) it is used a conditional distribution method, this method needs the conditional distribution function for V given U=u, which we denote as cu (v):<elsevierMultimedia ident="eq0045"></elsevierMultimedia>where CˆB is the bivariate Bernstein copula model, obtained by <a class="elsevierStyleCrossRef" href="#eq0030">(6)</a>.To simulate replications from the random vector with the dependence structure estimated from the observed data S:={(x11, x21), ..., (x1n, x2n)} it is applied the following algorithm:Sampling bivariate algorithm<ul class="elsevierStyleList" id="lis0005"><li class="elsevierStyleListItem" id="lsti0005">1Generate two independent and continuous Uniform (0,1) random variatesu andt</li><li class="elsevierStyleListItem" id="lsti0010">2Setv=cu-1 (t); wherecu is defined in <a class="elsevierStyleCrossRef" href="#eq0045">(9)</a>.</li><li class="elsevierStyleListItem" id="lsti0015">3The desired pair is (x1,x2)=(Q˜;n(u),R˜;n(v)) where Q˜;n and R˜;n according to <a class="elsevierStyleCrossRef" href="#eq0025">(5)</a>, are the estimated and smoothed quantile functionsX1 ofX2 and, respectively.</li></ul>For the multivariate case we must solve equations that represent conditional distribution functions for W given U=u,V=v.To simulate replications from the random vector with the dependence structure estimated from the observed data S:={(x11, x21, x31), ..., (x1n, x2n, x3n)} it is applied the next algorithm.Sampling trivariate algorithm<ul class="elsevierStyleList" id="lis0010"><li class="elsevierStyleListItem" id="lsti0020">1Generate three independent and continuous Uniform (0, 1) random variablesu andt1, t2</li><li class="elsevierStyleListItem" id="lsti0025">2Setv=cu-1 (t1) wherecu is defined in <a class="elsevierStyleCrossRef" href="#eq0045">(9)</a>.</li><li class="elsevierStyleListItem" id="lsti0030">3Setw=cuv-1 (t2) where <elsevierMultimedia ident="eq0050"></elsevierMultimedia></li></ul>where is the trivariate Bernstein copula model <a class="elsevierStyleCrossRef" href="#eq0040">(8)</a>. <ul class="elsevierStyleList" id="l1035"><li class="elsevierStyleListItem" id="lsti0035">4The desired vector is (x1,x2,x3)=(Q˜;n(u),R˜;n(v),H˜;n(w)) where Q˜;n(u),R˜;n(v) and, H˜;n(w) according to <a class="elsevierStyleCrossRef" href="#eq0025">(5)</a>, are the estimated and smoothed quantile functions ofX1,X2, yX3, respectively.</li></ul>Simulated Annealing methodTo perform stochastic simulations of a primary variable it is applied a method widely used in geostatistics, the Simulated Annealing method, which is a global optimization framework where we can add restrictions to simulation. In this case, in addition of using multivariate dependence by sampling the copula model, a spatial correlation function (a variogram) for the secondary variable is imposed. Additionally, if any raw data of the primary variable is available, it can be exactly honored, that is, the method can be conditional.A more detailed explanation of each single step of simulated annealing method, as well as technical terminology, can be found in <a class="elsevierStyleCrossRef" href="#bib0080">Dreo et al. (2006)</a>, Dafflon et. al., (<a class="elsevierStyleCrossRef" href="#bib0035">2009a</a>, <a class="elsevierStyleCrossRef" href="#bib0040">2009b</a>), <a class="elsevierStyleCrossRef" href="#bib0045">Dafflon and Barrash (2012)</a>, <a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado et al. (2012)</a>, Deutsch (2002). In this work, we will just briefly mention its steps.a. Generating an initial configuration, known as “seed”. An initial configuration is the starting point of the simulation, also, it can be <a name="p168"></a>considered as a “possible solution” of the spatial disposition of the primary variable. There are some ways of proposing this configuration, see Deutsch (2002), <a class="elsevierStyleCrossRef" href="#bib0080">Dréo et al. (2006)</a>. One of most popular of them is just proposing random numbers with uniform distribution. We decided to use the multivariate Bernstein copula sampling algorithm (section 3.1.1) to obtain this starting point in order to have a better approximation of the final solution.b. Defining the objective function. Since we use the Bernstein copula we do not have to include univariate histograms into the objective function because they are reproduced automatically; in the same way, as Bernstein copula models the dependence structure between variables, it is not necessary to include neither a correlation coefficient. In other words, in order to satisfy the spatial distribution of the primary variable, we are able to propose an objective function that consists only in one term <a class="elsevierStyleCrossRef" href="#eq0055">(11)</a> the variogram function<elsevierMultimedia ident="eq0055"></elsevierMultimedia>Where:<elsevierMultimedia ident="eq0060"></elsevierMultimedia><elsevierMultimedia ident="eq0065"></elsevierMultimedia>c. Measuring the initial configuration energy. By <a class="elsevierStyleCrossRef" href="#eq0060">(12)</a>, it is calculated the variogram of the initial configuration and it is compared to the variogram of the real configuration <a class="elsevierStyleCrossRef" href="#eq0055">(11)</a>, this value is the energy of the starting point, and the main idea is that it must be decreased during the simulation procedure, see Deutsch (2002). This step gives us an idea of the seed´s quality, if we use random numbers to generate it, its measured energy will be very high, however, as long as this method uses a sample given by the Bernstein copula the initial energy is significantly small.d. Obtaining the initial temperature and the annealing schedule. The initial temperature is obtained by <a class="elsevierStyleCrossRef" href="#eq0065">(13)</a>:The annealing schedule is obtained following the procedure proposed by <a class="elsevierStyleCrossRef" href="#bib0080">Dréo et al. (2006)</a>.e. Ending up a realization. The simulation ends whenever the objective function error is reached, or an accumulation of 3 stages without change occurs, or when the maximum number of attempted perturbations is reached.It is shown a diagram of the method to perform multivariate stochastic simulations in, <a class="elsevierStyleCrossRef" href="#f0005">Figure1</a>.<elsevierMultimedia ident="f0005"></elsevierMultimedia>Case studyIt will be modeled the Permeability of double porosity carbonate formations of a South Florida Aquifer in the western Hillsboro Basin of Palm Beach County, Florida. Based in the algorithm described in Section 3, we propose a copula-based approach to model the relationship between permeability, porosity and Shear Wave Velocity (S-waves). The dependency model obtained will be then used to stochastically simulate permeability using Porosity and Shear Wave velocity.Data descriptionThe characterization of this aquifer for the borehole and field scales is given in <a class="elsevierStyleCrossRef" href="#bib0160">Parra et al. (2001)</a>, <a class="elsevierStyleCrossRef" href="#bib0165">Parra and Hackert (2002)</a>, and a hydrogeological situation is described by <a class="elsevierStyleCrossRef" href="#bib0020">Bennett et al. (2002)</a>. The interpretation <a name="p169"></a>of the borehole data, the determination of the matrix and secondary porosity and secondary-pore types (shapes of spheroids approximating secondary pores) were presented by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>, where to determine the pore microstructure of aquifer carbonate formations the authors applied the petrophysical inversion technique that consists in minimizing a cost function that includes the sum of weighted square differences between the experimentally measured and theoretically calculated logs as in <a class="elsevierStyleCrossRef" href="#bib0135">Kazatchenko et al. (2004)</a>.In this case the following well logs were used for joint simultaneous inversion as input data: resistivity log, transit times of the Pand S-waves (acoustic log), total porosity (neutron log), and formation density (density log). To calculate the theoretical acoustic and resistivity logs the double-porosity model for describing carbonate formations was applied: <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>.This model treats carbonate rocks as a composite material that consists of a homogeneous isotropic matrix (solid skeleton and matrix pore system) where the secondary pores of different shapes are embedded. The secondary pores were approximated by spheroids with variable aspect ratios to represent different secondary porosity types: vugs (close-to-sphere shapes), quasi-vugs (oblate spheroids), channels (prolate spheroids), and microfractures (flattened spheroids) <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>.We used the results of inversion obtained by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a> for carbonate formations of South Florida Aquifer that includes the following petrophysical characteristics: matrix porosity, secondary vugular and crack porosities (<a class="elsevierStyleCrossRef" href="#f0010">Figure 2</a>).<elsevierMultimedia ident="f0010"></elsevierMultimedia>The secondary-porosity system of this formation has complex microstructure and corresponds to a model with two types of pore shapes: cracks (flattened ellipsoids) with the overall porosity of 2% and vugs (close to sphere) with the porosity variations in the range of 10-30%.Statistical data analysisGiven multivariate data, it is common to start choosing as explanatory variables those who exhibit higher dependence with the variable that is to be explained. We measured the dependence between the petrophysical properties in terms of the dependence index Φ proposed by <a class="elsevierStyleCrossRef" href="#bib0115">Hoeffding (1941)</a>, which satisfies all desirable properties for a dependence measure for continuous random variables, see <a class="elsevierStyleCrossRef" href="#bib0155">Nelsen (2006)</a>.From all possible explanatory variables (PHITOT, PHIMAT, PHICR, PHIV, VS Meas, VP meas) for permeability (K), for the first explanatory random variable it was chosen relative vugular porosity (PHIV) since it exhibited the highest dependence Φ (PHI, K)=0.71 (on a [0,1] scale). In choosing a second explanatory random variable we need, in addition, to have a high dependence with permeability, and the lowest possible dependence with the first explanatory variable (PHIV), otherwise it would mean that it is quite similar to it and it will add no significant information to what the first one already can provide. Under this criteria, the second best choice was Share Wave Velocity (VS meas), <a name="p170"></a>with Φ (VS meas, K)=0.60 and Φ (PHIV, VS meas)=0.55]<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a> shows all dependence index Φof explanatory variables that are the results of inversion obtained by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>; Total porosity (PHITOT), Matrix porosity (PHIMAT), Crack porosity (PHICR), Vugular porosity (PHIV), Share Wave Velocity (VS Meas), P-Wave velocity (VP meas).<elsevierMultimedia ident="tbl0005"></elsevierMultimedia>In <a class="elsevierStyleCrossRef" href="#tbl0005">table 1</a> explanatory variables are shown, using this information we confirm that PHIV has the highest dependence, but following the order of this table we should use as second explanatory variable PHITOT. We did not make this choice because, as we said before, a second explanatory variable needs to have a high dependence with permeability, and the lowest possible dependence with the first explanatory variable (PHIV). To see this behavior in table two we show the explanatory variables matrix dependence Hoeffding index.We have chosen PHIV as the first explanatory variable; it means we must pay attention to the first line of <a class="elsevierStyleCrossRef" href="#tbl0010">table 2</a>, and select the variable with the lowest dependence to PHIV. It should be PHICR but in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a> we see that this property offers poor information about K. Finally, we may choose VP meas but, VS Meas explains K in some better way (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>).<elsevierMultimedia ident="tbl0010"></elsevierMultimedia>In <a class="elsevierStyleCrossRef" href="#f0015">Figure 3</a> it is shown a scatter-plot distribution and histograms of PHIV–K sample Explanatory data taken from <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>, it shows their bivariate distribution. In <a class="elsevierStyleCrossRef" href="#f0020">Figure 4</a> it is shown a scatter-plot distribution and histograms of VS meas–K data values<elsevierMultimedia ident="f0015"></elsevierMultimedia><elsevierMultimedia ident="f0020"></elsevierMultimedia>Bernstein Copula ModelIt is modeled the trivariate Bernstein Copula using relative vugular porosities (PHIV), Shear Wave Velocity (VS meas=velocity of S-waves measured) and K as absolutely continuous random variables X1, X2 and X3, with unknown theoretical univariate distribution functions F1, F2and F3, see<a class="elsevierStyleCrossRef" href="#eq0010">(2)</a>. Observations of this random vector were obtained from<a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>.We use <a class="elsevierStyleCrossRef" href="#eq0010">(2)</a> to estimate the empirical distribution function for each variable; however; <a name="p171"></a>in order to model them as continuous random variables they should be smoothed using <a class="elsevierStyleCrossRef" href="#eq0025">(5)</a>, since <a class="elsevierStyleCrossRef" href="#eq0010">(2)</a> is not appropriate because Fˆj is a step function and therefore discontinuous. The trivariate Bernstein copula model <a class="elsevierStyleCrossRef" href="#eq0040">(8)</a> will be used to sample data to stochastically simulate permeability.Sampling from Bernstein copulaIn <a class="elsevierStyleCrossRef" href="#f0025">Figure 5</a> (left side) it is shown a scatterplot distribution and histograms of porosity– permeability real data (sample size of) taken from <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko et al. (2006)</a>. At the center of the same figure, it is shown a trivariate simulation of permeability values using vugular porosity and shear wave velocity as explicative variables. If we do not label these two scatter plots it would be difficult to know which one of them is the original one. In this case the multivariate Bernstein copula provides an acceptable the dependence structure for continuous variables with complex relationship. The right side of <a class="elsevierStyleCrossRef" href="#f0025">Figure 5</a> shows 3800 simulated values of permeability, according to some analysis, it shows more clearly percolation patterns.<elsevierMultimedia ident="f0025"></elsevierMultimedia>Permeability simulations and discussionIn this section spatial simulation will be performed using simulated annealing technique. Since multivariate dependence is modeled by the Bernstein copula, a spatial correlation function (the variogram) will be used as an objective function in the annealing algorithm.We propose a copula-based nonparametric approach and simulated annealing technique to model the spatial distribution of permeability using as secondary variables, porosity and shear wave velocity (S-waves). In <a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a> we can see a single nonconditional spatial simulation of permeability (K) in terms of depth. As we note the simulation (blue line) follows very close the real data behavior (red line), this is because we introduce more descriptive information into the Bernstein copula model. However, there is still a small-scale variability which does not exist in the original data values.<elsevierMultimedia ident="f0030"></elsevierMultimedia>As a descriptive measure of the goodness of fit of predicted values of log-permeability (Log K) given values of the explanatory variables PHIV and VS meas, we calculated the Mean Squared Error (MSE) in each case:<elsevierMultimedia ident="eq0070"></elsevierMultimedia><elsevierMultimedia ident="eq0075"></elsevierMultimedia><elsevierMultimedia ident="eq0080"></elsevierMultimedia>where a lower value of MSE represents a better goodness of fit.The simulated values of log K, given PHIV and VS meas (log K | PHIV, VS meas), show a MSE value of 0.5, which is very low compared to that obtained by <a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham’s method (1994)</a>, which MSE value is 7.76; also it is lower than the obtained one by Díaz-Viera and Casar-González method, where the MSE value is 5.5, details are shown in <a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado et al., (2012)</a>. Finally, the MSE value for the trivariate copula simulation is lower than the MSE value of the bivariate copula simulation (log K | PHIV)<a class="elsevierStyleCrossRef" href="#eq0070">(14)</a>, consequently it means that trivariate simulation has the best fitted values.<a name="p172"></a>In terms of histograms, scatterplots and variograms, the proposed method shows very good agreement with the original data values. As we are using three variables, the scatterplots of permeability can be seen in terms of K-PHIV (<a class="elsevierStyleCrossRef" href="#f0035">Figure 7</a>) and K-S-waves (<a class="elsevierStyleCrossRef" href="#f0040">Figure 8</a>). The scatterplots of the original data values are shown in <a class="elsevierStyleCrossRef" href="#f0015">Figure 3</a> and <a class="elsevierStyleCrossRef" href="#f0020">Figure 4</a> respectively.<elsevierMultimedia ident="f0035"></elsevierMultimedia><elsevierMultimedia ident="f0040"></elsevierMultimedia>The variogram of the simulation (<a class="elsevierStyleCrossRef" href="#f0045">Figure 9</a>) and its fitted model (<a class="elsevierStyleCrossRef" href="#tbl0015">Table 3</a>), also shows good agreement with respect to real data (<a class="elsevierStyleCrossRef" href="#f0050">Figure 10</a>). We can note that both variograms are almost equal, but in a detailed view, the empirical variogram of the simulated values is a little softer than the original one.<elsevierMultimedia ident="f0045"></elsevierMultimedia><elsevierMultimedia ident="tbl0015"></elsevierMultimedia><elsevierMultimedia ident="f0050"></elsevierMultimedia>As we could anticipate, the dependence structure between these three petrophysical properties is well represented. The extreme values, and in general, all of most important statistics are quite well reproduced, (<a class="elsevierStyleCrossRef" href="#f0055">Figure 11</a>). In order to reduce small-scale variability seen in <a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a>, a median of 10 nonconditional SA simulations was performed, <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a>. Once again, the spatial structure, the histograms and the Scatter-plot reproduction shows a very good agreement respecting to real data, see also <a class="elsevierStyleCrossRef" href="#f0065">Figure 13</a> and <a class="elsevierStyleCrossRef" href="#f0070">Figure 14</a>.<elsevierMultimedia ident="f0055"></elsevierMultimedia><elsevierMultimedia ident="f0060"></elsevierMultimedia><elsevierMultimedia ident="f0065"></elsevierMultimedia><elsevierMultimedia ident="f0070"></elsevierMultimedia>In <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a> we can see that for small values of K the simulation is not following real data values too closely, there are two <a name="p173"></a><a name="p174"></a>reasons. First, this “behavior” is observed in all simulation methods that were performed in order to compare the simulated values for the bivariate case of this proposal [<a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado et al., (2012)</a>]. In that work the most accurate results were shown by Bernstein copula. In this particular case is hard to have a simulation that follow too closely the behavior of the real data values, because there is very little information. And second, this is a median value of 10 different simulations, which means that each simulation has different degrees of accuracy, which is normal.<a name="p175"></a><a name="p176"></a><a class="elsevierStyleCrossRef" href="#f0075">Figure 15</a> shows the variogram of this simulation and its fitted model. In contrast to <a class="elsevierStyleCrossRef" href="#f0050">Figure 10</a> (real data values) we can see that this method represents very well, the spatial structure of permeability.<elsevierMultimedia ident="f0075"></elsevierMultimedia>The values of the variogram model for a single simulation and the median of 10 simulations are very close to the real data as well (<a class="elsevierStyleCrossRef" href="#tbl0015">Table 3</a>). The same happens with their most important statistics (<a class="elsevierStyleCrossRef" href="#f0080">Figure 16</a>).<elsevierMultimedia ident="f0080"></elsevierMultimedia>Using <a class="elsevierStyleCrossRef" href="#eq0080">(16)</a> we can see that the MSE values between the real data and these nonconditional simulations are very similar (<a class="elsevierStyleCrossRef" href="#tbl0020">Table 4</a>). It is an indicative that it is not necessary to perform <a name="p177"></a>many simulations to smooth the small-scale variability, and may be unnecessary.<elsevierMultimedia ident="tbl0020"></elsevierMultimedia>Finally we perform conditional simulations and compare the MSE values in order to study the effect of different percentage of conditioning values. <a class="elsevierStyleCrossRef" href="#f0085">Figure 17</a> shows the results of the spatial distribution of the permeability in various percentages of conditioning values, from top to bottom order, it is conditioned to 10 percent, 50 percent and 90 percent of the conditioning data values. Note that results are, graphically, very similar to real data and their MSE values as well (see <a class="elsevierStyleCrossRef" href="#tbl0025">Table 5</a>). This shows that the proposed method in a multivariate fashion does not require big amounts of conditioning data to have consistent results.<elsevierMultimedia ident="f0085"></elsevierMultimedia><elsevierMultimedia ident="tbl0025"></elsevierMultimedia>The Bernstein copula <a class="elsevierStyleCrossRef" href="#eq0040">(8)</a> requires that the empirical copula be ready to be used all the time, hence, it is necessary to propose an efficient and fast algorithm in terms of data storage and computing speed to implement a multidimensional copula. We propose a procedure to efficiently implement a multidimensional empirical copula. Shortly described, this procedure fundamentally solves the two main disadvantages of using nonparametric copulas: it avoids redundant calculation, which increases its speed and just needs to storage the most important values of the generated hypercube.Of course the method can be extended to more dimensions, but the purpose of the article is to show its application to a well log scale, where we have more control and better data quality. Its application to more dimensions will be the subject of future work.<a name="p178"></a><a name="p179"></a>ConclusionsFrom a methodological point of view, a multivariate copula approach provides a very flexible tool to model complex dependence relationships of petrophysical properties, such as porosity, permeability and shear wave velocity, without imposing strong assumptions of linearity or log-linearity, and/or normality when we are modeling them as random variables.All the information about the dependence structure is contained in the underlying copula, and its estimation is being used, instead of the extreme information reduction that is done by the use of numerical measures such as the linear correlation coefficient, which under the presence of nonlinear dependence may become useless and/or quite misleading, see Embrechts et al. (<a class="elsevierStyleCrossRef" href="#bib0090">1999</a>, <a class="elsevierStyleCrossRef" href="#bib0095">2003</a>). For this reason it is very important to extend the method to higher dimensions, because if we have more than one relevant explanatory variable, we may use this tool to take them into account.The main advantage of the methodology used in this work is that it represents a straightforward way to perform nonparametric simulations, conditional and nonconditional. As we have said in the results of the application of this method, including more than one explanatory petrophysical property to leads to more accurate results. Consequently, it will not be necessary to perform many simulations (<a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a> vs. <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a>) or conditioning with high quantities of data values (<a class="elsevierStyleCrossRef" href="#f0085">Figure 17</a>) to obtain acceptable results (<a class="elsevierStyleCrossRef" href="#tbl0025">Table 5</a>).This method has an important improvement with respect to the bivariate case, not only because it provides better results, but also it reduces the number of simulations to decrease the small scale variability. In other words it is not necessary to perform 10 or more simulations. Even more, it is not necessary to have a great number of conditioning values; the conditioned simulations did not show a lot of variation in terms of MSE values.For future work, we think we can use a Bernstein copula to explore complex relationships between petrophysical properties and then fit, if possible, a semiparametric gluing copula, <a class="elsevierStyleCrossRef" href="#bib0180">Siburg and Stoimenov (2008)</a>. Also, we consider we may use another optimization method which could give improved results, for example, we can implement genetic algorithms or Estimation of Distribution Algorithms (EDAs). It is very important to remark that the use of multivariate copulas leads to significant serious computational challenges, for this reason it is necessary to propose better ways to calculate them, for example parallel computing." "textoCompletoSecciones" => array:1 [ "secciones" => array:13 [ 0 => array:3 [ "identificador" => "xres498003" "titulo" => "Resumen" "secciones" => array:1 [ 0 => array:1 [ "identificador" => "abst0005" ] ] ] 1 => array:2 [ "identificador" => "xpalclavsec519550" "titulo" => "Palabras clave" ] 2 => array:3 [ "identificador" => "xres498004" "titulo" => "Abstract" "secciones" => array:1 [ 0 => array:1 [ "identificador" => "abst0010" ] ] ] 3 => array:2 [ "identificador" => "xpalclavsec519551" "titulo" => "Key words" ] 4 => array:2 [ "identificador" => "sec0005" "titulo" => "Introduction" ] 5 => array:2 [ "identificador" => "sec0010" "titulo" => "Brief introduction to multivariate copulas" ] 6 => array:3 [ "identificador" => "sec0015" "titulo" => "Joint stochastic simulation method using a multivariate Bernstein copula" "secciones" => array:2 [ 0 => array:2 [ "identificador" => "sec0020" "titulo" => "Multivariate Bernstein Copula Modeling" ] 1 => array:2 [ "identificador" => "sec0025" "titulo" => "Sampling algorithm of a Bernstein copula model" ] ] ] 7 => array:2 [ "identificador" => "sec1035" "titulo" => "Sampling bivariate algorithm" ] 8 => array:3 [ "identificador" => "sec2040" "titulo" => "Sampling trivariate algorithm" "secciones" => array:1 [ 0 => array:2 [ "identificador" => "sec0030" "titulo" => "Simulated Annealing method" ] ] ] 9 => array:3 [ "identificador" => "sec0035" "titulo" => "Case study" "secciones" => array:5 [ 0 => array:2 [ "identificador" => "sec0040" "titulo" => "Data description" ] 1 => array:2 [ "identificador" => "sec0045" "titulo" => "Statistical data analysis" ] 2 => array:2 [ "identificador" => "sec0050" "titulo" => "Bernstein Copula Model" ] 3 => array:2 [ "identificador" => "sec0055" "titulo" => "Sampling from Bernstein copula" ] 4 => array:2 [ "identificador" => "sec0060" "titulo" => "Permeability simulations and discussion" ] ] ] 10 => array:2 [ "identificador" => "sec0065" "titulo" => "Conclusions" ] 11 => array:2 [ "identificador" => "xack161016" "titulo" => "Acknowledgments" ] 12 => array:1 [ "titulo" => "References" ] ] ] "pdfFichero" => "main.pdf" "tienePdf" => true "fechaRecibido" => "2013-02-25" "fechaAceptado" => "2013-09-03" "PalabrasClave" => array:2 [ "es" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Palabras clave" "identificador" => "xpalclavsec519550" "palabras" => array:6 [ 0 => "permeabilidad" 1 => "porosidad" 2 => "velocidad de onda de corte" 3 => "dependencia multivariada" 4 => "cópula de Bernstein" 5 => "simulación geoestadística" ] ] ] "en" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Key words" "identificador" => "xpalclavsec519551" "palabras" => array:6 [ 0 => "permeability" 1 => "porosity" 2 => "shear wave velocity" 3 => "multivariate dependence" 4 => "Bernstein copula" 5 => "geostatistical simulation" ] ] ] ] "tieneResumen" => true "resumen" => array:2 [ "es" => array:2 [ "titulo" => "Resumen" "resumen" => "En este trabajo se presenta un método no paramétrico general de simulación estocástica conjunta de propiedades petrofísicas utilizando la cópula Bernstein. Este método consiste básicamente generar simulaciones estocásticas de una determinada propiedad petrofísica (variable primaria) modelando la dependencia empírica subyacente con otras propiedades petrofísicas (variables secundarias), mientras también es reproducida la dependencia espacial de la primera.Este enfoque multivariado provee una herramienta muy flexible para modelar las complejas relaciones de dependencia de las propiedades petrofísicas. Tiene varias ventajas sobre otros métodos tradicionales, ya que no se limita al caso de la dependencia lineal entre las variables, y tampoco requiere de la suposición de normalidad y/o existencia de momentos.En este trabajo este método es aplicado para simular un perfil de permeabilidad utilizando la porosidad vugular y velocidad de onda de corte (Ondas S) como covariables, en una formación carbonatada de doble porosidad a escala de pozo. Los valores simulados de la permeabilidad muestran un alto grado de precisión en comparación con los valores reales." ] "en" => array:2 [ "titulo" => "Abstract" "resumen" => "This paper introduces a general nonparametric method for joint stochastic simulation of petrophysical properties using the Bernstein copula. This method consists basically in generating stochastic simulations of a given petrophysical property (primary variable) modeling the underlying empirical dependence with other petrophysical properties (secondary variables) while reproducing the spatial dependence of the first one.This multivariate approach provides a very flexible tool to model the complex dependence relationships of petrophysical properties. It has several advantages over other traditional methods, since it is not restricted to the case of linear dependence among variables, it does not require the assumption of normality and/or existence of moments.In this paper this method is applied to simulate rock permeability using Vugular Porosity and Shear Wave Velocity (S-Waves) as covariates in a carbonate double-porosity formation at well log scale. Simulated permeability values show a high degree of accuracy compared to the actual values." ] ] "multimedia" => array:38 [ 0 => array:7 [ "identificador" => "f0005" "etiqueta" => "Figure 1" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr1.jpeg" "Alto" => 1239 "Ancho" => 878 "Tamanyo" => 138737 ] ] "descripcion" => array:1 [ "en" => "Diagram of the two-stage method to perform multivariate stochastic simulations using multivariate Bernstein copula." ] ] 1 => array:7 [ "identificador" => "f0010" "etiqueta" => "Figure 2" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr2.jpeg" "Alto" => 711 "Ancho" => 1417 "Tamanyo" => 152727 ] ] "descripcion" => array:1 [ "en" => "Crack porosity (PHICR), vuggy porosity (PHIV), total porosity (PHITOT) and permeability (K)." ] ] 2 => array:7 [ "identificador" => "f0015" "etiqueta" => "Figure 3" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr3.jpeg" "Alto" => 1078 "Ancho" => 1080 "Tamanyo" => 153450 ] ] "descripcion" => array:1 [ "en" => "Scatter plot distribution and histograms of PHIV and K of real data values." ] ] 3 => array:7 [ "identificador" => "f0020" "etiqueta" => "Figure 4" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr4.jpeg" "Alto" => 1068 "Ancho" => 1076 "Tamanyo" => 132907 ] ] "descripcion" => array:1 [ "en" => "Scatter plot distribution and histograms of PHIV and K of real data values." ] ] 4 => array:7 [ "identificador" => "f0025" "etiqueta" => "Figure 5" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr5.jpeg" "Alto" => 580 "Ancho" => 1788 "Tamanyo" => 164442 ] ] "descripcion" => array:1 [ "en" => "Scatter plot between Porosity and Permeability. Left, the 380 data values. Center, 380 simulated values using trivariate Bernstein copula. Right, 3800 simulated values." ] ] 5 => array:7 [ "identificador" => "f0030" "etiqueta" => "Figure 6" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr6.jpeg" "Alto" => 695 "Ancho" => 1441 "Tamanyo" => 110315 ] ] "descripcion" => array:1 [ "en" => "Single nonconditional simulation of permeability (K), using porosity and shear wave velocity (VS meas) as secondary variables." ] ] 6 => array:7 [ "identificador" => "f0035" "etiqueta" => "Figure 7" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr7.jpeg" "Alto" => 1073 "Ancho" => 1073 "Tamanyo" => 155056 ] ] "descripcion" => array:1 [ "en" => "Scatter plot and histograms of PHIV and K (Simulated values)." ] ] 7 => array:7 [ "identificador" => "f0040" "etiqueta" => "Figure 8" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr8.jpeg" "Alto" => 1081 "Ancho" => 1087 "Tamanyo" => 130538 ] ] "descripcion" => array:1 [ "en" => "Scatter plot and histograms of VS Meas and K (Simulated values)." ] ] 8 => array:7 [ "identificador" => "f0045" "etiqueta" => "Figure 9" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr9.jpeg" "Alto" => 1116 "Ancho" => 1113 "Tamanyo" => 105802 ] ] "descripcion" => array:1 [ "en" => "Estimated variogram and its fitted model for simulated K (Simulated values)." ] ] 9 => array:7 [ "identificador" => "f0050" "etiqueta" => "Figure 10" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr10.jpeg" "Alto" => 1118 "Ancho" => 1116 "Tamanyo" => 104997 ] ] "descripcion" => array:1 [ "en" => "Estimated variogram and its fitted model for original K (Data values)." ] ] 10 => array:7 [ "identificador" => "f0055" "etiqueta" => "Figure 11" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr11.jpeg" "Alto" => 679 "Ancho" => 1184 "Tamanyo" => 83653 ] ] "descripcion" => array:1 [ "en" => "Statistical comparison of to the original data and a single trivariate SA simulation." ] ] 11 => array:7 [ "identificador" => "f0060" "etiqueta" => "Figure 12" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr12.jpeg" "Alto" => 666 "Ancho" => 1449 "Tamanyo" => 104973 ] ] "descripcion" => array:1 [ "en" => "Median of 10 nonconditional simulation of permeability (K), using porosity and shear wave velocity (S-waves) as secondary variables." ] ] 12 => array:7 [ "identificador" => "f0065" "etiqueta" => "Figure 13" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr13.jpeg" "Alto" => 1080 "Ancho" => 1088 "Tamanyo" => 159213 ] ] "descripcion" => array:1 [ "en" => "Scatter plot and histograms of PHIV and K simulated values (median of 10 simulations)." ] ] 13 => array:7 [ "identificador" => "f0070" "etiqueta" => "Figure 14" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr14.jpeg" "Alto" => 1083 "Ancho" => 1085 "Tamanyo" => 133757 ] ] "descripcion" => array:1 [ "en" => "Scatter plot distribution and histograms of VS-meas and K simulated values (a median of 10 simulations)." ] ] 14 => array:7 [ "identificador" => "f0075" "etiqueta" => "Figure 15" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr15.jpeg" "Alto" => 1113 "Ancho" => 1107 "Tamanyo" => 107974 ] ] "descripcion" => array:1 [ "en" => "Estimated variogram and its fitted model for simulated K data values." ] ] 15 => array:7 [ "identificador" => "f0080" "etiqueta" => "Figure 16" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr16.jpeg" "Alto" => 648 "Ancho" => 1179 "Tamanyo" => 81663 ] ] "descripcion" => array:1 [ "en" => "Statistical comparison of to the original data, and the median of 10 SA simulations." ] ] 16 => array:7 [ "identificador" => "f0085" "etiqueta" => "Figure 17" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr17.jpeg" "Alto" => 1671 "Ancho" => 1449 "Tamanyo" => 269539 ] ] "descripcion" => array:1 [ "en" => "Spatial distribution of permeability in single conditional simulations using PHIV and VS meas as secondary variables. Conditioning levels (in top to bottom order) 10 percent, 50 percent and 90 percent." ] ] 17 => array:7 [ "identificador" => "tbl0005" "etiqueta" => "Table 1" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Explanatory variable \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Corresponding Hoeffding index Φ vs. K \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHIV \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.71 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHITOT \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.62 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHIMAT \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.50 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">VS Meas \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.60 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">VP meas \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.55 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHICR \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.16 \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab795073.png" ] ] ] ] "descripcion" => array:1 [ "en" => "Dependence index Φ, <a class="elsevierStyleCrossRef" href="#bib0115">Hoeffding (1941)</a>, for possible explanatory variables (PHITOT, PHIMAT, PHICR, PHIV, VS Meas, VP meas)." ] ] 18 => array:7 [ "identificador" => "tbl0010" "etiqueta" => "Table 2" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Hoeffding index Φ \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHIV \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHITOT \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHIMAT \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHICR \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">VS Meas \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">VP meas \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHIV \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.88 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.72 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.34 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.58 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.57 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHITOT \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.43 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.32 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.64 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.61 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHIMAT \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.24 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.30 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.36 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">PHICR \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.31 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.20 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">VS Meas \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.32 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">VP meas \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">1 \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab795072.png" ] ] ] ] "descripcion" => array:1 [ "en" => "Hoeffding dependence index Φ matrix, for possible explanatory variables (PHITOT, PHIMAT, PHICR, PHIV, VS Meas, VP meas)." ] ] 19 => array:7 [ "identificador" => "tbl0015" "etiqueta" => "Table 3" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Nugget \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Sill \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Range \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">Data values \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.00 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">2.61 x 106 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">85.07 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">Single SA nonconditional Simulation. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.00 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">2.57 x 106 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">87.52 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="center" valign="top">Median of 10 SA nonconditional Simulations. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.00 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">2.51 x 106 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">87.27 \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab795070.png" ] ] ] ] "descripcion" => array:1 [ "en" => "Variograms Models of the data values and single nonconditional SA simulation." ] ] 20 => array:7 [ "identificador" => "tbl0020" "etiqueta" => "Table 4" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">MSE \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Single SA nonconditional Sim. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.50 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Median of 10 nonconditionals Sims. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.35 \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab795071.png" ] ] ] ] "descripcion" => array:1 [ "en" => "MSE values of nonconditional simulations" ] ] 21 => array:7 [ "identificador" => "tbl0025" "etiqueta" => "Table 5" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">MSE \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Single SA nonconditional Sim. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.50 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Single SA 10% conditional Sim. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.41 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Single SA 50% conditional Sim. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.36 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Single SA 90% conditional Sim. \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="center" valign="top">0.27 \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab795069.png" ] ] ] ] "descripcion" => array:1 [ "en" => "MSE values of the single conditional simulations." ] ] 22 => array:6 [ "identificador" => "eq0005" "etiqueta" => "(1)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "H(x1,x2,…,xm)=C(F2(x1),F2(x),…,Fm(xm))" "Fichero" => "si1.jpeg" "Tamanyo" => 2642 "Alto" => 14 "Ancho" => 338 ] ] 23 => array:6 [ "identificador" => "eq0010" "etiqueta" => "(2)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "Fˆj(x)=1n∑k=1n1{xjk≤x}" "Fichero" => "si2.jpeg" "Tamanyo" => 2047 "Alto" => 46 "Ancho" => 160 ] ] 24 => array:6 [ "identificador" => "eq0015" "etiqueta" => "(3)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "Cn(v1n,…,vmn)=1nk=1n1{rank(x1k)£v1,…,rank(xmk)≤vm}" "Fichero" => "si6.jpeg" "Tamanyo" => 3587 "Alto" => 22 "Ancho" => 414 ] ] 25 => array:6 [ "identificador" => "eq0020" "etiqueta" => "(4)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "H(x1,x2,x3)=C(F1(x1),F2(x2),F3(x3))" "Fichero" => "si7.jpeg" "Tamanyo" => 2437 "Alto" => 14 "Ancho" => 276 ] ] 26 => array:6 [ "identificador" => "eq0025" "etiqueta" => "(5)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "Q˜;n(u)=∑k=0n12(xk+xk+1)nkuk(1-u)n-k" "Fichero" => "si9.jpeg" "Tamanyo" => 3586 "Alto" => 49 "Ancho" => 301 ] ] 27 => array:6 [ "identificador" => "eq0030" "etiqueta" => "(6)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "C˜;B(u1,…,um)=∑v1=0n1…∑vm=0nmαv1n1,…,vmnmPv1,m1(u1)…Pvm,nm(um)" "Fichero" => "si10.jpeg" "Tamanyo" => 4854 "Alto" => 49 "Ancho" => 484 ] ] 28 => array:6 [ "identificador" => "eq0035" "etiqueta" => "(7)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "Pvj,mj(uj)=njvjujvj(1-uj)nj-vj" "Fichero" => "si11.jpeg" "Tamanyo" => 2735 "Alto" => 49 "Ancho" => 243 ] ] 29 => array:6 [ "identificador" => "eq0040" "etiqueta" => "(8)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "C˜;B(u,v,w)=∑i=0n∑j=0n∑k=0nCnin,jn,knniui(1-u)n-injvj(1-v)n-jnkwk(1-w)n-k" "Fichero" => "si13.jpeg" "Tamanyo" => 8387 "Alto" => 52 "Ancho" => 647 ] ] 30 => array:6 [ "identificador" => "eq0045" "etiqueta" => "(9)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "cu(v)=∂C˜;B(u,v)∂u" "Fichero" => "si14.jpeg" "Tamanyo" => 1178 "Alto" => 24 "Ancho" => 106 ] ] 31 => array:6 [ "identificador" => "eq0050" "etiqueta" => "(10)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "cuv(w)=∂C˜;B(u,v,w)∂u∂v∂C˜;B(u,v,1)∂u∂v" "Fichero" => "si19.jpeg" "Tamanyo" => 1884 "Alto" => 42 "Ancho" => 122 ] ] 32 => array:6 [ "identificador" => "eq0055" "etiqueta" => "(11)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "FO=∑i[γ*(hi)-γ(hi)]2γ(hi)2" "Fichero" => "si23.jpeg" "Tamanyo" => 1911 "Alto" => 40 "Ancho" => 151 ] ] 33 => array:6 [ "identificador" => "eq0060" "etiqueta" => "(12)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "γ(h)=12N(h)∑i=1N(h)[Z(xi+h)-Z(xi)]2" "Fichero" => "si24.jpeg" "Tamanyo" => 2790 "Alto" => 50 "Ancho" => 248 ] ] 34 => array:6 [ "identificador" => "eq0065" "etiqueta" => "(13)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "T0=(-ΔElog(τ0))" "Fichero" => "si25.jpeg" "Tamanyo" => 1078 "Alto" => 23 "Ancho" => 93 ] ] 35 => array:6 [ "identificador" => "eq0070" "etiqueta" => "(14)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "MSE(logK|PHIV)=1n∑m=1n[zm-β0.5(xm)]2" "Fichero" => "si27.jpeg" "Tamanyo" => 3414 "Alto" => 46 "Ancho" => 311 ] ] 36 => array:6 [ "identificador" => "eq0075" "etiqueta" => "(15)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "MSE(logK|VSmeas)=1n∑m=1n[zm-β0.5(ym)]2" "Fichero" => "si28.jpeg" "Tamanyo" => 3613 "Alto" => 46 "Ancho" => 334 ] ] 37 => array:6 [ "identificador" => "eq0080" "etiqueta" => "(16)" "tipo" => "MULTIMEDIAFORMULA" "mostrarFloat" => false "mostrarDisplay" => true "Formula" => array:5 [ "Matematica" => "MSE(logK|PHIV,VSmeas)=1n∑m=1n[zm-β0.5(xm,ym)]2" "Fichero" => "si29.jpeg" "Tamanyo" => 4305 "Alto" => 46 "Ancho" => 399 ] ] ] "bibliografia" => array:2 [ "titulo" => "References" "seccion" => array:1 [ 0 => array:2 [ "identificador" => "bibs0005" "bibliografiaReferencia" => array:38 [ 0 => array:3 [ "identificador" => "bib0005" "etiqueta" => "Al-harthy et al., 2005" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Copulas: A new technique to model dependence in petroleum decision making" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "M. Al-harthy" 1 => "S. Begg" 2 => "R. Bratvold" ] ] ] ] ] "host" => array:1 [ 0 => array:2 [ "doi" => "10.1038/srep09602" "Revista" => array:6 [ "tituloSerie" => "Journal of Petroleum Science & Engineering" "fecha" => "2005" "volumen" => "57" "paginaInicial" => "195" "paginaFinal" => "208" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319" "web" => "Medline" ] ] ] ] ] ] ] ] 1 => array:3 [ "identificador" => "bib0010" "etiqueta" => "Armstrong et al., 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Incorporating technical uncertainty in real option valuation of oil projects" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "M. Armstrong" 1 => "A. Galli" 2 => "W. Bailey" 3 => "B. Couet" ] ] ] ] ] "host" => array:1 [ 0 => array:2 [ "doi" => "10.1038/srep09602" "Revista" => array:6 [ "tituloSerie" => "Journal of Petroleum Science & Engineering" "fecha" => "2004" "volumen" => "44" "paginaInicial" => "67" "paginaFinal" => "82" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319" "web" => "Medline" ] ] ] ] ] ] ] ] 2 => array:3 [ "identificador" => "bib0015" "etiqueta" => "Balan et al., 1995" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "State-of-The-Art in Permeability Determination From Well Log Data: Part 1.- A Comparative Study" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "B. Balan" 1 => "S. Mohaghegh" 2 => "S. Ameri" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "Model Development" "fecha" => "1995" ] ] ] ] ] ] 3 => array:3 [ "identificador" => "bib0020" "etiqueta" => "Bennett et al., 2002" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Hydrologic investigation of the Floridian aquifer system, western Hillsboro Basin, Palm Beach County, Florida, Technical Publication WS-8" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "W.M. Bennett" 1 => "P.F. Linton" 2 => "E.E. Rectenwald" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "South Florida Water Management District" "fecha" => "2002" ] ] ] ] ] ] 4 => array:3 [ "identificador" => "bib0025" "etiqueta" => "Billingsley, 1995" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "P. Billingsley" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Probability and Measure" "fecha" => "1995" "editorial" => "Wiley" "editorialLocalizacion" => "New York" ] ] ] ] ] ] 5 => array:3 [ "identificador" => "bib0030" "etiqueta" => "Cosentino, 2001" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "L. Cosentino" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "Integrated Reservoir Studies" "fecha" => "2001" ] ] ] ] ] ] 6 => array:3 [ "identificador" => "bib0035" "etiqueta" => "Dafflon et al., 2009a" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "B. Dafflon" 1 => "J. Irving" 2 => "K. Holliger" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Applied Geophysics" "fecha" => "2009" "volumen" => "68" "paginaInicial" => "60" "paginaFinal" => "70" ] ] ] ] ] ] 7 => array:3 [ "identificador" => "bib0040" "etiqueta" => "Dafflon et al., 2009b" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Quantitative Integration of High-Resolution Hydrogeophysical Data: A Novel Approach to Monte-Carlo-Type Conditional Stochastic Simulations and Implications for Hydrological Predictions" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "B. Dafflon" 1 => "J. Irving" 2 => "K. Holliger*" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Earth Science" "fecha" => "2009" "volumen" => "20" "paginaInicial" => "580" "paginaFinal" => "591" ] ] ] ] ] ] 8 => array:3 [ "identificador" => "bib0045" "etiqueta" => "Dafflon and Barrash, 2012" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Three-dimensional stochastic estimation of porosity distribution: Benefits of using ground-penetrating radar velocity tomograms in simulated-annealing-based or Bayesian sequential simulation approaches" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "B. Dafflon" 1 => "W. Barrash" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:4 [ "tituloSerie" => "Water Resources Research" "fecha" => "2012" "volumen" => "48" "paginaInicial" => "5" ] ] ] ] ] ] 9 => array:3 [ "identificador" => "bib0050" "etiqueta" => "Accioly and Chiyoshi, 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Modeling dependence with copulas: a useful tool for field development decision process" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "R. Accioly" 1 => "F. Chiyoshi" ] ] ] ] ] "host" => array:1 [ 0 => array:2 [ "doi" => "10.1038/srep09602" "Revista" => array:6 [ "tituloSerie" => "Journal of Petroleum Science & Engineering" "fecha" => "2004" "volumen" => "44" "paginaInicial" => "83" "paginaFinal" => "91" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319" "web" => "Medline" ] ] ] ] ] ] ] ] 10 => array:3 [ "identificador" => "bib0055" "etiqueta" => "Deheuvels, 1979" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "La fonction de dépendance empirique et ses porpriétés. Un test nonparamétrique d’independance, Acad" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "P. Deheuvels" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:6 [ "tituloSerie" => "Roy. Belg. Bull. Cl. Sci. (5)" "fecha" => "1979" "volumen" => "65" "numero" => "6" "paginaInicial" => "274" "paginaFinal" => "292" ] ] ] ] ] ] 11 => array:3 [ "identificador" => "bib0060" "etiqueta" => "Deutsch, 1992" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Annealing Techniques Applied to Reservoir Modeling and the Integration of Geological and Engineering (Well Test) Data" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "C.V. Deutsch" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "PhD thesis" "fecha" => "1992" "editorial" => "Stanford University" "editorialLocalizacion" => "Stanford, CA" ] ] ] ] ] ] 12 => array:3 [ "identificador" => "bib0065" "etiqueta" => "Deutsch and Cockerham, 1994" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "C.V. Deutsch" 1 => "P.W. Cockerham" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "Geostatistical Modeling of Permeability with Annealing Cosimulation (ACS)" "fecha" => "1994" ] ] ] ] ] ] 13 => array:3 [ "identificador" => "bib0070" "etiqueta" => "Díaz-Viera and Casar-González, 2005" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Stochastic simulation of complex dependency patterns of petrophysical properties using t-copulas" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "M. Díaz-Viera" 1 => "R. Casar-González" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Proc. IAMG’05: GIS and Spatial Analysis" "fecha" => "2005" "volumen" => "2" "paginaInicial" => "749" "paginaFinal" => "755" ] ] ] ] ] ] 14 => array:3 [ "identificador" => "bib0075" "etiqueta" => "Díaz-Viera et al., 2006" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Stochastic modeling of permeability in double porosity carbonates applying a Monte-Carlo simulation method with t-copula" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:5 [ 0 => "M. Díaz-Viera" 1 => "P. Anguiano-Rojas" 2 => "A. Mousatov" 3 => "E. Kazatchenko" 4 => "M. Markov" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "SPWLA 47th Annual Logging Symposium" "fecha" => "2006" ] ] ] ] ] ] 15 => array:3 [ "identificador" => "bib0080" "etiqueta" => "Dréo et al., 2006" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "J. Dréo" 1 => "A. Pétrowski" 2 => "P. Siarry" 3 => "T. Taillard" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Metaheuristics for Hard Optimization" "fecha" => "2006" "editorial" => "Springer-Verlag Berlin Heidelberg" "editorialLocalizacion" => "Germany" ] ] ] ] ] ] 16 => array:3 [ "identificador" => "bib0085" "etiqueta" => "Erdely and Díaz-Viera, 2009" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Nonparametric and semiparametric bivariate modeling of petrophysical porosity-permeability well log data" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "A. Erdely" 1 => "M. Díaz-Viera" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Copula Theory and Its Applications (P Jaworski et al. eds), Lecture Notes in Statistics – Springer" "fecha" => "2009" "volumen" => "198" "paginaInicial" => "267" "paginaFinal" => "278" ] ] ] ] ] ] 17 => array:3 [ "identificador" => "bib0090" "etiqueta" => "Embrechts et al., 1999" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Correlation: pitfalls and alternatives" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "P. Embrechts" 1 => "A. McNeil" 2 => "D. Straumann" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Risk Magazine" "fecha" => "1999" "volumen" => "5" "paginaInicial" => "69" "paginaFinal" => "71" ] ] ] ] ] ] 18 => array:3 [ "identificador" => "bib0095" "etiqueta" => "Embrechts et al., 2003" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Modeling dependence with copulas and applications to risk management" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "P. Embrechts" 1 => "F. Lindskog" 2 => "A.J. McNeil" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "LibroEditado" => array:4 [ "titulo" => "Handbook of Heavy-Tailed Distributions in Finance" "paginaInicial" => "329" "paginaFinal" => "384" "serieFecha" => "2003" ] ] ] ] ] ] 19 => array:3 [ "identificador" => "bib0100" "etiqueta" => "Fermanian et al., 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Weak convergence of empirical copula processes" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "J-D. Fermanian" 1 => "D. Radulovic" 2 => "M. Wegcamp" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Bernoulli" "fecha" => "2004" "volumen" => "10" "paginaInicial" => "547" "paginaFinal" => "560" ] ] ] ] ] ] 20 => array:3 [ "identificador" => "bib0105" "etiqueta" => "Haldorsen and Damsleth, 1990" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Stochastic modeling" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "H. Haldorsen" 1 => "E. Damsleth" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Journal of Petroleum Technology" "fecha" => "1990" "paginaInicial" => "404" "paginaFinal" => "412" ] ] ] ] ] ] 21 => array:3 [ "identificador" => "bib0110" "etiqueta" => "Hernández-Maldonado et al., 2012" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "A joint stochastic simulation method using the Bernstein copula as a flexible tool for modeling nonlinear dependence structures between petrophysical properties" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "V. Hernández-Maldonado" 1 => "M. Díaz-Viera" 2 => "A. Erdely" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Journal of Petroleum Science & Engineering" "fecha" => "2012" "paginaInicial" => "90" "paginaFinal" => "91" ] ] ] ] ] ] 22 => array:3 [ "identificador" => "bib0115" "etiqueta" => "Hoeffding, 1941" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Scale-invariant correlation theory" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "W. Hoeffding" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "LibroEditado" => array:4 [ "titulo" => "The Collected Works of Wassily Hoeffding" "paginaInicial" => "57" "paginaFinal" => "107" "serieFecha" => "1941" ] ] ] ] ] ] 23 => array:3 [ "identificador" => "bib0120" "etiqueta" => "Holden et al., 1995" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Use of Well Test Data in Stochastic Reservoir Modelling" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:6 [ 0 => "L. Holden" 1 => "R. Madsen" 2 => "A. Skorstad" 3 => "K.A. Jakobsen" 4 => "C.B. Tjolsen" 5 => "S. Vik" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "SPE annual Technical Conference and Exhibition held in Dallas" "fecha" => "1995" "paginaInicial" => "22" "paginaFinal" => "25" ] ] ] ] ] ] 24 => array:3 [ "identificador" => "bib0125" "etiqueta" => "Journel and Alabert, 1990" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "New method for reservoir mapping" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "A. Journel" 1 => "F. Alabert" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Journal of Petroleum Technology" "fecha" => "1990" "paginaInicial" => "212" "paginaFinal" => "218" ] ] ] ] ] ] 25 => array:3 [ "identificador" => "bib0130" "etiqueta" => "Journel and Zhu, 1990" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Integrating soft seismic data: Markov-Bayes updating, an alternative to cokriging and traditional regression" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "A. Journel" 1 => "H. Zhu" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:2 [ "titulo" => "In Report 3, Stanford Center for Reservoir Forecasting" "fecha" => "1990" ] ] ] ] ] ] 26 => array:3 [ "identificador" => "bib0135" "etiqueta" => "Kazatchenko et al., 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Joint inversion of acoustic and resistivity data for carbonate microstructure evaluation" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "E. Kazatchenko" 1 => "M. Markov" 2 => "A. Mousatov" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Petrophysics" "fecha" => "2004" "volumen" => "45" "paginaInicial" => "130" "paginaFinal" => "140" ] ] ] ] ] ] 27 => array:3 [ "identificador" => "bib0140" "etiqueta" => "Kazatchenko et al., 2006" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Carbonate microstructure determination by inversion of acoustic and electrical data: application to a South Florida Aquifer" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "E. Kazatchenko" 1 => "M. Markov" 2 => "A. Mousatov" 3 => "J. Parra" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "J. Appl. Geophys" "fecha" => "2006" "volumen" => "59" "paginaInicial" => "1" "paginaFinal" => "15" ] ] ] ] ] ] 28 => array:3 [ "identificador" => "bib0145" "etiqueta" => "Landa et al., 1996" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Reservoir Characterization Constrained to Well Test Data: A Field Example" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "J.L. Landa" 1 => "M.M. Kamal" 2 => "C.D. Jenkins" 3 => "R.N. Horne" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "SPE Annual Technical Conference and Exhibition held in Denver" "fecha" => "1996" "paginaInicial" => "6" "paginaFinal" => "9" ] ] ] ] ] ] 29 => array:3 [ "identificador" => "bib0150" "etiqueta" => "Muñoz-Perez and Fernández-Palacín, 1987" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Estimating the quantile function by Bernstein polynomials" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "J. Muñoz-Perez" 1 => "A. Fernández-Palacín" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Computational Statistics and Data Analysis" "fecha" => "1987" "volumen" => "5" "paginaInicial" => "391" "paginaFinal" => "397" ] ] ] ] ] ] 30 => array:3 [ "identificador" => "bib0155" "etiqueta" => "Nelsen, 2006" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "R.B. Nelsen" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "An Introduction to Copulas" "fecha" => "2006" "editorial" => "Springer" "editorialLocalizacion" => "New York" ] ] ] ] ] ] 31 => array:3 [ "identificador" => "bib0160" "etiqueta" => "Parra et al., 2001" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "A methodology to integrate magnetic resonance and acoustic measurements for reservoir characterization" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "J.O. Parra" 1 => "C.L. Hackert" 2 => "H.A. Collier" 3 => "M. Bennett" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "Report DOE/BC/15203-3" "fecha" => "2001" "editorial" => "National Petroleum Technology Office, Department of Energy" "editorialLocalizacion" => "Tulsa, Oklahoma" ] ] ] ] ] ] 32 => array:3 [ "identificador" => "bib0165" "etiqueta" => "Parra and Hackert, 2002" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Permeability and porosity images based on crosswell reflection seismic measurements of a vuggy carbonate aquifer at the Hillsboro site" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "J.O. Parra" 1 => "C.L. Hackert" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:3 [ "titulo" => "South Florida 72nd Annual Meeting of the Society of Exploration Geophysicists" "fecha" => "2002" "editorial" => "Paper VCD P1.2" ] ] ] ] ] ] 33 => array:3 [ "identificador" => "bib0170" "etiqueta" => "Sancetta and Satchell, 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "The Bernstein copula and its applications to modeling and approximations of multivariate distributions" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "A. Sancetta" 1 => "S. Satchell" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Econometric Theory" "fecha" => "2004" "volumen" => "20" "paginaInicial" => "535" "paginaFinal" => "562" ] ] ] ] ] ] 34 => array:3 [ "identificador" => "bib0175" "etiqueta" => "Sancetta, 2007" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Nonparametric estimation of distributions with given marginals via Bernstein-Kantorovic polynomials: L1 and pointwise convergence theory" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "A. Sancetta" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:6 [ "tituloSerie" => "Journal of Multivariate Analysis" "fecha" => "2007" "volumen" => "98" "numero" => "7" "paginaInicial" => "1376" "paginaFinal" => "1390" ] ] ] ] ] ] 35 => array:3 [ "identificador" => "bib0180" "etiqueta" => "Siburg and Stoimenov, 2008" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "A scalar product for copulas" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "K. Siburg" 1 => "P. Stoimenov" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Mathematical Analysis and Applications" "fecha" => "2008" "volumen" => "344" "paginaInicial" => "429" "paginaFinal" => "439" ] ] ] ] ] ] 36 => array:3 [ "identificador" => "bib0185" "etiqueta" => "Suro-Perez, 1992" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Indicator Kriging and Simulation Based on Principal Component Analysis" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "V. Suro-Perez" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:4 [ "titulo" => "PhD thesis" "fecha" => "1992" "editorial" => "Stanford University" "editorialLocalizacion" => "Stanford, CA" ] ] ] ] ] ] 37 => array:3 [ "identificador" => "bib0190" "etiqueta" => "Sklar, 1959" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:1 [ "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "A. Sklar" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "LibroEditado" => array:5 [ "paginaInicial" => "229" "paginaFinal" => "231" "serieVolumen" => "8" "serieTitulo" => "Fonctions de répartition á n dimensions et leurs marges" "serieFecha" => "1959" ] ] ] ] ] ] ] ] ] ] "agradecimientos" => array:1 [ 0 => array:4 [ "identificador" => "xack161016" "titulo" => "Acknowledgments" "texto" => "The present work was partially supported by the two following projects: IN110311 from Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (PAPIIT) at Universidad Nacional Autónoma de México, and by 116606 (Y.00102) from Fondo Sectorial CONACYT-SENER-Hidrocarburos." "vista" => "all" ] ] ] "idiomaDefecto" => "en" "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714989/v2_201505081358/en/main.assets" "Apartado" => array:4 [ "identificador" => "36047" "tipo" => "SECCION" "es" => array:2 [ "titulo" => "Original paper" "idiomaDefecto" => true ] "idiomaDefecto" => "es" ] "PDF" => "https://static.elsevier.es/multimedia/00167169/0000005300000002/v2_201505081358/S0016716914714989/v2_201505081358/en/main.pdf?idApp=UINPBA00004N&text.app=https://www.elsevier.es/" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714989?idApp=UINPBA00004N"]

Statistics

Victor Hernández-Maldonado

Corresponding author

vmhernann@yahoo.com.mx

Corresponding author:

, Martín Díaz-Viera

Instituto Mexicano del Petróleo, Eje Central Lázaro Cárdenas 152, 07730, México, D.F.

Arturo Erdely

Programa de Actuaria, Facultad de Estudios Superiores Acatlán, Universidad Nacional Autónoma de México, México D.F.

Read

2461

Times
was read the article

765

Total PDF

1696

Total HTML

Share statistics

 array:23 [
  "pii" => "S0016716914714989"
  "issn" => "00167169"
  "doi" => "10.1016/S0016-7169(14)71498-9"
  "estado" => "S300"
  "fechaPublicacion" => "2014-04-01"
  "aid" => "71498"
  "copyright" => "Universidad Nacional Aut&#243;noma de M&#233;xico"
  "copyrightAnyo" => "2014"
  "documento" => "article"
  "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/"
  "subdocumento" => "fla"
  "cita" => "Geofisica Internacional. 2014;53:163-81"
  "abierto" => array:3 [
    "ES" => true
    "ES2" => true
    "LATM" => true
  ]
  "gratuito" => true
  "lecturas" => array:2 [
    "total" => 1298
    "formatos" => array:3 [
      "EPUB" => 25
      "HTML" => 843
      "PDF" => 430
    ]
  ]
  "itemSiguiente" => array:18 [
    "pii" => "S0016716914714990"
    "issn" => "00167169"
    "doi" => "10.1016/S0016-7169(14)71499-0"
    "estado" => "S300"
    "fechaPublicacion" => "2014-04-01"
    "aid" => "71499"
    "copyright" => "Universidad Nacional Aut&#243;noma de M&#233;xico"
    "documento" => "article"
    "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/"
    "subdocumento" => "fla"
    "cita" => "Geofisica Internacional. 2014;53:183-98"
    "abierto" => array:3 [
      "ES" => true
      "ES2" => true
      "LATM" => true
    ]
    "gratuito" => true
    "lecturas" => array:2 [
      "total" => 1383
      "formatos" => array:3 [
        "EPUB" => 39
        "HTML" => 936
        "PDF" => 408
      ]
    ]
    "es" => array:11 [
      "idiomaDefecto" => true
      "titulo" => "Joint interpretation of geoelectrical and volatile organic compounds data&#58; a case study in a hydrocarbons contaminated urban site"
      "tienePdf" => "es"
      "tieneTextoCompleto" => "es"
      "tieneResumen" => array:2 [
        0 => "es"
        1 => "en"
      ]
      "paginas" => array:1 [
        0 => array:2 [
          "paginaInicial" => "183"
          "paginaFinal" => "198"
        ]
      ]
      "contieneResumen" => array:2 [
        "es" => true
        "en" => true
      ]
      "contieneTextoCompleto" => array:1 [
        "es" => true
      ]
      "contienePdf" => array:1 [
        "es" => true
      ]
      "resumenGrafico" => array:2 [
        "original" => 0
        "multimedia" => array:7 [
          "identificador" => "f0045"
          "etiqueta" => "Figure 9"
          "tipo" => "MULTIMEDIAFIGURA"
          "mostrarFloat" => true
          "mostrarDisplay" => false
          "figura" => array:1 [
            0 => array:4 [
              "imagen" => "gr9.jpeg"
              "Alto" => 854
              "Ancho" => 1221
              "Tamanyo" => 74157
            ]
          ]
          "descripcion" => array:1 [
            "es" => "<p id="sp0045" class="elsevierStyleSimplePara elsevierViewall">BTEX ratio of free-phase gasoline samples&#46;</p>"
          ]
        ]
      ]
      "autores" => array:2 [
        0 => array:2 [
          "autoresLista" => "Omar Delgado-Rodr&#237;guez, David Flores-Hern&#225;ndez, Myriam A&#46; Amezcua-Allieri, Andr&#233;s Rosas-Molina, Salvador Mar&#237;n-C&#243;rdova"
          "autores" => array:5 [
            0 => array:2 [
              "nombre" => "Omar"
              "apellidos" => "Delgado-Rodr&#237;guez"
            ]
            1 => array:2 [
              "nombre" => "David"
              "apellidos" => "Flores-Hern&#225;ndez"
            ]
            2 => array:2 [
              "nombre" => "Myriam A&#46;"
              "apellidos" => "Amezcua-Allieri"
            ]
            3 => array:2 [
              "nombre" => "Andr&#233;s"
              "apellidos" => "Rosas-Molina"
            ]
            4 => array:2 [
              "nombre" => "Salvador"
              "apellidos" => "Mar&#237;n-C&#243;rdova"
            ]
          ]
        ]
        1 => array:2 [
          "autoresLista" => "Vladimir Shevnin"
          "autores" => array:1 [
            0 => array:2 [
              "nombre" => "Vladimir"
              "apellidos" => "Shevnin"
            ]
          ]
        ]
      ]
    ]
    "idiomaDefecto" => "es"
    "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714990?idApp=UINPBA00004N"
    "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714990/v2_201505081358/es/main.assets"
  ]
  "itemAnterior" => array:18 [
    "pii" => "S0016716914714977"
    "issn" => "00167169"
    "doi" => "10.1016/S0016-7169(14)71497-7"
    "estado" => "S300"
    "fechaPublicacion" => "2014-04-01"
    "aid" => "71497"
    "copyright" => "Universidad Nacional Aut&#243;noma de M&#233;xico"
    "documento" => "article"
    "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/"
    "subdocumento" => "fla"
    "cita" => "Geofisica Internacional. 2014;53:153-62"
    "abierto" => array:3 [
      "ES" => true
      "ES2" => true
      "LATM" => true
    ]
    "gratuito" => true
    "lecturas" => array:2 [
      "total" => 825
      "formatos" => array:3 [
        "EPUB" => 24
        "HTML" => 507
        "PDF" => 294
      ]
    ]
    "en" => array:11 [
      "idiomaDefecto" => true
      "titulo" => "Sensitivity of the surface temperature to changes in total solar irradiance calculated with the WRF model"
      "tienePdf" => "en"
      "tieneTextoCompleto" => "en"
      "tieneResumen" => array:2 [
        0 => "es"
        1 => "en"
      ]
      "paginas" => array:1 [
        0 => array:2 [
          "paginaInicial" => "153"
          "paginaFinal" => "162"
        ]
      ]
      "contieneResumen" => array:2 [
        "es" => true
        "en" => true
      ]
      "contieneTextoCompleto" => array:1 [
        "en" => true
      ]
      "contienePdf" => array:1 [
        "en" => true
      ]
      "resumenGrafico" => array:2 [
        "original" => 0
        "multimedia" => array:7 [
          "identificador" => "f0035"
          "etiqueta" => "Figure 7"
          "tipo" => "MULTIMEDIAFIGURA"
          "mostrarFloat" => true
          "mostrarDisplay" => false
          "figura" => array:1 [
            0 => array:4 [
              "imagen" => "gr7.jpeg"
              "Alto" => 1022
              "Ancho" => 1879
              "Tamanyo" => 195044
            ]
          ]
          "descripcion" => array:1 [
            "en" => "<p id="sp0035" class="elsevierStyleSimplePara elsevierViewall">Daily temperature sensitivity in&#58; a&#41; The whole area of study&#46; b&#41; Zone A&#46; c&#41; Zone B&#46; d&#41; Zone C&#46; e&#41; Zone D&#46; f&#41; Zone E&#46;</p>"
          ]
        ]
      ]
      "autores" => array:3 [
        0 => array:2 [
          "autoresLista" => "Carolina Cipagauta"
          "autores" => array:1 [
            0 => array:2 [
              "nombre" => "Carolina"
              "apellidos" => "Cipagauta"
            ]
          ]
        ]
        1 => array:2 [
          "autoresLista" => "Blanca Mendoza"
          "autores" => array:1 [
            0 => array:2 [
              "nombre" => "Blanca"
              "apellidos" => "Mendoza"
            ]
          ]
        ]
        2 => array:2 [
          "autoresLista" => "Jorge Zavala-Hidalgo"
          "autores" => array:1 [
            0 => array:2 [
              "nombre" => "Jorge"
              "apellidos" => "Zavala-Hidalgo"
            ]
          ]
        ]
      ]
    ]
    "idiomaDefecto" => "en"
    "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714977?idApp=UINPBA00004N"
    "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714977/v2_201505081358/en/main.assets"
  ]
  "en" => array:18 [
    "idiomaDefecto" => true
    "titulo" => "A multivariate Bernstein copula model for permeability stochastic simulation"
    "tieneTextoCompleto" => true
    "paginas" => array:1 [
      0 => array:2 [
        "paginaInicial" => "163"
        "paginaFinal" => "181"
      ]
    ]
    "autores" => array:2 [
      0 => array:4 [
        "autoresLista" => "Victor Hern&#225;ndez-Maldonado, Mart&#237;n D&#237;az-Viera"
        "autores" => array:2 [
          0 => array:4 [
            "nombre" => "Victor"
            "apellidos" => "Hern&#225;ndez-Maldonado"
            "email" => array:1 [
              0 => "vmhernann&#64;yahoo&#46;com&#46;mx"
            ]
            "referencia" => array:1 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">&#42;</span>"
                "identificador" => "cor0005"
              ]
            ]
          ]
          1 => array:2 [
            "nombre" => "Mart&#237;n"
            "apellidos" => "D&#237;az-Viera"
          ]
        ]
        "afiliaciones" => array:1 [
          0 => array:2 [
            "entidad" => "Instituto Mexicano del Petr&#243;leo&#44; Eje Central L&#225;zaro C&#225;rdenas 152&#44; 07730&#44; M&#233;xico&#44; D&#46;F&#46;"
            "identificador" => "aff0005"
          ]
        ]
        "correspondencia" => array:1 [
          0 => array:3 [
            "identificador" => "cor0005"
            "etiqueta" => "&#42;"
            "correspondencia" => "Corresponding author&#58;"
          ]
        ]
      ]
      1 => array:3 [
        "autoresLista" => "Arturo Erdely"
        "autores" => array:1 [
          0 => array:2 [
            "nombre" => "Arturo"
            "apellidos" => "Erdely"
          ]
        ]
        "afiliaciones" => array:1 [
          0 => array:2 [
            "entidad" => "Programa de Actuaria&#44; Facultad de Estudios Superiores Acatl&#225;n&#44; Universidad Nacional Aut&#243;noma de M&#233;xico&#44; M&#233;xico D&#46;F&#46;"
            "identificador" => "aff0010"
          ]
        ]
      ]
    ]
    "resumenGrafico" => array:2 [
      "original" => 0
      "multimedia" => array:7 [
        "identificador" => "f0040"
        "etiqueta" => "Figure 8"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr8.jpeg"
            "Alto" => 1081
            "Ancho" => 1087
            "Tamanyo" => 130538
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0040" class="elsevierStyleSimplePara elsevierViewall">Scatter plot and histograms of VS Meas and K &#40;Simulated values&#41;&#46;</p>"
        ]
      ]
    ]
    "textoCompleto" => "<span class="elsevierStyleSections"><span id="sec0005" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0025"><a name="p164"></a>Introduction</span><p id="par0005" class="elsevierStylePara elsevierViewall">Integrated Reservoir Modeling &#40;IRM&#41; is the most accepted way to obtain the spatial distribution of petrophysical properties in oilfields &#40;<a class="elsevierStyleCrossRef" href="#bib0030">Cosentino&#44; 2001</a>&#41;&#46; An important and common task performed in this method is the estimation of permeability&#44; because it is well known that this petrophysical property is quite informative about the oil flux patterns in a reservoir&#46; However&#44; it is difficult to obtain direct information about permeability&#44; and therefore it is necessary to find models of dependence with another petrophysical property &#40;such as porosity&#44; water saturation&#44; etc&#46;&#41; in order to have an estimation of its profile &#40;<a class="elsevierStyleCrossRef" href="#bib0145">Landa <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1996</a>&#41;&#46;</p><p id="par0010" class="elsevierStylePara elsevierViewall">The linear regression approach is the most common way to model permeability values using other petrophysical properties as covariates &#40;<a class="elsevierStyleCrossRef" href="#bib0015">Balan <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1995</a>&#41;&#46; To meet the requirement of linearity it is common to perform transformations which imply that the final result could be biased when it is back transformed&#59; there are approaches that&#44; in order to fit linear models&#44; apply logarithmic transformations to induce this behavior&#44; for example&#44; to use the Cokriging method it is necessary to have a linear relationship since it requires the linear corregionalization model &#40;Sanjay and Journel 1994&#41;&#46;</p><p id="par0015" class="elsevierStylePara elsevierViewall">The main disadvantage of linear dependency models is their lack of ability to capture and model the dependence structure or pattern &#40;<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy <span class="elsevierStyleItalic">et al</span>&#46; 2005</a>&#41;&#46; In other words&#44; traditional methods cannot capture the complex variability of data&#44; in terms of variance or standard deviation&#59; hence&#44; the predicted permeability will not reproduce extreme values of the real data&#46; In other words&#44; these approaches will not be able to represent impermeable barriers or high permeability zones&#44; and from a fluid flow point of view this aspects are the most important characteristics that determine the patterns of fluid movement&#46; In this context&#44; the predicted permeability profile&#44; using linear estimators&#44; will not be an effective approximation due to its oversmoothing nature&#46;</p><p id="par0020" class="elsevierStylePara elsevierViewall">On the other hand&#44; model-free function estimators like artificial neural networks are very flexible tools that have been used to model permeability&#46; However&#44; neural networks have some disadvantages too&#46; First&#44; the training process has to be done with caution and can be a lengthy process&#46; The good results obtained by this technique are reached using a comprehensive training data set&#44; which is not always available&#46; On the other hand&#44; failing in correctly calibrating the network may result in aberrant results&#46; Another point to take into consideration is that the methodology is not yet an &#8220;off the shelf&#8221; application and requires expertise by the geoscientist &#40;<a class="elsevierStyleCrossRef" href="#bib0030">Cosentino&#44; 2001</a>&#41;&#46;</p><p id="par0025" class="elsevierStylePara elsevierViewall">Another alternative are Bayesian methods&#59; however&#44; the traditional framework of the Bayesian analysis is based on the multivariate normal distribution where the lower and upper tails are symmetrical&#46; <a class="elsevierStyleCrossRef" href="#bib0010">Armstrong <span class="elsevierStyleItalic">et al</span>&#46; &#40;2004&#41;</a> proposed an alternative Bayesian analysis that is based on Archimedean copulas where the joint distribution does not have to be normal and there is flexibility to have a lower tail or upper tail dependence based on the specific type of copula&#46;</p><p id="par0030" class="elsevierStylePara elsevierViewall">Constructing numerical models of the reservoir that honor all available data &#40;core measurements&#44; well logs&#44; seismic and geological interpretations&#44; etc&#46;&#41; having sparse knowledge of rock properties&#44; leads us to consider the stochastic simulation approach &#40;<a class="elsevierStyleCrossRef" href="#bib0060">Deutsch&#44; 1992</a>&#41;&#46; This is not a new concept &#40;<a class="elsevierStyleCrossRef" href="#bib0105">Haldorsen and Damsleth&#44; 1990</a>&#59; <a class="elsevierStyleCrossRef" href="#bib0125">Journel and Alabert&#44; 1990</a>&#41;&#44; stochastic models of physical systems are used extensively in many disciplines&#46;</p><p id="par0035" class="elsevierStylePara elsevierViewall">Stochastic simulation is the process of building alternate&#44; equally probable models of the spatial distribution of a random function&#46; It is said that a simulation is conditional if the resulting realizations honor the raw data values at their locations&#46; The most straightforward algorithms for generating realizations of a multivariate Gaussian field is provided by Sequential Gaussian Simulation &#40;SGS&#41; and Sequential Indicator Simulation &#40;SIS&#41;&#44; which are extensively used to perform permeability simulations &#40;<a class="elsevierStyleCrossRef" href="#bib0120">Holden <span class="elsevierStyleItalic">et al</span>&#46; 1995</a>&#41;&#46; However despite of their improvements &#40;<a class="elsevierStyleCrossRef" href="#bib0130">Journel and Zhu&#44; 1990</a>&#59; Suro-Perez&#44; 1990&#41; these methods are limited to cases when the spatial continuity is characterized by stationary two-point statistics and to data that is defined on the same support &#40;<a class="elsevierStyleCrossRef" href="#bib0060">Deutsch&#44; 1992</a>&#41;&#46;</p><p id="par0040" class="elsevierStylePara elsevierViewall">A competitive and more systematic method for predicting permeability may be achieved by applying stochastic joint simulations&#44; in which the correct specification of dependence pattern between petrophysical properties is crucial &#40;<a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham&#44; 1994</a>&#41;&#46; According to Deutsch this approach basically consists of an annealing geostatistical cosimulation of porosity-permeability using their empirical joint distribution&#46;<a name="p165"></a></p><p id="par0045" class="elsevierStylePara elsevierViewall">A modification of Deutsch&#8217;s methodology was presented by <a class="elsevierStyleCrossRef" href="#bib0070">D&#237;az-Viera and Casar-Gonz&#225;lez&#44; &#40;2005&#41;</a>&#46; Here&#44; it was proposed the use of a bivariate t-copula to construct the joint distribution function between porosity and permeability rather than use the joint distribution function of the sample data&#44; the dependence structure was specified by matching concordance measures such as Kendall&#8217;s &#964; and Spearman&#8217;s &#961;&#46; The above methodology was applied to simulate permeability from the porosity profile in a double porosity carbonate systems&#44; restricted to one-dimensional case at well-log scale &#40;<a class="elsevierStyleCrossRef" href="#bib0075">D&#237;az-Viera <span class="elsevierStyleItalic">et al</span>&#46; 2006</a>&#41;&#46;</p><p id="par0050" class="elsevierStylePara elsevierViewall">While D&#237;az-Viera and Casar-Gonz&#225;lez proposal can reproduce adequately the observed data&#44; and also their extreme values using a t-copula&#44; the critical problem is that the applied copula is parametric&#44; and consequently&#44; it is based on a given distribution function&#44; the student t&#46; In that sense&#44; expecting a single copula family to be able to model any kind of dependency relationship seems to be too restrictive&#44; at least for the petrophysical properties under consideration&#46; It does not mean that parametric copulas have not practical use&#46; It means that petrophysical properties modeling&#44; where dependence structure can be complex or non linear&#44; using parametric copulas could be available but in a very complicated form&#44; <a class="elsevierStyleCrossRef" href="#bib0170">Sancetta and Satchell &#40;2004&#41;</a>&#46;</p><p id="par0055" class="elsevierStylePara elsevierViewall">The copula approach has been successfully used to model dependence patterns in few areas of oil industry&#44; for example&#44; for field development decision process &#40;<a class="elsevierStyleCrossRef" href="#bib0050">Accioly and Chiyoshi 2004</a>&#41; or to model dependence in petroleum decision making &#40;<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy <span class="elsevierStyleItalic">et al</span>&#46; 2005</a>&#41;&#46;</p><p id="par0060" class="elsevierStylePara elsevierViewall">In a recent work we use the copula approach to model the bivariate dependence between petrophysical properties in a complete nonparametric fashion &#40;<a class="elsevierStyleCrossRef" href="#bib0110">Hern&#225;ndez-Maldonado <span class="elsevierStyleItalic">et al</span>&#46; 2012</a>&#41;&#44; it means we did not impose a priori a parametric joint distribution function to characterize the dependence structure of the sample&#59; instead&#44; we used a bivariate Bernstein copula function to model the intrinsic distribution of data values&#46; In this work we introduce a multivariate method for estimating empirical dependence among several petrophysical properties using a nonparametric copula&#46; The dependency model obtained is then used to stochastically simulate one property &#40;primary variable&#41; given other ones &#40;secondary variables&#41;&#46; Here&#44; we apply a method widely used in geostatistics&#44; the Simulated Annealing method&#44; which is a global optimization framework where we can add restrictions to simulation&#46; In this case&#44; in addition of using multivariate dependence by sampling the empirical copula model&#44; a spatial correlation function &#40;a variogram&#41; for the secondary variable is imposed&#46; Additionally&#44; if any raw data of the primary variable is available&#44; it can be exactly honored&#44; that is&#44; the method can be conditional&#46;</p><p id="par0065" class="elsevierStylePara elsevierViewall">In this paper this method is applied to simulate rock permeability using a trivariate copula model where permeability is described by Vugular Porosity and Shear Wave Velocity&#46; It is introduced a two stage simulation method in a multivariate fashion to model stochastically the spatial distribution of permeability at well log scale&#46; Conditional and nonconditional simulations&#44; and the median of simulations are performed in order to show the results that this method provides&#59; also it is established that&#44; as far as more descriptive variables integrate the Bernstein copula it will reproduce permeability values in a very precise form&#44; consequently it will not be necessary to perform many simulations to reduce small scale variability&#46;</p></span><span id="sec0010" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0030">Brief introduction to multivariate copulas</span><p id="par0070" class="elsevierStylePara elsevierViewall">The main disadvantage of dependency models based on a linear regression approach is their lack of ability to capture and model more general dependence structure&#46; As an alternative a copula can model the joint distribution of petrophysical properties in better ways&#46; The essence of the copula approach is that a joint distribution of random variables can be expressed as a function of the marginal distributions &#40;<a class="elsevierStyleCrossRef" href="#bib0005">Al-Harthy <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005</a>&#41;&#46;</p><p id="par0075" class="elsevierStylePara elsevierViewall">According to Sklar&#8217;s theorem &#40;<a class="elsevierStyleCrossRef" href="#bib0190">Sklar&#44; 1959</a>&#41;&#44; the underlying copula associated to a multivariate random vector &#40;<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">1</span>&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">2</span>&#44; &#46;&#46;&#46;&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span>&#41; represents a functional link between the joint probability distribution and the univariate marginal distributions <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf">1</span>&#44; <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf">2</span>&#44; &#46;&#46;&#46;&#44; <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">M</span></span> respectively&#58;</p><p id="par0080" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0005"></elsevierMultimedia></p><p id="par0085" class="elsevierStylePara elsevierViewall">For all <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">1</span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">2</span>&#44; &#46;&#46;&#46;&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span> in the extended real numbers system&#44; where <span class="elsevierStyleItalic">C&#58;</span>&#91;0&#44; 1&#93;<span class="elsevierStyleSup"><span class="elsevierStyleItalic">m</span></span> &#8594;&#91;0&#44; 1&#93; the underlying copula is unique whenever <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span><span class="elsevierStyleItalic">&#44; X</span><span class="elsevierStyleInf">2</span>&#44; and <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span> are continuous random variables&#46; Therefore&#44; all the information about the dependence between continuous random variables is contained in their corresponding copula&#46; Several properties may be derived for copulas &#40;<a class="elsevierStyleCrossRef" href="#bib0155">Nelsen&#44; 2006</a>&#41;&#44; and among them <a name="p166"></a>we have an immediate corollary from Sklar&#8217;s theorem&#58; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span>&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">2</span></span> and <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span> are independent continuous random variables if and only if their underlying copula is <span class="elsevierStyleItalic">C&#40;u</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span><span class="elsevierStyleItalic">&#44;&#46;&#46;&#46;&#44;u</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">n</span></span><span class="elsevierStyleItalic">&#41;</span>&#61;<span class="elsevierStyleItalic">u&#44;&#46;&#46;&#46;&#44;u</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">n</span></span>&#46;</p><p id="par0090" class="elsevierStylePara elsevierViewall">Let S&#61;&#123;&#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">11</span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">21</span>&#44; &#46;&#46;&#46; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span>1</span>&#41;&#44; &#46;&#46;&#46;&#44; &#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">1<span class="elsevierStyleItalic">n</span></span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">2<span class="elsevierStyleItalic">n</span></span>&#44; &#46;&#46;&#46;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">mn</span></span>&#41;&#125; be <span class="elsevierStyleItalic">n</span> observations of a random vector &#40;<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">1</span>&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">2</span>&#44; &#46;&#46;&#46;&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span>&#41;&#46; We may obtain empirical estimates for the marginal distributions <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">1</span>&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">2</span>&#44; &#46;&#46;&#46;&#44; <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span> by means of&#58;</p><p id="par0095" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0010"></elsevierMultimedia></p><p id="par0100" class="elsevierStylePara elsevierViewall">where stands for an indicator function which takes value whenever its argument is true&#44; and otherwise&#46; It is well-known &#40;<a class="elsevierStyleCrossRef" href="#bib0025">Billingsley&#44; 1995</a>&#41; that the empirical distribution F&#710;j is a consistent estimator of <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">j</span></span> that is&#44; <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">j</span></span> that is F&#710;j&#40;t&#41;<span class="elsevierStyleItalic">i</span>s a converges almost surely to <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">j</span></span>&#40;<span class="elsevierStyleItalic">t</span>&#41; <span class="elsevierStyleItalic">as n &#8594;&#8734;</span> all <span class="elsevierStyleItalic">t&#46;</span></p><p id="par0105" class="elsevierStylePara elsevierViewall">Similarly&#44; we have the empirical copula &#40;<a class="elsevierStyleCrossRef" href="#bib0055">Deheuvels&#44; 1979</a>&#41;&#44; a function <span class="elsevierStyleItalic">C</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">n</span></span> with domain &#123;1n&#59;i&#61;0&#44;1&#44;&#8230;n&#125;n &#40;where n is the sample size&#41; defined as</p><p id="par0110" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0015"></elsevierMultimedia></p><p id="par0115" class="elsevierStylePara elsevierViewall">and its convergence to the true copula has also been proved &#40;<a class="elsevierStyleCrossRef" href="#bib0100">Fermanian <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2004</a>&#41;&#46; The empirical copula is not a copula&#44; since it is only defined on a finite grid&#44; not in the whole unit hypercube &#91;0&#44; 1&#93;<span class="elsevierStyleSup"><span class="elsevierStyleItalic">n</span></span>&#44; but by Sklar&#8217;s Theorem &#40;<a class="elsevierStyleCrossRef" href="#bib0190">Sklar&#44; 1959</a>&#41; may be extended to a copula&#46;<span class="elsevierStyleItalic">C</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">n</span></span> may be extended to a copula&#46;</p><p id="par0120" class="elsevierStylePara elsevierViewall">Sklar&#8217;s theorem is completely general and a joint distribution function can be constructed using the copula function&#46; The copula separates the marginal distributions from correlation and the copula itself can capture the dependence structure&#46; This is an essential property of copulas&#46;</p><p id="par0125" class="elsevierStylePara elsevierViewall">A general example of a multivariate copula model that describes the relationship of three random variables &#40;a trivariate copula&#41; according to Sklar&#8217;s theorem is shown&#58;</p><p id="par0130" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0020"></elsevierMultimedia></p><p id="par0135" class="elsevierStylePara elsevierViewall">There are different classes and families of copulas but&#44; for this study&#44; we will work with the Bernstein copula&#44; which is a nonparametric copula since it is a nonparametric smoothing based on the empirical copula&#46;</p></span><span id="sec0015" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0035">Joint stochastic simulation method using a multivariate Bernstein copula</span><p id="par0140" class="elsevierStylePara elsevierViewall">The method presented here basically consists in generating a stochastic simulation of a primary variable&#44; where the joint distribution function of its covariates is modeled using a nonparametric copula&#46; In other words&#44; the spatial dependence and dependence pattern between variables are decoupled&#59; in this context&#44; the Bernstein copula only models the dependence between variables and the spatial dependence is modeled by a variogram applying the Simulated Annealing method&#46;</p><p id="par0145" class="elsevierStylePara elsevierViewall">Basically this proposal is a two-stage modeling algorithm method&#59; first&#44; the pattern relationship between petrophysical properties under study&#59; and second&#44; the spatial structure&#46; In the first stage it is modeled the dependence structure of petrophysical properties using a multivariate nonparametric copula &#40;the Bernstein copula&#41;&#44; then a geostatistical simulation of primary variable is performed using simulated annealing technique&#44; whose objective function is the variogram model &#40;<a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham&#44; 1994</a>&#41;&#46; A detailed description of each step of the algorithm will be described&#46;</p><span id="sec0020" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0040">Multivariate Bernstein Copula Modeling</span><p id="par0150" class="elsevierStylePara elsevierViewall">Each petrophysical property is modeled as an absolutely continuous random variable <span class="elsevierStyleItalic">X</span> with unknown marginal distribution function <span class="elsevierStyleItalic">F&#46;</span> For simulation of continuous random variables&#44; the use of the empirical distribution function <a class="elsevierStyleCrossRef" href="#eq0010">&#40;2&#41;</a> is not appropriate since F&#710;j is a step function&#44; and therefore discontinuous&#44; so a smoothing technique is needed&#46; Since our main goal is to simulate a primary variable using more than one descriptive variable&#44; it will be better to have a smooth estimation of the marginal quantile function <span class="elsevierStyleItalic">Q</span>&#40;<span class="elsevierStyleItalic">u</span>&#41;<span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">F</span><span class="elsevierStyleSup">-1</span>&#40;<span class="elsevierStyleItalic">u</span>&#41;<span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">inf</span>&#123;<span class="elsevierStyleItalic">x</span>&#58;<span class="elsevierStyleItalic">F</span>&#40;<span class="elsevierStyleItalic">x</span>&#41;&#8805;<span class="elsevierStyleItalic">u</span>&#125;&#44;0&#8804;<span class="elsevierStyleItalic">u</span>&#8804;1 which is possible by means of Bernstein polynomials as in <a class="elsevierStyleCrossRef" href="#bib0150">Munoz-Perez and Fernandez-Palac&#237;n &#40;1987&#41;</a>&#46;</p><p id="par0155" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0025"></elsevierMultimedia></p><p id="par0160" class="elsevierStylePara elsevierViewall">For a smooth estimation of the underlying copula we make use of the Bernstein copula &#40;<a class="elsevierStyleCrossRef" href="#bib0170">Sancetta and Satchell&#44; 2004</a>&#59; <a class="elsevierStyleCrossRef" href="#bib0175">Sancetta&#44; 2007</a>&#41;&#58;</p><p id="par0165" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0030"></elsevierMultimedia></p><p id="par0170" class="elsevierStylePara elsevierViewall"><a name="p167"></a>where&#59;</p><p id="par0175" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0035"></elsevierMultimedia></p><p id="par0180" class="elsevierStylePara elsevierViewall">For every &#40;<span class="elsevierStyleItalic">u</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span>&#44; &#8230;&#44; <span class="elsevierStyleItalic">u</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">m</span></span>&#41; in the unit hypercube &#91;0&#44; 1&#93;<span class="elsevierStyleSup"><span class="elsevierStyleItalic">m</span></span>&#44; and &#945;v1n1&#44;&#8230;&#44;vmnm is the empirical copula&#44; defined in <a class="elsevierStyleCrossRef" href="#eq0015">&#40;3&#41;</a>&#46;</p><p id="par0185" class="elsevierStylePara elsevierViewall">For example&#44; the trivariate Bernstein copula model derived from <a class="elsevierStyleCrossRef" href="#eq0030">&#40;6&#41;</a> is&#58;</p><p id="par0190" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0040"></elsevierMultimedia></p></span><span id="sec0025" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0045">Sampling algorithm of a Bernstein copula model</span><p id="par0195" class="elsevierStylePara elsevierViewall">For a pair of random variables &#40;<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span><span class="elsevierStyleItalic">&#44; X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">2</span></span>&#41; with joint distribution function <span class="elsevierStyleItalic">H</span> and underlying copula <span class="elsevierStyleItalic">C</span> we need to generate an observation of uniform &#40;0&#44; 1&#41; random variables &#40;<span class="elsevierStyleItalic">U&#44; V</span>&#41; whose joint distribution function is <span class="elsevierStyleItalic">C</span> and then transform those uniform variables as in step 3 of the sampling bivariate algorithm&#46; For generating such pair <span class="elsevierStyleItalic">&#40;u&#44; v&#41;</span> it is used a conditional distribution method&#44; this method needs the conditional distribution function for <span class="elsevierStyleItalic">V</span> given <span class="elsevierStyleItalic">U</span><span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">u</span>&#44; which we denote as <span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">u</span></span> &#40;<span class="elsevierStyleItalic">v</span>&#41;&#58;</p><p id="par0200" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0045"></elsevierMultimedia></p><p id="par0205" class="elsevierStylePara elsevierViewall">where C&#710;B is the bivariate Bernstein copula model&#44; obtained by <a class="elsevierStyleCrossRef" href="#eq0030">&#40;6&#41;</a>&#46;</p><p id="par0210" class="elsevierStylePara elsevierViewall">To simulate replications from the random vector with the dependence structure estimated from the observed data S&#58;&#61;&#123;&#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">11</span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">21</span>&#41;&#44; &#46;&#46;&#46;&#44; &#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">1<span class="elsevierStyleItalic">n</span></span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">2<span class="elsevierStyleItalic">n</span></span>&#41;&#125; it is applied the following algorithm&#58;</p></span></span><span id="sec1035" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect1055">Sampling bivariate algorithm</span><p id="par3220" class="elsevierStylePara elsevierViewall"><ul class="elsevierStyleList" id="lis0005"><li class="elsevierStyleListItem" id="lsti0005"><span class="elsevierStyleLabel">1</span><p id="par0220" class="elsevierStylePara elsevierViewall">Generate two independent and continuous Uniform &#40;0&#44;1&#41; random variates<span class="elsevierStyleItalic">u</span> and<span class="elsevierStyleItalic">t</span></p></li><li class="elsevierStyleListItem" id="lsti0010"><span class="elsevierStyleLabel">2</span><p id="par0225" class="elsevierStylePara elsevierViewall">Set<span class="elsevierStyleItalic">v</span><span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">u</span></span><span class="elsevierStyleSup">-1</span> &#40;<span class="elsevierStyleItalic">t</span>&#41;&#59; where<span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">u</span></span> is defined in <a class="elsevierStyleCrossRef" href="#eq0045">&#40;9&#41;</a>&#46;</p></li><li class="elsevierStyleListItem" id="lsti0015"><span class="elsevierStyleLabel">3</span><p id="par0230" class="elsevierStylePara elsevierViewall">The desired pair is &#40;x1&#44;x2&#41;&#61;&#40;Q&#732;&#59;n&#40;u&#41;&#44;R&#732;&#59;n&#40;v&#41;&#41; where Q&#732;&#59;n and R&#732;&#59;n according to <a class="elsevierStyleCrossRef" href="#eq0025">&#40;5&#41;</a>&#44; are the estimated and smoothed quantile functions<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">1</span> of<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">2</span> and&#44; respectively&#46;</p></li></ul></p><p id="par0235" class="elsevierStylePara elsevierViewall">For the multivariate case we must solve equations that represent conditional distribution functions for <span class="elsevierStyleItalic">W</span> given <span class="elsevierStyleItalic">U&#61;u&#44;V&#61;v</span>&#46;</p><p id="par0240" class="elsevierStylePara elsevierViewall">To simulate replications from the random vector with the dependence structure estimated from the observed data S&#58;&#61;&#123;&#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">11</span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">21</span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">31</span>&#41;&#44; &#46;&#46;&#46;&#44; &#40;<span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">1<span class="elsevierStyleItalic">n</span></span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">2<span class="elsevierStyleItalic">n</span></span>&#44; <span class="elsevierStyleItalic">x</span><span class="elsevierStyleInf">3<span class="elsevierStyleItalic">n</span></span>&#41;&#125; it is applied the next algorithm&#46;</p></span><span id="sec2040" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect1060">Sampling trivariate algorithm</span><p id="par1250" class="elsevierStylePara elsevierViewall"><ul class="elsevierStyleList" id="lis0010"><li class="elsevierStyleListItem" id="lsti0020"><span class="elsevierStyleLabel">1</span><p id="par0250" class="elsevierStylePara elsevierViewall">Generate three independent and continuous Uniform &#40;0&#44; 1&#41; random variables<span class="elsevierStyleItalic">u</span> and<span class="elsevierStyleItalic">t</span><span class="elsevierStyleInf">1</span><span class="elsevierStyleItalic">&#44; t</span><span class="elsevierStyleInf">2</span></p></li><li class="elsevierStyleListItem" id="lsti0025"><span class="elsevierStyleLabel">2</span><p id="par0255" class="elsevierStylePara elsevierViewall">Set<span class="elsevierStyleItalic">v</span><span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">u</span></span><span class="elsevierStyleSup">-1</span> &#40;<span class="elsevierStyleItalic">t</span><span class="elsevierStyleInf">1</span>&#41; where<span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">u</span></span> is defined in <a class="elsevierStyleCrossRef" href="#eq0045">&#40;9&#41;</a>&#46;</p></li><li class="elsevierStyleListItem" id="lsti0030"><span class="elsevierStyleLabel">3</span><p id="par0260" class="elsevierStylePara elsevierViewall">Set<span class="elsevierStyleItalic">w</span><span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span><span class="elsevierStyleItalic">c</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">uv</span></span><span class="elsevierStyleSup">-1</span> &#40;<span class="elsevierStyleItalic">t</span><span class="elsevierStyleInf">2</span>&#41; where <elsevierMultimedia ident="eq0050"></elsevierMultimedia></p></li></ul></p><p id="par0270" class="elsevierStylePara elsevierViewall">where is the trivariate Bernstein copula model <a class="elsevierStyleCrossRef" href="#eq0040">&#40;8&#41;</a>&#46; <ul class="elsevierStyleList" id="l1035"><li class="elsevierStyleListItem" id="lsti0035"><span class="elsevierStyleLabel">4</span><p id="par0275" class="elsevierStylePara elsevierViewall">The desired vector is &#40;x1&#44;x2&#44;x3&#41;&#61;&#40;Q&#732;&#59;n&#40;u&#41;&#44;R&#732;&#59;n&#40;v&#41;&#44;H&#732;&#59;n&#40;w&#41;&#41; where Q&#732;&#59;n&#40;u&#41;&#44;R&#732;&#59;n&#40;v&#41; and&#44; H&#732;&#59;n&#40;w&#41; according to <a class="elsevierStyleCrossRef" href="#eq0025">&#40;5&#41;</a>&#44; are the estimated and smoothed quantile functions of<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">1</span></span>&#44;<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">2</span></span>&#44; y<span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf"><span class="elsevierStyleItalic">3</span></span>&#44; respectively&#46;</p></li></ul></p><span id="sec0030" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0050">Simulated Annealing method</span><p id="par0280" class="elsevierStylePara elsevierViewall">To perform stochastic simulations of a primary variable it is applied a method widely used in geostatistics&#44; the Simulated Annealing method&#44; which is a global optimization framework where we can add restrictions to simulation&#46; In this case&#44; in addition of using multivariate dependence by sampling the copula model&#44; a spatial correlation function &#40;a variogram&#41; for the secondary variable is imposed&#46; Additionally&#44; if any raw data of the primary variable is available&#44; it can be exactly honored&#44; that is&#44; the method can be conditional&#46;</p><p id="par0285" class="elsevierStylePara elsevierViewall">A more detailed explanation of each single step of simulated annealing method&#44; as well as technical terminology&#44; can be found in <a class="elsevierStyleCrossRef" href="#bib0080">Dreo <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#44; Dafflon et&#46; al&#46;&#44; &#40;<a class="elsevierStyleCrossRef" href="#bib0035">2009a</a>&#44; <a class="elsevierStyleCrossRef" href="#bib0040">2009b</a>&#41;&#44; <a class="elsevierStyleCrossRef" href="#bib0045">Dafflon and Barrash &#40;2012&#41;</a>&#44; <a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado <span class="elsevierStyleItalic">et al</span>&#46; &#40;2012&#41;</a>&#44; Deutsch &#40;2002&#41;&#46; In this work&#44; we will just briefly mention its steps&#46;</p><p id="par0290" class="elsevierStylePara elsevierViewall">a&#46; Generating an initial configuration&#44; known as &#8220;seed&#8221;&#46; An initial configuration is the starting point of the simulation&#44; also&#44; it can be <a name="p168"></a>considered as a &#8220;possible solution&#8221; of the spatial disposition of the primary variable&#46; There are some ways of proposing this configuration&#44; see Deutsch &#40;2002&#41;&#44; <a class="elsevierStyleCrossRef" href="#bib0080">Dr&#233;o et al&#46; &#40;2006&#41;</a>&#46; One of most popular of them is just proposing random numbers with uniform distribution&#46; We decided to use the multivariate Bernstein copula sampling algorithm &#40;section 3&#46;1&#46;1&#41; to obtain this starting point in order to have a better approximation of the final solution&#46;</p><p id="par0295" class="elsevierStylePara elsevierViewall">b&#46; Defining the objective function&#46; Since we use the Bernstein copula we do not have to include univariate histograms into the objective function because they are reproduced automatically&#59; in the same way&#44; as Bernstein copula models the dependence structure between variables&#44; it is not necessary to include neither a correlation coefficient&#46; In other words&#44; in order to satisfy the spatial distribution of the primary variable&#44; we are able to propose an objective function that consists only in one term <a class="elsevierStyleCrossRef" href="#eq0055">&#40;11&#41;</a> the variogram function</p><p id="par0300" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0055"></elsevierMultimedia></p><p id="par0305" class="elsevierStylePara elsevierViewall">Where&#58;</p><p id="par0310" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0060"></elsevierMultimedia></p><p id="par0315" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0065"></elsevierMultimedia></p><p id="par0320" class="elsevierStylePara elsevierViewall">c&#46; Measuring the initial configuration energy&#46; By <a class="elsevierStyleCrossRef" href="#eq0060">&#40;12&#41;</a>&#44; it is calculated the variogram of the initial configuration and it is compared to the variogram of the real configuration <a class="elsevierStyleCrossRef" href="#eq0055">&#40;11&#41;</a>&#44; this value is the energy of the starting point&#44; and the main idea is that it must be decreased during the simulation procedure&#44; see Deutsch &#40;2002&#41;&#46; This step gives us an idea of the seed&#180;s quality&#44; if we use random numbers to generate it&#44; its measured energy will be very high&#44; however&#44; as long as this method uses a sample given by the Bernstein copula the initial energy is significantly small&#46;</p><p id="par0325" class="elsevierStylePara elsevierViewall">d&#46; Obtaining the initial temperature and the annealing schedule&#46; The initial temperature is obtained by <a class="elsevierStyleCrossRef" href="#eq0065">&#40;13&#41;</a>&#58;</p><p id="par0330" class="elsevierStylePara elsevierViewall">The annealing schedule is obtained following the procedure proposed by <a class="elsevierStyleCrossRef" href="#bib0080">Dr&#233;o et al&#46; &#40;2006&#41;</a>&#46;</p><p id="par0335" class="elsevierStylePara elsevierViewall">e&#46; Ending up a realization&#46; The simulation ends whenever the objective function error is reached&#44; or an accumulation of 3 stages without change occurs&#44; or when the maximum number of attempted perturbations is reached&#46;</p><p id="par0340" class="elsevierStylePara elsevierViewall">It is shown a diagram of the method to perform multivariate stochastic simulations in&#44; <a class="elsevierStyleCrossRef" href="#f0005">Figure1</a>&#46;</p><elsevierMultimedia ident="f0005"></elsevierMultimedia></span></span><span id="sec0035" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0055">Case study</span><p id="par0345" class="elsevierStylePara elsevierViewall">It will be modeled the Permeability of double porosity carbonate formations of a South Florida Aquifer in the western Hillsboro Basin of Palm Beach County&#44; Florida&#46; Based in the algorithm described in Section 3&#44; we propose a copula-based approach to model the relationship between permeability&#44; porosity and Shear Wave Velocity &#40;S-waves&#41;&#46; The dependency model obtained will be then used to stochastically simulate permeability using Porosity and Shear Wave velocity&#46;</p><span id="sec0040" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0060">Data description</span><p id="par0350" class="elsevierStylePara elsevierViewall">The characterization of this aquifer for the borehole and field scales is given in <a class="elsevierStyleCrossRef" href="#bib0160">Parra <span class="elsevierStyleItalic">et al</span>&#46; &#40;2001&#41;</a>&#44; <a class="elsevierStyleCrossRef" href="#bib0165">Parra and Hackert &#40;2002&#41;</a>&#44; and a hydrogeological situation is described by <a class="elsevierStyleCrossRef" href="#bib0020">Bennett <span class="elsevierStyleItalic">et al</span>&#46; &#40;2002&#41;</a>&#46; The interpretation <a name="p169"></a>of the borehole data&#44; the determination of the matrix and secondary porosity and secondary-pore types &#40;shapes of spheroids approximating secondary pores&#41; were presented by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#44; where to determine the pore microstructure of aquifer carbonate formations the authors applied the petrophysical inversion technique that consists in minimizing a cost function that includes the sum of weighted square differences between the experimentally measured and theoretically calculated logs as in <a class="elsevierStyleCrossRef" href="#bib0135">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2004&#41;</a>&#46;</p><p id="par0355" class="elsevierStylePara elsevierViewall">In this case the following well logs were used for joint simultaneous inversion as input data&#58; resistivity log&#44; transit times of the Pand S-waves &#40;acoustic log&#41;&#44; total porosity &#40;neutron log&#41;&#44; and formation density &#40;density log&#41;&#46; To calculate the theoretical acoustic and resistivity logs the double-porosity model for describing carbonate formations was applied&#58; <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#46;</p><p id="par0360" class="elsevierStylePara elsevierViewall">This model treats carbonate rocks as a composite material that consists of a homogeneous isotropic matrix &#40;solid skeleton and matrix pore system&#41; where the secondary pores of different shapes are embedded&#46; The secondary pores were approximated by spheroids with variable aspect ratios to represent different secondary porosity types&#58; vugs &#40;close-to-sphere shapes&#41;&#44; quasi-vugs &#40;oblate spheroids&#41;&#44; channels &#40;prolate spheroids&#41;&#44; and microfractures &#40;flattened spheroids&#41; <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#46;</p><p id="par0365" class="elsevierStylePara elsevierViewall">We used the results of inversion obtained by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a> for carbonate formations of South Florida Aquifer that includes the following petrophysical characteristics&#58; matrix porosity&#44; secondary vugular and crack porosities &#40;<a class="elsevierStyleCrossRef" href="#f0010">Figure 2</a>&#41;&#46;</p><elsevierMultimedia ident="f0010"></elsevierMultimedia><p id="par0370" class="elsevierStylePara elsevierViewall">The secondary-porosity system of this formation has complex microstructure and corresponds to a model with two types of pore shapes&#58; cracks &#40;flattened ellipsoids&#41; with the overall porosity of 2&#37; and vugs &#40;close to sphere&#41; with the porosity variations in the range of 10-30&#37;&#46;</p></span><span id="sec0045" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0065">Statistical data analysis</span><p id="par0375" class="elsevierStylePara elsevierViewall">Given multivariate data&#44; it is common to start choosing as explanatory variables those who exhibit higher dependence with the variable that is to be explained&#46; We measured the dependence between the petrophysical properties in terms of the dependence index &#934; proposed by <a class="elsevierStyleCrossRef" href="#bib0115">Hoeffding &#40;1941&#41;</a>&#44; which satisfies all desirable properties for a dependence measure for continuous random variables&#44; see <a class="elsevierStyleCrossRef" href="#bib0155">Nelsen &#40;2006&#41;</a>&#46;</p><p id="par0380" class="elsevierStylePara elsevierViewall">From all possible explanatory variables &#40;PHITOT&#44; PHIMAT&#44; PHICR&#44; PHIV&#44; VS Meas&#44; VP meas&#41; for permeability &#40;K&#41;&#44; for the first explanatory random variable it was chosen relative vugular porosity &#40;PHIV&#41; since it exhibited the highest dependence <span class="elsevierStyleItalic">&#934;</span> &#40;<span class="elsevierStyleItalic">PHI&#44; K</span>&#41;&#61;0&#46;71 &#40;on a &#91;0&#44;1&#93; scale&#41;&#46; In choosing a second explanatory random variable we need&#44; in addition&#44; to have a high dependence with permeability&#44; and the lowest possible dependence with the first explanatory variable &#40;PHIV&#41;&#44; otherwise it would mean that it is quite similar to it and it will add no significant information to what the first one already can provide&#46; Under this criteria&#44; the second best choice was Share Wave Velocity &#40;VS meas&#41;&#44; <a name="p170"></a>with <span class="elsevierStyleItalic">&#934;</span> &#40;<span class="elsevierStyleItalic">VS meas&#44; K</span>&#41;<span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span>0&#46;60 and <span class="elsevierStyleItalic">&#934;</span> &#40;<span class="elsevierStyleItalic">PHIV&#44; VS meas</span>&#41;<span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span>0&#46;55</p><p id="par0385" class="elsevierStylePara elsevierViewall">&#93;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a> shows all dependence index <span class="elsevierStyleItalic">&#934;</span>of explanatory variables that are the results of inversion obtained by <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#59; Total porosity &#40;PHITOT&#41;&#44; Matrix porosity &#40;PHIMAT&#41;&#44; Crack porosity &#40;PHICR&#41;&#44; Vugular porosity &#40;PHIV&#41;&#44; Share Wave Velocity &#40;VS Meas&#41;&#44; P-Wave velocity &#40;VP meas&#41;&#46;</p><elsevierMultimedia ident="tbl0005"></elsevierMultimedia><p id="par0390" class="elsevierStylePara elsevierViewall">In <a class="elsevierStyleCrossRef" href="#tbl0005">table 1</a> explanatory variables are shown&#44; using this information we confirm that PHIV has the highest dependence&#44; but following the order of this table we should use as second explanatory variable PHITOT&#46; We did not make this choice because&#44; as we said before&#44; a second explanatory variable needs to have a high dependence with permeability&#44; and the lowest possible dependence with the first explanatory variable &#40;PHIV&#41;&#46; To see this behavior in table two we show the explanatory variables matrix dependence Hoeffding index&#46;</p><p id="par0395" class="elsevierStylePara elsevierViewall">We have chosen PHIV as the first explanatory variable&#59; it means we must pay attention to the first line of <a class="elsevierStyleCrossRef" href="#tbl0010">table 2</a>&#44; and select the variable with the lowest dependence to PHIV&#46; It should be PHICR but in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a> we see that this property offers poor information about K&#46; Finally&#44; we may choose VP meas but&#44; VS Meas explains K in some better way &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41;&#46;</p><elsevierMultimedia ident="tbl0010"></elsevierMultimedia><p id="par0400" class="elsevierStylePara elsevierViewall">In <a class="elsevierStyleCrossRef" href="#f0015">Figure 3</a> it is shown a scatter-plot distribution and histograms of PHIV&#8211;K sample Explanatory data taken from <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#44; it shows their bivariate distribution&#46; In <a class="elsevierStyleCrossRef" href="#f0020">Figure 4</a> it is shown a scatter-plot distribution and histograms of VS meas&#8211;K data values</p><elsevierMultimedia ident="f0015"></elsevierMultimedia><elsevierMultimedia ident="f0020"></elsevierMultimedia></span><span id="sec0050" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0070">Bernstein Copula Model</span><p id="par0405" class="elsevierStylePara elsevierViewall">It is modeled the trivariate Bernstein Copula using relative vugular porosities &#40;PHIV&#41;&#44; Shear Wave Velocity &#40;VS meas<span class="elsevierStyleHsp" style=""></span>&#61;<span class="elsevierStyleHsp" style=""></span>velocity of S-waves measured&#41; and K as absolutely continuous random variables <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">1</span><span class="elsevierStyleItalic">&#44; X</span><span class="elsevierStyleInf">2</span> and <span class="elsevierStyleItalic">X</span><span class="elsevierStyleInf">3</span>&#44; with unknown theoretical univariate distribution functions <span class="elsevierStyleItalic">F</span><span class="elsevierStyleInf">1</span><span class="elsevierStyleItalic">&#44; F</span><span class="elsevierStyleInf">2</span><span class="elsevierStyleItalic">and F</span><span class="elsevierStyleInf">3</span><span class="elsevierStyleItalic">&#44; see</span><a class="elsevierStyleCrossRef" href="#eq0010"><span class="elsevierStyleItalic">&#40;2&#41;</span></a><span class="elsevierStyleItalic">&#46; Observations of this random vector were obtained from</span><a class="elsevierStyleCrossRef" href="#bib0140"><span class="elsevierStyleItalic">Kazatchenko et al&#46; &#40;2006&#41;</span></a>&#46;</p><p id="par0410" class="elsevierStylePara elsevierViewall">We use <a class="elsevierStyleCrossRef" href="#eq0010">&#40;2&#41;</a> to estimate the empirical distribution function for each variable&#59; however&#59; <a name="p171"></a>in order to model them as continuous random variables they should be smoothed using <a class="elsevierStyleCrossRef" href="#eq0025">&#40;5&#41;</a>&#44; since <a class="elsevierStyleCrossRef" href="#eq0010">&#40;2&#41;</a> is not appropriate because F&#710;j is a step function and therefore discontinuous&#46; The trivariate Bernstein copula model <a class="elsevierStyleCrossRef" href="#eq0040">&#40;8&#41;</a> will be used to sample data to stochastically simulate permeability&#46;</p></span><span id="sec0055" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0075">Sampling from Bernstein copula</span><p id="par0415" class="elsevierStylePara elsevierViewall">In <a class="elsevierStyleCrossRef" href="#f0025">Figure 5</a> &#40;left side&#41; it is shown a scatterplot distribution and histograms of porosity&#8211; permeability real data &#40;sample size of&#41; taken from <a class="elsevierStyleCrossRef" href="#bib0140">Kazatchenko <span class="elsevierStyleItalic">et al</span>&#46; &#40;2006&#41;</a>&#46; At the center of the same figure&#44; it is shown a trivariate simulation of permeability values using vugular porosity and shear wave velocity as explicative variables&#46; If we do not label these two scatter plots it would be difficult to know which one of them is the original one&#46; In this case the multivariate Bernstein copula provides an acceptable the dependence structure for continuous variables with complex relationship&#46; The right side of <a class="elsevierStyleCrossRef" href="#f0025">Figure 5</a> shows 3800 simulated values of permeability&#44; according to some analysis&#44; it shows more clearly percolation patterns&#46;</p><elsevierMultimedia ident="f0025"></elsevierMultimedia></span><span id="sec0060" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0080">Permeability simulations and discussion</span><p id="par0420" class="elsevierStylePara elsevierViewall">In this section spatial simulation will be performed using simulated annealing technique&#46; Since multivariate dependence is modeled by the Bernstein copula&#44; a spatial correlation function &#40;the variogram&#41; will be used as an objective function in the annealing algorithm&#46;</p><p id="par0425" class="elsevierStylePara elsevierViewall">We propose a copula-based nonparametric approach and simulated annealing technique to model the spatial distribution of permeability using as secondary variables&#44; porosity and shear wave velocity &#40;S-waves&#41;&#46; In <a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a> we can see a single nonconditional spatial simulation of permeability &#40;K&#41; in terms of depth&#46; As we note the simulation &#40;blue line&#41; follows very close the real data behavior &#40;red line&#41;&#44; this is because we introduce more descriptive information into the Bernstein copula model&#46; However&#44; there is still a small-scale variability which does not exist in the original data values&#46;</p><elsevierMultimedia ident="f0030"></elsevierMultimedia><p id="par0430" class="elsevierStylePara elsevierViewall">As a descriptive measure of the goodness of fit of predicted values of log-permeability &#40;Log K&#41; given values of the explanatory variables PHIV and VS meas&#44; we calculated the Mean Squared Error &#40;MSE&#41; in each case&#58;</p><p id="par0435" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0070"></elsevierMultimedia></p><p id="par0440" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0075"></elsevierMultimedia></p><p id="par0445" class="elsevierStylePara elsevierViewall"><elsevierMultimedia ident="eq0080"></elsevierMultimedia></p><p id="par0450" class="elsevierStylePara elsevierViewall">where a lower value of MSE represents a better goodness of fit&#46;</p><p id="par0455" class="elsevierStylePara elsevierViewall">The simulated values of log K&#44; given PHIV and VS meas &#40;log <span class="elsevierStyleItalic">K &#124; PHIV&#44; VS meas&#41;&#44;</span> show a MSE value of 0&#46;5&#44; which is very low compared to that obtained by <a class="elsevierStyleCrossRef" href="#bib0065">Deutsch and Cockerham&#8217;s method &#40;1994&#41;</a>&#44; which MSE value is 7&#46;76&#59; also it is lower than the obtained one by D&#237;az-Viera and Casar-Gonz&#225;lez method&#44; where the MSE value is 5&#46;5&#44; details are shown in <a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado <span class="elsevierStyleItalic">et al</span>&#46;&#44; &#40;2012&#41;</a>&#46; Finally&#44; the MSE value for the trivariate copula simulation is lower than the MSE value of the bivariate copula simulation &#40;log <span class="elsevierStyleItalic">K &#124; PHIV&#41;</span><a class="elsevierStyleCrossRef" href="#eq0070">&#40;14&#41;</a>&#44; consequently it means that trivariate simulation has the best fitted values&#46;<a name="p172"></a></p><p id="par0460" class="elsevierStylePara elsevierViewall">In terms of histograms&#44; scatterplots and variograms&#44; the proposed method shows very good agreement with the original data values&#46; As we are using three variables&#44; the scatterplots of permeability can be seen in terms of K-PHIV &#40;<a class="elsevierStyleCrossRef" href="#f0035">Figure 7</a>&#41; and K-S-waves &#40;<a class="elsevierStyleCrossRef" href="#f0040">Figure 8</a>&#41;&#46; The scatterplots of the original data values are shown in <a class="elsevierStyleCrossRef" href="#f0015">Figure 3</a> and <a class="elsevierStyleCrossRef" href="#f0020">Figure 4</a> respectively&#46;</p><elsevierMultimedia ident="f0035"></elsevierMultimedia><elsevierMultimedia ident="f0040"></elsevierMultimedia><p id="par0465" class="elsevierStylePara elsevierViewall">The variogram of the simulation &#40;<a class="elsevierStyleCrossRef" href="#f0045">Figure 9</a>&#41; and its fitted model &#40;<a class="elsevierStyleCrossRef" href="#tbl0015">Table 3</a>&#41;&#44; also shows good agreement with respect to real data &#40;<a class="elsevierStyleCrossRef" href="#f0050">Figure 10</a>&#41;&#46; We can note that both variograms are almost equal&#44; but in a detailed view&#44; the empirical variogram of the simulated values is a little softer than the original one&#46;</p><elsevierMultimedia ident="f0045"></elsevierMultimedia><elsevierMultimedia ident="tbl0015"></elsevierMultimedia><elsevierMultimedia ident="f0050"></elsevierMultimedia><p id="par0470" class="elsevierStylePara elsevierViewall">As we could anticipate&#44; the dependence structure between these three petrophysical properties is well represented&#46; The extreme values&#44; and in general&#44; all of most important statistics are quite well reproduced&#44; &#40;<a class="elsevierStyleCrossRef" href="#f0055">Figure 11</a>&#41;&#46; In order to reduce small-scale variability seen in <a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a>&#44; a median of 10 nonconditional SA simulations was performed&#44; <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a>&#46; Once again&#44; the spatial structure&#44; the histograms and the Scatter-plot reproduction shows a very good agreement respecting to real data&#44; see also <a class="elsevierStyleCrossRef" href="#f0065">Figure 13</a> and <a class="elsevierStyleCrossRef" href="#f0070">Figure 14</a>&#46;</p><elsevierMultimedia ident="f0055"></elsevierMultimedia><elsevierMultimedia ident="f0060"></elsevierMultimedia><elsevierMultimedia ident="f0065"></elsevierMultimedia><elsevierMultimedia ident="f0070"></elsevierMultimedia><p id="par0475" class="elsevierStylePara elsevierViewall">In <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a> we can see that for small values of K the simulation is not following real data values too closely&#44; there are two <a name="p173"></a><a name="p174"></a>reasons&#46; First&#44; this &#8220;behavior&#8221; is observed in all simulation methods that were performed in order to compare the simulated values for the bivariate case of this proposal &#91;<a class="elsevierStyleCrossRef" href="#bib0110">Hernandez-Maldonado <span class="elsevierStyleItalic">et al</span>&#46;&#44; &#40;2012&#41;</a>&#93;&#46; In that work the most accurate results were shown by Bernstein copula&#46; In this particular case is hard to have a simulation that follow too closely the behavior of the real data values&#44; because there is very little information&#46; And second&#44; this is a median value of 10 different simulations&#44; which means that each simulation has different degrees of accuracy&#44; which is normal&#46;<a name="p175"></a><a name="p176"></a></p><p id="par0480" class="elsevierStylePara elsevierViewall"><a class="elsevierStyleCrossRef" href="#f0075">Figure 15</a> shows the variogram of this simulation and its fitted model&#46; In contrast to <a class="elsevierStyleCrossRef" href="#f0050">Figure 10</a> &#40;real data values&#41; we can see that this method represents very well&#44; the spatial structure of permeability&#46;</p><elsevierMultimedia ident="f0075"></elsevierMultimedia><p id="par0485" class="elsevierStylePara elsevierViewall">The values of the variogram model for a single simulation and the median of 10 simulations are very close to the real data as well &#40;<a class="elsevierStyleCrossRef" href="#tbl0015">Table 3</a>&#41;&#46; The same happens with their most important statistics &#40;<a class="elsevierStyleCrossRef" href="#f0080">Figure 16</a>&#41;&#46;</p><elsevierMultimedia ident="f0080"></elsevierMultimedia><p id="par0490" class="elsevierStylePara elsevierViewall">Using <a class="elsevierStyleCrossRef" href="#eq0080">&#40;16&#41;</a> we can see that the MSE values between the real data and these nonconditional simulations are very similar &#40;<a class="elsevierStyleCrossRef" href="#tbl0020">Table 4</a>&#41;&#46; It is an indicative that it is not necessary to perform <a name="p177"></a>many simulations to smooth the small-scale variability&#44; and may be unnecessary&#46;</p><elsevierMultimedia ident="tbl0020"></elsevierMultimedia><p id="par0495" class="elsevierStylePara elsevierViewall">Finally we perform conditional simulations and compare the MSE values in order to study the effect of different percentage of conditioning values&#46; <a class="elsevierStyleCrossRef" href="#f0085">Figure 17</a> shows the results of the spatial distribution of the permeability in various percentages of conditioning values&#44; from top to bottom order&#44; it is conditioned to 10 percent&#44; 50 percent and 90 percent of the conditioning data values&#46; Note that results are&#44; graphically&#44; very similar to real data and their MSE values as well &#40;see <a class="elsevierStyleCrossRef" href="#tbl0025">Table 5</a>&#41;&#46; This shows that the proposed method in a multivariate fashion does not require big amounts of conditioning data to have consistent results&#46;</p><elsevierMultimedia ident="f0085"></elsevierMultimedia><elsevierMultimedia ident="tbl0025"></elsevierMultimedia><p id="par0500" class="elsevierStylePara elsevierViewall">The Bernstein copula <a class="elsevierStyleCrossRef" href="#eq0040">&#40;8&#41;</a> requires that the empirical copula be ready to be used all the time&#44; hence&#44; it is necessary to propose an efficient and fast algorithm in terms of data storage and computing speed to implement a multidimensional copula&#46; We propose a procedure to efficiently implement a multidimensional empirical copula&#46; Shortly described&#44; this procedure fundamentally solves the two main disadvantages of using nonparametric copulas&#58; it avoids redundant calculation&#44; which increases its speed and just needs to storage the most important values of the generated hypercube&#46;</p><p id="par0505" class="elsevierStylePara elsevierViewall">Of course the method can be extended to more dimensions&#44; but the purpose of the article is to show its application to a well log scale&#44; where we have more control and better data quality&#46; Its application to more dimensions will be the subject of future work&#46;<a name="p178"></a><a name="p179"></a></p></span></span><span id="sec0065" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0085">Conclusions</span><p id="par0510" class="elsevierStylePara elsevierViewall">From a methodological point of view&#44; a multivariate copula approach provides a very flexible tool to model complex dependence relationships of petrophysical properties&#44; such as porosity&#44; permeability and shear wave velocity&#44; without imposing strong assumptions of linearity or log-linearity&#44; and&#47;or normality when we are modeling them as random variables&#46;</p><p id="par0515" class="elsevierStylePara elsevierViewall">All the information about the dependence structure is contained in the underlying copula&#44; and its estimation is being used&#44; instead of the extreme information reduction that is done by the use of numerical measures such as the linear correlation coefficient&#44; which under the presence of nonlinear dependence may become useless and&#47;or quite misleading&#44; see Embrechts <span class="elsevierStyleItalic">et al</span>&#46; &#40;<a class="elsevierStyleCrossRef" href="#bib0090">1999</a>&#44; <a class="elsevierStyleCrossRef" href="#bib0095">2003</a>&#41;&#46; For this reason it is very important to extend the method to higher dimensions&#44; because if we have more than one relevant explanatory variable&#44; we may use this tool to take them into account&#46;</p><p id="par0520" class="elsevierStylePara elsevierViewall">The main advantage of the methodology used in this work is that it represents a straightforward way to perform nonparametric simulations&#44; conditional and nonconditional&#46; As we have said in the results of the application of this method&#44; including more than one explanatory petrophysical property to leads to more accurate results&#46; Consequently&#44; it will not be necessary to perform many simulations &#40;<a class="elsevierStyleCrossRef" href="#f0030">Figure 6</a> vs&#46; <a class="elsevierStyleCrossRef" href="#f0060">Figure 12</a>&#41; or conditioning with high quantities of data values &#40;<a class="elsevierStyleCrossRef" href="#f0085">Figure 17</a>&#41; to obtain acceptable results &#40;<a class="elsevierStyleCrossRef" href="#tbl0025">Table 5</a>&#41;&#46;</p><p id="par0525" class="elsevierStylePara elsevierViewall">This method has an important improvement with respect to the bivariate case&#44; not only because it provides better results&#44; but also it reduces the number of simulations to decrease the small scale variability&#46; In other words it is not necessary to perform 10 or more simulations&#46; Even more&#44; it is not necessary to have a great number of conditioning values&#59; the conditioned simulations did not show a lot of variation in terms of MSE values&#46;</p><p id="par0530" class="elsevierStylePara elsevierViewall">For future work&#44; we think we can use a Bernstein copula to explore complex relationships between petrophysical properties and then fit&#44; if possible&#44; a semiparametric gluing copula&#44; <a class="elsevierStyleCrossRef" href="#bib0180">Siburg and Stoimenov &#40;2008&#41;</a>&#46; Also&#44; we consider we may use another optimization method which could give improved results&#44; for example&#44; we can implement genetic algorithms or Estimation of Distribution Algorithms &#40;EDAs&#41;&#46; It is very important to remark that the use of multivariate copulas leads to significant serious computational challenges&#44; for this reason it is necessary to propose better ways to calculate them&#44; for example parallel computing&#46;</p></span></span>"
    "textoCompletoSecciones" => array:1 [
      "secciones" => array:13 [
        0 => array:3 [
          "identificador" => "xres498003"
          "titulo" => "Resumen"
          "secciones" => array:1 [
            0 => array:1 [
              "identificador" => "abst0005"
            ]
          ]
        ]
        1 => array:2 [
          "identificador" => "xpalclavsec519550"
          "titulo" => "Palabras clave"
        ]
        2 => array:3 [
          "identificador" => "xres498004"
          "titulo" => "Abstract"
          "secciones" => array:1 [
            0 => array:1 [
              "identificador" => "abst0010"
            ]
          ]
        ]
        3 => array:2 [
          "identificador" => "xpalclavsec519551"
          "titulo" => "Key words"
        ]
        4 => array:2 [
          "identificador" => "sec0005"
          "titulo" => "Introduction"
        ]
        5 => array:2 [
          "identificador" => "sec0010"
          "titulo" => "Brief introduction to multivariate copulas"
        ]
        6 => array:3 [
          "identificador" => "sec0015"
          "titulo" => "Joint stochastic simulation method using a multivariate Bernstein copula"
          "secciones" => array:2 [
            0 => array:2 [
              "identificador" => "sec0020"
              "titulo" => "Multivariate Bernstein Copula Modeling"
            ]
            1 => array:2 [
              "identificador" => "sec0025"
              "titulo" => "Sampling algorithm of a Bernstein copula model"
            ]
          ]
        ]
        7 => array:2 [
          "identificador" => "sec1035"
          "titulo" => "Sampling bivariate algorithm"
        ]
        8 => array:3 [
          "identificador" => "sec2040"
          "titulo" => "Sampling trivariate algorithm"
          "secciones" => array:1 [
            0 => array:2 [
              "identificador" => "sec0030"
              "titulo" => "Simulated Annealing method"
            ]
          ]
        ]
        9 => array:3 [
          "identificador" => "sec0035"
          "titulo" => "Case study"
          "secciones" => array:5 [
            0 => array:2 [
              "identificador" => "sec0040"
              "titulo" => "Data description"
            ]
            1 => array:2 [
              "identificador" => "sec0045"
              "titulo" => "Statistical data analysis"
            ]
            2 => array:2 [
              "identificador" => "sec0050"
              "titulo" => "Bernstein Copula Model"
            ]
            3 => array:2 [
              "identificador" => "sec0055"
              "titulo" => "Sampling from Bernstein copula"
            ]
            4 => array:2 [
              "identificador" => "sec0060"
              "titulo" => "Permeability simulations and discussion"
            ]
          ]
        ]
        10 => array:2 [
          "identificador" => "sec0065"
          "titulo" => "Conclusions"
        ]
        11 => array:2 [
          "identificador" => "xack161016"
          "titulo" => "Acknowledgments"
        ]
        12 => array:1 [
          "titulo" => "References"
        ]
      ]
    ]
    "pdfFichero" => "main.pdf"
    "tienePdf" => true
    "fechaRecibido" => "2013-02-25"
    "fechaAceptado" => "2013-09-03"
    "PalabrasClave" => array:2 [
      "es" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Palabras clave"
          "identificador" => "xpalclavsec519550"
          "palabras" => array:6 [
            0 => "permeabilidad"
            1 => "porosidad"
            2 => "velocidad de onda de corte"
            3 => "dependencia multivariada"
            4 => "c&#243;pula de Bernstein"
            5 => "simulaci&#243;n geoestad&#237;stica"
          ]
        ]
      ]
      "en" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Key words"
          "identificador" => "xpalclavsec519551"
          "palabras" => array:6 [
            0 => "permeability"
            1 => "porosity"
            2 => "shear wave velocity"
            3 => "multivariate dependence"
            4 => "Bernstein copula"
            5 => "geostatistical simulation"
          ]
        ]
      ]
    ]
    "tieneResumen" => true
    "resumen" => array:2 [
      "es" => array:2 [
        "titulo" => "Resumen"
        "resumen" => "<span id="abst0005" class="elsevierStyleSection elsevierViewall"><p id="spar0030" class="elsevierStyleSimplePara elsevierViewall">En este trabajo se presenta un m&#233;todo no param&#233;trico general de simulaci&#243;n estoc&#225;stica conjunta de propiedades petrof&#237;sicas utilizando la c&#243;pula Bernstein&#46; Este m&#233;todo consiste b&#225;sicamente generar simulaciones estoc&#225;sticas de una determinada propiedad petrof&#237;sica &#40;variable primaria&#41; modelando la dependencia emp&#237;rica subyacente con otras propiedades petrof&#237;sicas &#40;variables secundarias&#41;&#44; mientras tambi&#233;n es reproducida la dependencia espacial de la primera&#46;</p><p id="spar0035" class="elsevierStyleSimplePara elsevierViewall">Este enfoque multivariado provee una herramienta muy flexible para modelar las complejas relaciones de dependencia de las propiedades petrof&#237;sicas&#46; Tiene varias ventajas sobre otros m&#233;todos tradicionales&#44; ya que no se limita al caso de la dependencia lineal entre las variables&#44; y tampoco requiere de la suposici&#243;n de normalidad y&#47;o existencia de momentos&#46;</p><p id="spar0040" class="elsevierStyleSimplePara elsevierViewall">En este trabajo este m&#233;todo es aplicado para simular un perfil de permeabilidad utilizando la porosidad vugular y velocidad de onda de corte &#40;Ondas S&#41; como covariables&#44; en una formaci&#243;n carbonatada de doble porosidad a escala de pozo&#46; Los valores simulados de la permeabilidad muestran un alto grado de precisi&#243;n en comparaci&#243;n con los valores reales&#46;</p></span>"
      ]
      "en" => array:2 [
        "titulo" => "Abstract"
        "resumen" => "<span id="abst0010" class="elsevierStyleSection elsevierViewall"><p id="spar0045" class="elsevierStyleSimplePara elsevierViewall">This paper introduces a general nonparametric method for joint stochastic simulation of petrophysical properties using the Bernstein copula&#46; This method consists basically in generating stochastic simulations of a given petrophysical property &#40;primary variable&#41; modeling the underlying empirical dependence with other petrophysical properties &#40;secondary variables&#41; while reproducing the spatial dependence of the first one&#46;</p><p id="spar0050" class="elsevierStyleSimplePara elsevierViewall">This multivariate approach provides a very flexible tool to model the complex dependence relationships of petrophysical properties&#46; It has several advantages over other traditional methods&#44; since it is not restricted to the case of linear dependence among variables&#44; it does not require the assumption of normality and&#47;or existence of moments&#46;</p><p id="spar0055" class="elsevierStyleSimplePara elsevierViewall">In this paper this method is applied to simulate rock permeability using Vugular Porosity and Shear Wave Velocity &#40;S-Waves&#41; as covariates in a carbonate double-porosity formation at well log scale&#46; Simulated permeability values show a high degree of accuracy compared to the actual values&#46;</p></span>"
      ]
    ]
    "multimedia" => array:38 [
      0 => array:7 [
        "identificador" => "f0005"
        "etiqueta" => "Figure 1"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr1.jpeg"
            "Alto" => 1239
            "Ancho" => 878
            "Tamanyo" => 138737
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0005" class="elsevierStyleSimplePara elsevierViewall">Diagram of the two-stage method to perform multivariate stochastic simulations using multivariate Bernstein copula&#46;</p>"
        ]
      ]
      1 => array:7 [
        "identificador" => "f0010"
        "etiqueta" => "Figure 2"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr2.jpeg"
            "Alto" => 711
            "Ancho" => 1417
            "Tamanyo" => 152727
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0010" class="elsevierStyleSimplePara elsevierViewall">Crack porosity &#40;PHICR&#41;&#44; vuggy porosity &#40;PHIV&#41;&#44; total porosity &#40;PHITOT&#41; and permeability &#40;K&#41;&#46;</p>"
        ]
      ]
      2 => array:7 [
        "identificador" => "f0015"
        "etiqueta" => "Figure 3"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr3.jpeg"
            "Alto" => 1078
            "Ancho" => 1080
            "Tamanyo" => 153450
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0015" class="elsevierStyleSimplePara elsevierViewall">Scatter plot distribution and histograms of PHIV and K of real data values&#46;</p>"
        ]
      ]
      3 => array:7 [
        "identificador" => "f0020"
        "etiqueta" => "Figure 4"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr4.jpeg"
            "Alto" => 1068
            "Ancho" => 1076
            "Tamanyo" => 132907
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0020" class="elsevierStyleSimplePara elsevierViewall">Scatter plot distribution and histograms of PHIV and K of real data values&#46;</p>"
        ]
      ]
      4 => array:7 [
        "identificador" => "f0025"
        "etiqueta" => "Figure 5"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr5.jpeg"
            "Alto" => 580
            "Ancho" => 1788
            "Tamanyo" => 164442
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0025" class="elsevierStyleSimplePara elsevierViewall">Scatter plot between Porosity and Permeability&#46; Left&#44; the 380 data values&#46; Center&#44; 380 simulated values using trivariate Bernstein copula&#46; Right&#44; 3800 simulated values&#46;</p>"
        ]
      ]
      5 => array:7 [
        "identificador" => "f0030"
        "etiqueta" => "Figure 6"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr6.jpeg"
            "Alto" => 695
            "Ancho" => 1441
            "Tamanyo" => 110315
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0030" class="elsevierStyleSimplePara elsevierViewall">Single nonconditional simulation of permeability &#40;K&#41;&#44; using porosity and shear wave velocity &#40;VS meas&#41; as secondary variables&#46;</p>"
        ]
      ]
      6 => array:7 [
        "identificador" => "f0035"
        "etiqueta" => "Figure 7"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr7.jpeg"
            "Alto" => 1073
            "Ancho" => 1073
            "Tamanyo" => 155056
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0035" class="elsevierStyleSimplePara elsevierViewall">Scatter plot and histograms of PHIV and K &#40;Simulated values&#41;&#46;</p>"
        ]
      ]
      7 => array:7 [
        "identificador" => "f0040"
        "etiqueta" => "Figure 8"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr8.jpeg"
            "Alto" => 1081
            "Ancho" => 1087
            "Tamanyo" => 130538
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0040" class="elsevierStyleSimplePara elsevierViewall">Scatter plot and histograms of VS Meas and K &#40;Simulated values&#41;&#46;</p>"
        ]
      ]
      8 => array:7 [
        "identificador" => "f0045"
        "etiqueta" => "Figure 9"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr9.jpeg"
            "Alto" => 1116
            "Ancho" => 1113
            "Tamanyo" => 105802
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0045" class="elsevierStyleSimplePara elsevierViewall">Estimated variogram and its fitted model for simulated K &#40;Simulated values&#41;&#46;</p>"
        ]
      ]
      9 => array:7 [
        "identificador" => "f0050"
        "etiqueta" => "Figure 10"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr10.jpeg"
            "Alto" => 1118
            "Ancho" => 1116
            "Tamanyo" => 104997
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0050" class="elsevierStyleSimplePara elsevierViewall">Estimated variogram and its fitted model for original K &#40;Data values&#41;&#46;</p>"
        ]
      ]
      10 => array:7 [
        "identificador" => "f0055"
        "etiqueta" => "Figure 11"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr11.jpeg"
            "Alto" => 679
            "Ancho" => 1184
            "Tamanyo" => 83653
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0055" class="elsevierStyleSimplePara elsevierViewall">Statistical comparison of to the original data and a single trivariate SA simulation&#46;</p>"
        ]
      ]
      11 => array:7 [
        "identificador" => "f0060"
        "etiqueta" => "Figure 12"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr12.jpeg"
            "Alto" => 666
            "Ancho" => 1449
            "Tamanyo" => 104973
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0060" class="elsevierStyleSimplePara elsevierViewall">Median of 10 nonconditional simulation of permeability &#40;K&#41;&#44; using porosity and shear wave velocity &#40;S-waves&#41; as secondary variables&#46;</p>"
        ]
      ]
      12 => array:7 [
        "identificador" => "f0065"
        "etiqueta" => "Figure 13"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr13.jpeg"
            "Alto" => 1080
            "Ancho" => 1088
            "Tamanyo" => 159213
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0065" class="elsevierStyleSimplePara elsevierViewall">Scatter plot and histograms of PHIV and K simulated values &#40;median of 10 simulations&#41;&#46;</p>"
        ]
      ]
      13 => array:7 [
        "identificador" => "f0070"
        "etiqueta" => "Figure 14"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr14.jpeg"
            "Alto" => 1083
            "Ancho" => 1085
            "Tamanyo" => 133757
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0070" class="elsevierStyleSimplePara elsevierViewall">Scatter plot distribution and histograms of VS-meas and K simulated values &#40;a median of 10 simulations&#41;&#46;</p>"
        ]
      ]
      14 => array:7 [
        "identificador" => "f0075"
        "etiqueta" => "Figure 15"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr15.jpeg"
            "Alto" => 1113
            "Ancho" => 1107
            "Tamanyo" => 107974
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0075" class="elsevierStyleSimplePara elsevierViewall">Estimated variogram and its fitted model for simulated K data values&#46;</p>"
        ]
      ]
      15 => array:7 [
        "identificador" => "f0080"
        "etiqueta" => "Figure 16"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr16.jpeg"
            "Alto" => 648
            "Ancho" => 1179
            "Tamanyo" => 81663
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0080" class="elsevierStyleSimplePara elsevierViewall">Statistical comparison of to the original data&#44; and the median of 10 SA simulations&#46;</p>"
        ]
      ]
      16 => array:7 [
        "identificador" => "f0085"
        "etiqueta" => "Figure 17"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr17.jpeg"
            "Alto" => 1671
            "Ancho" => 1449
            "Tamanyo" => 269539
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="sp0085" class="elsevierStyleSimplePara elsevierViewall">Spatial distribution of permeability in single conditional simulations using PHIV and VS meas as secondary variables&#46; Conditioning levels &#40;in top to bottom order&#41; 10 percent&#44; 50 percent and 90 percent&#46;</p>"
        ]
      ]
      17 => array:7 [
        "identificador" => "tbl0005"
        "etiqueta" => "Table 1"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Explanatory variable&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Corresponding Hoeffding index <span class="elsevierStyleItalic">&#934;</span> vs&#46; K&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHIV&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;71&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHITOT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;62&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHIMAT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;50&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">VS Meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;60&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">VP meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;55&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHICR&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;16&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab795073.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0005" class="elsevierStyleSimplePara elsevierViewall">Dependence index <span class="elsevierStyleItalic">&#934;</span>&#44; <a class="elsevierStyleCrossRef" href="#bib0115">Hoeffding &#40;1941&#41;</a>&#44; for possible explanatory variables &#40;PHITOT&#44; PHIMAT&#44; PHICR&#44; PHIV&#44; VS Meas&#44; VP meas&#41;&#46;</p>"
        ]
      ]
      18 => array:7 [
        "identificador" => "tbl0010"
        "etiqueta" => "Table 2"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Hoeffding index <span class="elsevierStyleItalic">&#934;</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHIV&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHITOT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHIMAT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">PHICR&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">VS Meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">VP meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHIV&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;88&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;72&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;34&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;58&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;57&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHITOT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;43&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;32&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;64&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;61&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHIMAT&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;24&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;30&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;36&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">PHICR&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;31&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;20&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">VS Meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;32&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">VP meas&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">1&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab795072.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0010" class="elsevierStyleSimplePara elsevierViewall">Hoeffding dependence index <span class="elsevierStyleItalic">&#934;</span> matrix&#44; for possible explanatory variables &#40;PHITOT&#44; PHIMAT&#44; PHICR&#44; PHIV&#44; VS Meas&#44; VP meas&#41;&#46;</p>"
        ]
      ]
      19 => array:7 [
        "identificador" => "tbl0015"
        "etiqueta" => "Table 3"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Nugget&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Sill&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">Range&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">Data values&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">2&#46;61 x 106&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">85&#46;07&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">Single SA nonconditional Simulation&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">2&#46;57 x 106&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">87&#46;52&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="center" valign="top">Median of 10 SA nonconditional Simulations&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">2&#46;51 x 106&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">87&#46;27&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab795070.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0015" class="elsevierStyleSimplePara elsevierViewall">Variograms Models of the data values and single nonconditional SA simulation&#46;</p>"
        ]
      ]
      20 => array:7 [
        "identificador" => "tbl0020"
        "etiqueta" => "Table 4"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">MSE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Single SA nonconditional Sim&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;50&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Median of 10 nonconditionals Sims&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;35&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab795071.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0020" class="elsevierStyleSimplePara elsevierViewall">MSE values of nonconditional simulations</p>"
        ]
      ]
      21 => array:7 [
        "identificador" => "tbl0025"
        "etiqueta" => "Table 5"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Configuration&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="center" valign="top" scope="col" style="border-bottom: 2px solid black">MSE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Single SA nonconditional Sim&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;50&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Single SA 10&#37; conditional Sim&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;41&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Single SA 50&#37; conditional Sim&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;36&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Single SA 90&#37; conditional Sim&#46;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="center" valign="top">0&#46;27&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab795069.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0025" class="elsevierStyleSimplePara elsevierViewall">MSE values of the single conditional simulations&#46;</p>"
        ]
      ]
      22 => array:6 [
        "identificador" => "eq0005"
        "etiqueta" => "&#40;1&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "H&#40;x1&#44;x2&#44;&#8230;&#44;xm&#41;&#61;C&#40;F2&#40;x1&#41;&#44;F2&#40;x&#41;&#44;&#8230;&#44;Fm&#40;xm&#41;&#41;"
          "Fichero" => "si1.jpeg"
          "Tamanyo" => 2642
          "Alto" => 14
          "Ancho" => 338
        ]
      ]
      23 => array:6 [
        "identificador" => "eq0010"
        "etiqueta" => "&#40;2&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "F&#710;j&#40;x&#41;&#61;1n&#8721;k&#61;1n1&#123;xjk&#8804;x&#125;"
          "Fichero" => "si2.jpeg"
          "Tamanyo" => 2047
          "Alto" => 46
          "Ancho" => 160
        ]
      ]
      24 => array:6 [
        "identificador" => "eq0015"
        "etiqueta" => "&#40;3&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "Cn&#40;v1n&#44;&#8230;&#44;vmn&#41;&#61;1nk&#61;1n1&#123;rank&#40;x1k&#41;&#163;v1&#44;&#8230;&#44;rank&#40;xmk&#41;&#8804;vm&#125;"
          "Fichero" => "si6.jpeg"
          "Tamanyo" => 3587
          "Alto" => 22
          "Ancho" => 414
        ]
      ]
      25 => array:6 [
        "identificador" => "eq0020"
        "etiqueta" => "&#40;4&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "H&#40;x1&#44;x2&#44;x3&#41;&#61;C&#40;F1&#40;x1&#41;&#44;F2&#40;x2&#41;&#44;F3&#40;x3&#41;&#41;"
          "Fichero" => "si7.jpeg"
          "Tamanyo" => 2437
          "Alto" => 14
          "Ancho" => 276
        ]
      ]
      26 => array:6 [
        "identificador" => "eq0025"
        "etiqueta" => "&#40;5&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "Q&#732;&#59;n&#40;u&#41;&#61;&#8721;k&#61;0n12&#40;xk&#43;xk&#43;1&#41;nkuk&#40;1-u&#41;n-k"
          "Fichero" => "si9.jpeg"
          "Tamanyo" => 3586
          "Alto" => 49
          "Ancho" => 301
        ]
      ]
      27 => array:6 [
        "identificador" => "eq0030"
        "etiqueta" => "&#40;6&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "C&#732;&#59;B&#40;u1&#44;&#8230;&#44;um&#41;&#61;&#8721;v1&#61;0n1&#8230;&#8721;vm&#61;0nm&#945;v1n1&#44;&#8230;&#44;vmnmPv1&#44;m1&#40;u1&#41;&#8230;Pvm&#44;nm&#40;um&#41;"
          "Fichero" => "si10.jpeg"
          "Tamanyo" => 4854
          "Alto" => 49
          "Ancho" => 484
        ]
      ]
      28 => array:6 [
        "identificador" => "eq0035"
        "etiqueta" => "&#40;7&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "Pvj&#44;mj&#40;uj&#41;&#61;njvjujvj&#40;1-uj&#41;nj-vj"
          "Fichero" => "si11.jpeg"
          "Tamanyo" => 2735
          "Alto" => 49
          "Ancho" => 243
        ]
      ]
      29 => array:6 [
        "identificador" => "eq0040"
        "etiqueta" => "&#40;8&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "C&#732;&#59;B&#40;u&#44;v&#44;w&#41;&#61;&#8721;i&#61;0n&#8721;j&#61;0n&#8721;k&#61;0nCnin&#44;jn&#44;knniui&#40;1-u&#41;n-injvj&#40;1-v&#41;n-jnkwk&#40;1-w&#41;n-k"
          "Fichero" => "si13.jpeg"
          "Tamanyo" => 8387
          "Alto" => 52
          "Ancho" => 647
        ]
      ]
      30 => array:6 [
        "identificador" => "eq0045"
        "etiqueta" => "&#40;9&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "cu&#40;v&#41;&#61;&#8706;C&#732;&#59;B&#40;u&#44;v&#41;&#8706;u"
          "Fichero" => "si14.jpeg"
          "Tamanyo" => 1178
          "Alto" => 24
          "Ancho" => 106
        ]
      ]
      31 => array:6 [
        "identificador" => "eq0050"
        "etiqueta" => "&#40;10&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "cuv&#40;w&#41;&#61;&#8706;C&#732;&#59;B&#40;u&#44;v&#44;w&#41;&#8706;u&#8706;v&#8706;C&#732;&#59;B&#40;u&#44;v&#44;1&#41;&#8706;u&#8706;v"
          "Fichero" => "si19.jpeg"
          "Tamanyo" => 1884
          "Alto" => 42
          "Ancho" => 122
        ]
      ]
      32 => array:6 [
        "identificador" => "eq0055"
        "etiqueta" => "&#40;11&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "FO&#61;&#8721;i&#91;&#947;&#42;&#40;hi&#41;-&#947;&#40;hi&#41;&#93;2&#947;&#40;hi&#41;2"
          "Fichero" => "si23.jpeg"
          "Tamanyo" => 1911
          "Alto" => 40
          "Ancho" => 151
        ]
      ]
      33 => array:6 [
        "identificador" => "eq0060"
        "etiqueta" => "&#40;12&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "&#947;&#40;h&#41;&#61;12N&#40;h&#41;&#8721;i&#61;1N&#40;h&#41;&#91;Z&#40;xi&#43;h&#41;-Z&#40;xi&#41;&#93;2"
          "Fichero" => "si24.jpeg"
          "Tamanyo" => 2790
          "Alto" => 50
          "Ancho" => 248
        ]
      ]
      34 => array:6 [
        "identificador" => "eq0065"
        "etiqueta" => "&#40;13&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "T0&#61;&#40;-&#916;Elog&#40;&#964;0&#41;&#41;"
          "Fichero" => "si25.jpeg"
          "Tamanyo" => 1078
          "Alto" => 23
          "Ancho" => 93
        ]
      ]
      35 => array:6 [
        "identificador" => "eq0070"
        "etiqueta" => "&#40;14&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "MSE&#40;logK&#124;PHIV&#41;&#61;1n&#8721;m&#61;1n&#91;zm-&#946;0&#46;5&#40;xm&#41;&#93;2"
          "Fichero" => "si27.jpeg"
          "Tamanyo" => 3414
          "Alto" => 46
          "Ancho" => 311
        ]
      ]
      36 => array:6 [
        "identificador" => "eq0075"
        "etiqueta" => "&#40;15&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "MSE&#40;logK&#124;VSmeas&#41;&#61;1n&#8721;m&#61;1n&#91;zm-&#946;0&#46;5&#40;ym&#41;&#93;2"
          "Fichero" => "si28.jpeg"
          "Tamanyo" => 3613
          "Alto" => 46
          "Ancho" => 334
        ]
      ]
      37 => array:6 [
        "identificador" => "eq0080"
        "etiqueta" => "&#40;16&#41;"
        "tipo" => "MULTIMEDIAFORMULA"
        "mostrarFloat" => false
        "mostrarDisplay" => true
        "Formula" => array:5 [
          "Matematica" => "MSE&#40;logK&#124;PHIV&#44;VSmeas&#41;&#61;1n&#8721;m&#61;1n&#91;zm-&#946;0&#46;5&#40;xm&#44;ym&#41;&#93;2"
          "Fichero" => "si29.jpeg"
          "Tamanyo" => 4305
          "Alto" => 46
          "Ancho" => 399
        ]
      ]
    ]
    "bibliografia" => array:2 [
      "titulo" => "References"
      "seccion" => array:1 [
        0 => array:2 [
          "identificador" => "bibs0005"
          "bibliografiaReferencia" => array:38 [
            0 => array:3 [
              "identificador" => "bib0005"
              "etiqueta" => "Al-harthy et al&#46;&#44; 2005"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Copulas&#58; A new technique to model dependence in petroleum decision making"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "M&#46; Al-harthy"
                            1 => "S&#46; Begg"
                            2 => "R&#46; Bratvold"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:2 [
                      "doi" => "10.1038/srep09602"
                      "Revista" => array:6 [
                        "tituloSerie" => "Journal of Petroleum Science &#38; Engineering"
                        "fecha" => "2005"
                        "volumen" => "57"
                        "paginaInicial" => "195"
                        "paginaFinal" => "208"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            1 => array:3 [
              "identificador" => "bib0010"
              "etiqueta" => "Armstrong et al&#46;&#44; 2004"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Incorporating technical uncertainty in real option valuation of oil projects"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "M&#46; Armstrong"
                            1 => "A&#46; Galli"
                            2 => "W&#46; Bailey"
                            3 => "B&#46; Couet"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:2 [
                      "doi" => "10.1038/srep09602"
                      "Revista" => array:6 [
                        "tituloSerie" => "Journal of Petroleum Science &#38; Engineering"
                        "fecha" => "2004"
                        "volumen" => "44"
                        "paginaInicial" => "67"
                        "paginaFinal" => "82"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            2 => array:3 [
              "identificador" => "bib0015"
              "etiqueta" => "Balan et al&#46;&#44; 1995"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "State-of-The-Art in Permeability Determination From Well Log Data&#58; Part 1&#46;- A Comparative Study"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "B&#46; Balan"
                            1 => "S&#46; Mohaghegh"
                            2 => "S&#46; Ameri"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "Model Development"
                        "fecha" => "1995"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            3 => array:3 [
              "identificador" => "bib0020"
              "etiqueta" => "Bennett et al&#46;&#44; 2002"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Hydrologic investigation of the Floridian aquifer system&#44; western Hillsboro Basin&#44; Palm Beach County&#44; Florida&#44; Technical Publication WS-8"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "W&#46;M&#46; Bennett"
                            1 => "P&#46;F&#46; Linton"
                            2 => "E&#46;E&#46; Rectenwald"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "South Florida Water Management District"
                        "fecha" => "2002"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            4 => array:3 [
              "identificador" => "bib0025"
              "etiqueta" => "Billingsley&#44; 1995"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "P&#46; Billingsley"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Probability and Measure"
                        "fecha" => "1995"
                        "editorial" => "Wiley"
                        "editorialLocalizacion" => "New York"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            5 => array:3 [
              "identificador" => "bib0030"
              "etiqueta" => "Cosentino&#44; 2001"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "L&#46; Cosentino"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "Integrated Reservoir Studies"
                        "fecha" => "2001"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            6 => array:3 [
              "identificador" => "bib0035"
              "etiqueta" => "Dafflon et al&#46;&#44; 2009a"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "B&#46; Dafflon"
                            1 => "J&#46; Irving"
                            2 => "K&#46; Holliger"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Applied Geophysics"
                        "fecha" => "2009"
                        "volumen" => "68"
                        "paginaInicial" => "60"
                        "paginaFinal" => "70"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            7 => array:3 [
              "identificador" => "bib0040"
              "etiqueta" => "Dafflon et al&#46;&#44; 2009b"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Quantitative Integration of High-Resolution Hydrogeophysical Data&#58; A Novel Approach to Monte-Carlo-Type Conditional Stochastic Simulations and Implications for Hydrological Predictions"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "B&#46; Dafflon"
                            1 => "J&#46; Irving"
                            2 => "K&#46; Holliger&#42;"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Earth Science"
                        "fecha" => "2009"
                        "volumen" => "20"
                        "paginaInicial" => "580"
                        "paginaFinal" => "591"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            8 => array:3 [
              "identificador" => "bib0045"
              "etiqueta" => "Dafflon and Barrash&#44; 2012"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Three-dimensional stochastic estimation of porosity distribution&#58; Benefits of using ground-penetrating radar velocity tomograms in simulated-annealing-based or Bayesian sequential simulation approaches"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "B&#46; Dafflon"
                            1 => "W&#46; Barrash"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:4 [
                        "tituloSerie" => "Water Resources Research"
                        "fecha" => "2012"
                        "volumen" => "48"
                        "paginaInicial" => "5"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            9 => array:3 [
              "identificador" => "bib0050"
              "etiqueta" => "Accioly and Chiyoshi&#44; 2004"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Modeling dependence with copulas&#58; a useful tool for field development decision process"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "R&#46; Accioly"
                            1 => "F&#46; Chiyoshi"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:2 [
                      "doi" => "10.1038/srep09602"
                      "Revista" => array:6 [
                        "tituloSerie" => "Journal of Petroleum Science &#38; Engineering"
                        "fecha" => "2004"
                        "volumen" => "44"
                        "paginaInicial" => "83"
                        "paginaFinal" => "91"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/25946319"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            10 => array:3 [
              "identificador" => "bib0055"
              "etiqueta" => "Deheuvels&#44; 1979"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "La fonction de d&#233;pendance empirique et ses porpri&#233;t&#233;s&#46; Un test nonparam&#233;trique d&#8217;independance&#44; Acad"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "P&#46; Deheuvels"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:6 [
                        "tituloSerie" => "Roy&#46; Belg&#46; Bull&#46; Cl&#46; Sci&#46; &#40;5&#41;"
                        "fecha" => "1979"
                        "volumen" => "65"
                        "numero" => "6"
                        "paginaInicial" => "274"
                        "paginaFinal" => "292"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            11 => array:3 [
              "identificador" => "bib0060"
              "etiqueta" => "Deutsch&#44; 1992"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Annealing Techniques Applied to Reservoir Modeling and the Integration of Geological and Engineering &#40;Well Test&#41; Data"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "C&#46;V&#46; Deutsch"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "PhD thesis"
                        "fecha" => "1992"
                        "editorial" => "Stanford University"
                        "editorialLocalizacion" => "Stanford&#44; CA"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            12 => array:3 [
              "identificador" => "bib0065"
              "etiqueta" => "Deutsch and Cockerham&#44; 1994"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "C&#46;V&#46; Deutsch"
                            1 => "P&#46;W&#46; Cockerham"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "Geostatistical Modeling of Permeability with Annealing Cosimulation &#40;ACS&#41;"
                        "fecha" => "1994"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            13 => array:3 [
              "identificador" => "bib0070"
              "etiqueta" => "D&#237;az-Viera and Casar-Gonz&#225;lez&#44; 2005"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Stochastic simulation of complex dependency patterns of petrophysical properties using t-copulas"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "M&#46; D&#237;az-Viera"
                            1 => "R&#46; Casar-Gonz&#225;lez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Proc&#46; IAMG&#8217;05&#58; GIS and Spatial Analysis"
                        "fecha" => "2005"
                        "volumen" => "2"
                        "paginaInicial" => "749"
                        "paginaFinal" => "755"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            14 => array:3 [
              "identificador" => "bib0075"
              "etiqueta" => "D&#237;az-Viera et al&#46;&#44; 2006"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Stochastic modeling of permeability in double porosity carbonates applying a Monte-Carlo simulation method with t-copula"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:5 [
                            0 => "M&#46; D&#237;az-Viera"
                            1 => "P&#46; Anguiano-Rojas"
                            2 => "A&#46; Mousatov"
                            3 => "E&#46; Kazatchenko"
                            4 => "M&#46; Markov"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "SPWLA 47th Annual Logging Symposium"
                        "fecha" => "2006"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            15 => array:3 [
              "identificador" => "bib0080"
              "etiqueta" => "Dr&#233;o et al&#46;&#44; 2006"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "J&#46; Dr&#233;o"
                            1 => "A&#46; P&#233;trowski"
                            2 => "P&#46; Siarry"
                            3 => "T&#46; Taillard"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Metaheuristics for Hard Optimization"
                        "fecha" => "2006"
                        "editorial" => "Springer-Verlag Berlin Heidelberg"
                        "editorialLocalizacion" => "Germany"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            16 => array:3 [
              "identificador" => "bib0085"
              "etiqueta" => "Erdely and D&#237;az-Viera&#44; 2009"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Nonparametric and semiparametric bivariate modeling of petrophysical porosity-permeability well log data"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "A&#46; Erdely"
                            1 => "M&#46; D&#237;az-Viera"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Copula Theory and Its Applications &#40;P Jaworski <span class="elsevierStyleItalic">et al</span>&#46; eds&#41;&#44; Lecture Notes in Statistics &#8211; Springer"
                        "fecha" => "2009"
                        "volumen" => "198"
                        "paginaInicial" => "267"
                        "paginaFinal" => "278"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            17 => array:3 [
              "identificador" => "bib0090"
              "etiqueta" => "Embrechts et al&#46;&#44; 1999"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Correlation&#58; pitfalls and alternatives"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "P&#46; Embrechts"
                            1 => "A&#46; McNeil"
                            2 => "D&#46; Straumann"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Risk Magazine"
                        "fecha" => "1999"
                        "volumen" => "5"
                        "paginaInicial" => "69"
                        "paginaFinal" => "71"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            18 => array:3 [
              "identificador" => "bib0095"
              "etiqueta" => "Embrechts et al&#46;&#44; 2003"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Modeling dependence with copulas and applications to risk management"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "P&#46; Embrechts"
                            1 => "F&#46; Lindskog"
                            2 => "A&#46;J&#46; McNeil"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "LibroEditado" => array:4 [
                        "titulo" => "Handbook of Heavy-Tailed Distributions in Finance"
                        "paginaInicial" => "329"
                        "paginaFinal" => "384"
                        "serieFecha" => "2003"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            19 => array:3 [
              "identificador" => "bib0100"
              "etiqueta" => "Fermanian et al&#46;&#44; 2004"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Weak convergence of empirical copula processes"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "J-D&#46; Fermanian"
                            1 => "D&#46; Radulovic"
                            2 => "M&#46; Wegcamp"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Bernoulli"
                        "fecha" => "2004"
                        "volumen" => "10"
                        "paginaInicial" => "547"
                        "paginaFinal" => "560"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            20 => array:3 [
              "identificador" => "bib0105"
              "etiqueta" => "Haldorsen and Damsleth&#44; 1990"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Stochastic modeling"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "H&#46; Haldorsen"
                            1 => "E&#46; Damsleth"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Journal of Petroleum Technology"
                        "fecha" => "1990"
                        "paginaInicial" => "404"
                        "paginaFinal" => "412"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            21 => array:3 [
              "identificador" => "bib0110"
              "etiqueta" => "Hern&#225;ndez-Maldonado et al&#46;&#44; 2012"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "A joint stochastic simulation method using the Bernstein copula as a flexible tool for modeling nonlinear dependence structures between petrophysical properties"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "V&#46; Hern&#225;ndez-Maldonado"
                            1 => "M&#46; D&#237;az-Viera"
                            2 => "A&#46; Erdely"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Journal of Petroleum Science &#38; Engineering"
                        "fecha" => "2012"
                        "paginaInicial" => "90"
                        "paginaFinal" => "91"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            22 => array:3 [
              "identificador" => "bib0115"
              "etiqueta" => "Hoeffding&#44; 1941"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Scale-invariant correlation theory"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "W&#46; Hoeffding"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "LibroEditado" => array:4 [
                        "titulo" => "The Collected Works of Wassily Hoeffding"
                        "paginaInicial" => "57"
                        "paginaFinal" => "107"
                        "serieFecha" => "1941"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            23 => array:3 [
              "identificador" => "bib0120"
              "etiqueta" => "Holden et al&#46;&#44; 1995"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Use of Well Test Data in Stochastic Reservoir Modelling"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:6 [
                            0 => "L&#46; Holden"
                            1 => "R&#46; Madsen"
                            2 => "A&#46; Skorstad"
                            3 => "K&#46;A&#46; Jakobsen"
                            4 => "C&#46;B&#46; Tjolsen"
                            5 => "S&#46; Vik"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "SPE annual Technical Conference and Exhibition held in Dallas"
                        "fecha" => "1995"
                        "paginaInicial" => "22"
                        "paginaFinal" => "25"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            24 => array:3 [
              "identificador" => "bib0125"
              "etiqueta" => "Journel and Alabert&#44; 1990"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "New method for reservoir mapping"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "A&#46; Journel"
                            1 => "F&#46; Alabert"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Journal of Petroleum Technology"
                        "fecha" => "1990"
                        "paginaInicial" => "212"
                        "paginaFinal" => "218"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            25 => array:3 [
              "identificador" => "bib0130"
              "etiqueta" => "Journel and Zhu&#44; 1990"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Integrating soft seismic data&#58; Markov-Bayes updating&#44; an alternative to cokriging and traditional regression"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "A&#46; Journel"
                            1 => "H&#46; Zhu"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:2 [
                        "titulo" => "In Report 3&#44; Stanford Center for Reservoir Forecasting"
                        "fecha" => "1990"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            26 => array:3 [
              "identificador" => "bib0135"
              "etiqueta" => "Kazatchenko et al&#46;&#44; 2004"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Joint inversion of acoustic and resistivity data for carbonate microstructure evaluation"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "E&#46; Kazatchenko"
                            1 => "M&#46; Markov"
                            2 => "A&#46; Mousatov"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Petrophysics"
                        "fecha" => "2004"
                        "volumen" => "45"
                        "paginaInicial" => "130"
                        "paginaFinal" => "140"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            27 => array:3 [
              "identificador" => "bib0140"
              "etiqueta" => "Kazatchenko et al&#46;&#44; 2006"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Carbonate microstructure determination by inversion of acoustic and electrical data&#58; application to a South Florida Aquifer"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "E&#46; Kazatchenko"
                            1 => "M&#46; Markov"
                            2 => "A&#46; Mousatov"
                            3 => "J&#46; Parra"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "J&#46; Appl&#46; Geophys"
                        "fecha" => "2006"
                        "volumen" => "59"
                        "paginaInicial" => "1"
                        "paginaFinal" => "15"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            28 => array:3 [
              "identificador" => "bib0145"
              "etiqueta" => "Landa et al&#46;&#44; 1996"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Reservoir Characterization Constrained to Well Test Data&#58; A Field Example"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "J&#46;L&#46; Landa"
                            1 => "M&#46;M&#46; Kamal"
                            2 => "C&#46;D&#46; Jenkins"
                            3 => "R&#46;N&#46; Horne"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "SPE Annual Technical Conference and Exhibition held in Denver"
                        "fecha" => "1996"
                        "paginaInicial" => "6"
                        "paginaFinal" => "9"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            29 => array:3 [
              "identificador" => "bib0150"
              "etiqueta" => "Mu&#241;oz-Perez and Fern&#225;ndez-Palac&#237;n&#44; 1987"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Estimating the quantile function by Bernstein polynomials"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "J&#46; Mu&#241;oz-Perez"
                            1 => "A&#46; Fern&#225;ndez-Palac&#237;n"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Computational Statistics and Data Analysis"
                        "fecha" => "1987"
                        "volumen" => "5"
                        "paginaInicial" => "391"
                        "paginaFinal" => "397"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            30 => array:3 [
              "identificador" => "bib0155"
              "etiqueta" => "Nelsen&#44; 2006"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "R&#46;B&#46; Nelsen"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "An Introduction to Copulas"
                        "fecha" => "2006"
                        "editorial" => "Springer"
                        "editorialLocalizacion" => "New York"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            31 => array:3 [
              "identificador" => "bib0160"
              "etiqueta" => "Parra et al&#46;&#44; 2001"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "A methodology to integrate magnetic resonance and acoustic measurements for reservoir characterization"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "J&#46;O&#46; Parra"
                            1 => "C&#46;L&#46; Hackert"
                            2 => "H&#46;A&#46; Collier"
                            3 => "M&#46; Bennett"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "Report DOE&#47;BC&#47;15203-3"
                        "fecha" => "2001"
                        "editorial" => "National Petroleum Technology Office&#44; Department of Energy"
                        "editorialLocalizacion" => "Tulsa&#44; Oklahoma"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            32 => array:3 [
              "identificador" => "bib0165"
              "etiqueta" => "Parra and Hackert&#44; 2002"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Permeability and porosity images based on crosswell reflection seismic measurements of a vuggy carbonate aquifer at the Hillsboro site"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "J&#46;O&#46; Parra"
                            1 => "C&#46;L&#46; Hackert"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:3 [
                        "titulo" => "South Florida 72nd Annual Meeting of the Society of Exploration Geophysicists"
                        "fecha" => "2002"
                        "editorial" => "Paper VCD P1&#46;2"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            33 => array:3 [
              "identificador" => "bib0170"
              "etiqueta" => "Sancetta and Satchell&#44; 2004"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "The Bernstein copula and its applications to modeling and approximations of multivariate distributions"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "A&#46; Sancetta"
                            1 => "S&#46; Satchell"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Econometric Theory"
                        "fecha" => "2004"
                        "volumen" => "20"
                        "paginaInicial" => "535"
                        "paginaFinal" => "562"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            34 => array:3 [
              "identificador" => "bib0175"
              "etiqueta" => "Sancetta&#44; 2007"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Nonparametric estimation of distributions with given marginals via Bernstein-Kantorovic polynomials&#58; L1 and pointwise convergence theory"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "A&#46; Sancetta"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:6 [
                        "tituloSerie" => "Journal of Multivariate Analysis"
                        "fecha" => "2007"
                        "volumen" => "98"
                        "numero" => "7"
                        "paginaInicial" => "1376"
                        "paginaFinal" => "1390"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            35 => array:3 [
              "identificador" => "bib0180"
              "etiqueta" => "Siburg and Stoimenov&#44; 2008"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "A scalar product for copulas"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "K&#46; Siburg"
                            1 => "P&#46; Stoimenov"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Mathematical Analysis and Applications"
                        "fecha" => "2008"
                        "volumen" => "344"
                        "paginaInicial" => "429"
                        "paginaFinal" => "439"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            36 => array:3 [
              "identificador" => "bib0185"
              "etiqueta" => "Suro-Perez&#44; 1992"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Indicator Kriging and Simulation Based on Principal Component Analysis"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "V&#46; Suro-Perez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:4 [
                        "titulo" => "PhD thesis"
                        "fecha" => "1992"
                        "editorial" => "Stanford University"
                        "editorialLocalizacion" => "Stanford&#44; CA"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            37 => array:3 [
              "identificador" => "bib0190"
              "etiqueta" => "Sklar&#44; 1959"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:1 [
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "A&#46; Sklar"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "LibroEditado" => array:5 [
                        "paginaInicial" => "229"
                        "paginaFinal" => "231"
                        "serieVolumen" => "8"
                        "serieTitulo" => "Fonctions de r&#233;partition &#225; n dimensions et leurs marges"
                        "serieFecha" => "1959"
                      ]
                    ]
                  ]
                ]
              ]
            ]
          ]
        ]
      ]
    ]
    "agradecimientos" => array:1 [
      0 => array:4 [
        "identificador" => "xack161016"
        "titulo" => "Acknowledgments"
        "texto" => "<p id="par0535" class="elsevierStylePara elsevierViewall">The present work was partially supported by the two following projects&#58; IN110311 from Programa de Apoyo a Proyectos de Investigaci&#243;n e Innovaci&#243;n Tecnol&#243;gica &#40;PAPIIT&#41; at Universidad Nacional Aut&#243;noma de M&#233;xico&#44; and by 116606 &#40;Y&#46;00102&#41; from Fondo Sectorial CONACYT-SENER-Hidrocarburos&#46;</p>"
        "vista" => "all"
      ]
    ]
  ]
  "idiomaDefecto" => "en"
  "url" => "/00167169/0000005300000002/v2_201505081358/S0016716914714989/v2_201505081358/en/main.assets"
  "Apartado" => array:4 [
    "identificador" => "36047"
    "tipo" => "SECCION"
    "es" => array:2 [
      "titulo" => "Original paper"
      "idiomaDefecto" => true
    ]
    "idiomaDefecto" => "es"
  ]
  "PDF" => "https://static.elsevier.es/multimedia/00167169/0000005300000002/v2_201505081358/S0016716914714989/v2_201505081358/en/main.pdf?idApp=UINPBA00004N&text.app=https://www.elsevier.es/"
  "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716914714989?idApp=UINPBA00004N"
]

Article information

ISSN: 00167169
Original language: English

DOI:10.1016/S0016-7169(14)71498-9

The statistics are updated each day

Year/Month	Html	Pdf	Total

2024 November	1	0	1
2024 October	11	3	14
2024 September	10	5	15
2024 August	7	4	11
2024 July	13	2	15
2024 June	15	5	20
2024 May	18	4	22
2024 April	15	5	20
2024 March	19	2	21
2024 February	34	3	37
2024 January	15	6	21
2023 December	19	7	26
2023 November	11	4	15
2023 October	14	8	22
2023 September	12	0	12
2023 August	10	5	15
2023 July	6	2	8
2023 June	7	4	11
2023 May	24	3	27
2023 April	13	4	17
2023 March	25	3	28
2023 February	19	7	26
2023 January	9	5	14
2022 December	11	5	16
2022 November	14	6	20
2022 October	14	8	22
2022 September	17	17	34
2022 August	24	4	28
2022 July	16	5	21
2022 June	14	11	25
2022 May	9	9	18
2022 April	13	10	23
2022 March	17	8	25
2022 February	6	3	9
2022 January	32	8	40
2021 December	25	7	32
2021 November	10	10	20
2021 October	29	6	35
2021 September	16	12	28
2021 August	17	9	26
2021 July	15	3	18
2021 June	9	6	15
2021 May	11	4	15
2021 April	26	6	32
2021 March	13	8	21
2021 February	16	14	30
2021 January	10	9	19
2020 December	9	11	20
2020 November	13	3	16
2020 October	9	8	17
2020 September	12	9	21
2020 August	23	7	30
2020 July	13	0	13
2020 June	16	4	20
2020 May	15	3	18
2020 April	9	4	13
2020 March	21	7	28
2020 February	7	2	9
2020 January	15	1	16
2019 December	13	6	19
2019 November	14	9	23
2019 October	9	3	12
2019 September	14	6	20
2019 August	13	3	16
2019 July	20	7	27
2019 June	31	21	52
2019 May	54	60	114
2019 April	38	28	66
2019 March	8	3	11
2019 February	6	6	12
2019 January	7	7	14
2018 December	9	6	15
2018 November	8	1	9
2018 October	13	7	20
2018 September	6	4	10
2018 August	5	10	15
2018 July	11	3	14
2018 June	8	8	16
2018 May	7	8	15
2018 April	7	11	18
2018 March	8	2	10
2018 February	8	4	12
2018 January	7	0	7
2017 December	6	0	6
2017 November	11	0	11
2017 October	11	4	15
2017 September	7	9	16
2017 August	6	5	11
2017 July	9	5	14
2017 June	15	2	17
2017 May	16	4	20
2017 April	13	2	15
2017 March	9	59	68
2017 February	12	3	15
2017 January	13	0	13
2016 December	17	6	23
2016 November	22	2	24
2016 October	21	8	29
2016 September	31	5	36
2016 August	9	3	12
2016 July	7	0	7
2016 June	14	3	17
2016 May	11	19	30
2016 April	14	4	18
2016 March	19	9	28
2016 February	18	9	27
2016 January	10	7	17
2015 December	13	0	13
2015 November	12	4	16
2015 October	28	7	35
2015 September	14	8	22
2015 August	36	2	38
2015 July	27	6	33
2015 June	17	1	18
2015 May	17	4	21
2015 April	12	4	16
2015 March	2	0	2

Show all

Indexed in:

Follow us:

Statistics

Indexed in:

Follow us:

Statistics

Subscribe to our newsletter