anh_poro

The poro-elastic scaling applies a correction for poro-elasticity on top of the anharmonic scaling. The scaling method follows Appendix A of Takei, 2002, “Effect of pore geometry on VP/VS: From equilibrium geometry to crack”, JGR Solid Earth, DOI.

Requires

  • VBR.out.elastic.anharmonic result of anharmonic calculation.
  • VBR.in.SV.phi melt fraction or porosity (volume fraction)

Calling Procedure

% set required state variables
clear
VBR.in.SV.T_K = linspace(800,1200,10)+273; % temperature [K]
VBR.in.SV.P_GPa = 2 * ones(size(VBR.in.SV.T_K)); % pressure [GPa]
VBR.in.SV.rho = 3300 * ones(size(VBR.in.SV.T_K)); % density [kg m^-3]
VBR.in.SV.phi = 0.01 * ones(size(VBR.in.SV.T_K)); % melt fraction / porosity

% add to elastic methods list
VBR.in.elastic.methods_list={'anharmonic','anh_poro'};

% call VBR_spine
[VBR] = VBR_spine(VBR) ;

Output

Output is stored in VBR.out.elastic.anh_poro:

>> disp(fieldnames(VBR.out.elastic.anharmonic))

{
  [1,1] = Gu % unrelaxed shear modulus at desired P,T,phi
  [2,1] = Ku % unrelaxed bulk modulus at desired P,T,phi
  [3,1] = Vpu % unrelaxed compressional wave velocity
  [4,1] = Vsu % unrelaxed shear wave velocity
}

Examples

See Projects/vbr_core_examples/CB_009_anhporo.m:

'CB_009_anhporo'