Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Wrong name in petrophysics module #87

Open
efpeel opened this issue Aug 16, 2021 · 1 comment
Open

Wrong name in petrophysics module #87

efpeel opened this issue Aug 16, 2021 · 1 comment

Comments

@efpeel
Copy link

efpeel commented Aug 16, 2021

Dear all,

It was pointed out to me that the porosity_to_density function doesn't work. I traced it back to a from name in the module (phi_rho should be replaced by phi):

def porosity_to_density(phi, rho_matrix, rho_fluid):
    """
    Get density from a porosity log. Typical values:
        rho_matrix (sandstone) : 2650 kg/m^3
        rho_matrix (limestome): 2710 kg/m^3
        rho_matrix (dolomite): 2876 kg/m^3
        rho_matrix (anyhydrite): 2977 kg/m^3
        rho_matrix (salt): 20320 kg/m^3
        rho_fluid (fresh water): 1000 kg/m^3
        rho_fluid (salt water): 1100 kg/m^3
    See wiki.aapg.org/Density-neutron_log_porosity.
    Args:
        phi (ndarray): The porosity log.
        rho_matrix (float)
        rho_fluid (float)
    Returns:
        Estimate of bulk density, rho.
    """
    return rho_matrix * (1 - phi) + rho_fluid * phi_rhob
@kwinkunks
Copy link
Member

@EvanBianco If you get a chance, can you take a look at this?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants