Assessment of rock formation permeability is a complicated and challenging problem that plays a key role in oil reservoir modeling, production forecast, and the optimal exploitation management. Generally, permeability evaluation is performed using porosity-permeability relationships obtained by integrated analysis of various petrophysical measurements taken from cores and wireline well logs. Dependence relationships between pairs of petrophysical variables, such as permeability and porosity, are usually nonlinear and complex, and therefore those statistical tools that rely on assumptions of linearity and/or normality and/or existence of moments are commonly not suitable in this case. But even expecting a single copula family to be able to model a complex bivariate dependency seems to be still too restrictive, at least for the petrophysical variables under consideration in this work. Therefore, we explore the use of the Bernstein copula, and we also look for an appropriate partition of the data into subsets for which the dependence strucure was simpler to model, and then a conditional gluing copula technique is applied to build the bivariate joint distribution for the whole data set.