Accept a COMPVD phase change away from the type-1 gas-oil contact

[GitPaean]

With EQUIL item 10 = 1 the COMPVD rows give the composition versus depth. A two-zone table whose last vapour row and first liquid row do not straddle the gas-oil contact stopped the run. The gas zone now ends at the nearest edge of the gap between those rows, with a warning, and a cell on the last vapour row stays in the gas zone. Type 3 keeps the error, as it uses the contact as reference depth.

When the gas zone reaches the water-oil contact, the water now takes its pressure from the gas column there rather than from the liquid one.

test_compequil covers these cases; the existing equilibration regressions are unchanged. The new regression equil_1d_compvd_water_gascap_contact_mismatch runs the case from OPM/opm-tests#1638 and needs its reference data.

Initialization flow

EQUIL supplies contact depths and reference pressure. COMPVD supplies composition versus depth, with each row marked as vapour or liquid. The following flow applies to a COMPVD table containing both vapour and liquid rows.

flowchart TD
    A["COMPVD contains both vapour and liquid rows"] --> B{"Are all vapour rows<br/>above the liquid rows?"}

    B -->|No| C["Stop: interleaved phase rows"]
    B -->|Yes| D{"Is the EQUIL gas-oil contact<br/>between the last vapour row<br/>and first liquid row?"}

    D -->|Yes| E["Use the EQUIL contact<br/>as the zone boundary"]
    D -->|No| F{"EQUIL initialization type?"}

    F -->|Type 1| G["Move the boundary to the<br/>nearest edge of the row gap<br/>and issue a warning"]
    F -->|Type 3| H["Stop: type 3 requires the<br/>contact as its reference depth"]

    E --> I["Use the selected boundary for<br/>gas and liquid pressure integration"]
    G --> I

    I --> J["Initialize water pressure when needed"]
    J --> K["Assign each cell its zone composition<br/>and nominal phase"]
    K --> L["For type 1, a cell exactly on the<br/>last vapour row remains gas"]
    L --> M["Downstream flash computes<br/>the final hydrocarbon phase split"]

Loading

For example, with the last vapour row at 2049 m and first liquid row at 2051 m, type 1 moves a contact at 2060 m to 2051 m, or 2030 m to 2049 m. Depth increases downward.

Water-pressure anchoring

For a two-zone COMPVD region, when water pressure needs anchoring from the hydrocarbon column, the phase just above the water contact determines which pressure to use.

flowchart TD
    A["Water pressure needs anchoring<br/>from the hydrocarbon column"] --> B{"Does the gas zone reach<br/>the water contact?"}
    B -->|Yes| C["Take gas pressure<br/>at the water contact"]
    B -->|No| D["Take liquid pressure<br/>at the water contact"]
    C --> E["Subtract the specified capillary pressure<br/>and integrate the water column"]
    D --> E

Loading

This prevents a pressure discontinuity when moving the gas/liquid boundary makes gas touch water.

[GitPaean]

jenkins build this opm-tests=1638 failure_report please

https://ci.opm-project.org/job/opm-simulators-PR-builder/10822/

[GitPaean]

jenkins build this opm-tests=1638 failure_report please