Commit ·
079252b
1
Parent(s): 9d81a04
v21: revert Rwa correction (hurt score), keep RMFT Fahrenheit fix
Browse files- pipeline.py +0 -21
pipeline.py
CHANGED
|
@@ -255,26 +255,6 @@ def get_rw(las, gr=None):
|
|
| 255 |
return estimate_rw(las, gr=gr)
|
| 256 |
|
| 257 |
|
| 258 |
-
def rwa_correct_rw(rw_sp, phie, rt, vsh):
|
| 259 |
-
"""
|
| 260 |
-
Cross-check Rw via Rwa method. SP can be suppressed in thin beds,
|
| 261 |
-
inflating Rw and causing missed pay. If Rwa is significantly lower than
|
| 262 |
-
SP_Rw, use geometric mean as a more conservative estimate.
|
| 263 |
-
"""
|
| 264 |
-
if rt is None or phie is None or vsh is None:
|
| 265 |
-
return rw_sp
|
| 266 |
-
clean = (vsh < 0.40) & (phie > 0.08) & np.isfinite(rt) & (rt > 0.05) & (rt < 200)
|
| 267 |
-
if clean.sum() < 5:
|
| 268 |
-
return rw_sp
|
| 269 |
-
rwa = rt[clean] * phie[clean]**2
|
| 270 |
-
rwa_ok = rwa[(rwa > 0.003) & np.isfinite(rwa)]
|
| 271 |
-
if len(rwa_ok) < 5:
|
| 272 |
-
return rw_sp
|
| 273 |
-
rw_rwa = float(np.nanpercentile(rwa_ok, 10))
|
| 274 |
-
# Only correct when Rwa is strongly lower (SP suppressed by >3x)
|
| 275 |
-
if rw_rwa < rw_sp / 3.0:
|
| 276 |
-
return float(np.clip(np.sqrt(rw_sp * rw_rwa), 0.005, 2.0))
|
| 277 |
-
return rw_sp
|
| 278 |
|
| 279 |
|
| 280 |
def compute_sw(phie, rt, rw):
|
|
@@ -539,7 +519,6 @@ def process_well(input_path, output_path):
|
|
| 539 |
phie = compute_phie(phit, vsh)
|
| 540 |
|
| 541 |
rw = get_rw(las, gr=gr)
|
| 542 |
-
rw = rwa_correct_rw(rw, phie, rt, vsh)
|
| 543 |
if rt is not None:
|
| 544 |
sw = compute_sw(phie, rt, rw)
|
| 545 |
else:
|
|
|
|
| 255 |
return estimate_rw(las, gr=gr)
|
| 256 |
|
| 257 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 258 |
|
| 259 |
|
| 260 |
def compute_sw(phie, rt, rw):
|
|
|
|
| 519 |
phie = compute_phie(phit, vsh)
|
| 520 |
|
| 521 |
rw = get_rw(las, gr=gr)
|
|
|
|
| 522 |
if rt is not None:
|
| 523 |
sw = compute_sw(phie, rt, rw)
|
| 524 |
else:
|