100 REM ****************************************************************** 110 COLOR 15:LOCATE 9,15 200 PRINT" NEW SHALY SAND EQUATION FOR RESERVOIRS WITH DISPERSED CLAY" 210 LOCATE 11,15 300 PRINT" by G. B. ASQUITH 310 LOCATE 20,15:PRINT " Press any key to continue" 320 Q$=INKEY$: IF Q$="" THEN 320 330 CLS 350 REM GEOBYTE - JUNE 1990 - pp.43-50 400 REM THIS SHALY SAND TECHNIQUE DOES NOT REQUIRE VALUES FOR VOLUME OF CLAY 500 REM TOTAL SHALE POROSITY,OR ADJACENT SHALE RESISTIVITY. THE TECHNIQUE 600 REM USES DENSITY POROSITY FOR EFFECTIVE POROSITY AND RESISTIVITY POROSITY 700 REM FOR TOTAL POROSITY. THE PERCENT WATER FILLED CLAY BOUND MICROPOROSITY 800 REM IS CALCULATED BY SUBSTRACTING PHIe FROM PHIt. EFFECTIVE WATER 900 REM SATURATION(Swe) IS DETERMINED BY FIRST CALCULATING EFFECTIVE BULK 1000 REM VOLUME WATER(BVWe=BVWt-% WATER FILLED CLAY BOUND MICROPOROSITY) 1100 REM AND THEN BY DIVIDING BVWe BY PHIe. 1200 REM THE ADVANTAGE OF THIS METHOD IS THAT IT DOES NOT REQUIRE VALUES FOR 1300 REM Vcl,PHItsh, AND Rsh AS DOES THE DUAL WATER METHOD, OR CEC VALUES 1400 REM NEEDED BY THE WAXMAN-SMIT METHOD. 1500 REM ***************************************************************** 1600 DIM SW(50),RHOB(50),PHIT(50),RT(50),RXO(50),PHIR(50),PHIE(50) 1700 DIM SWE(50),BVWT(50),BVWE(50),SXO(50) 1800 DIM A$(50),B$(50),C$(50),D$(50),F$(50),DEPTH(50) 1900 LET Y$ = "Y" 2000 LET YY$="y" 2100 LET N$="N" 2200 LET NN$="n" 2300 CLS 2350 COLOR 15:LOCATE 2,14: PRINT "EFFECTIVE WATER SATURATION IN SHALY SANDS" 2400 COLOR 12:LOCATE 5,10: PRINT "INPUT NUMBER OF DEPTHS INPUT ZERO TO STOP":LOCATE 5,55 :COLOR 11:INPUT "";N 2500 IF N=0 THEN 22900 ELSE 2600 2600 CLS 2700 PRINT "INPUT WELL NAME"; 2800 INPUT A$ 2900 CLS 3000 PRINT "INPUT OPERATOR"; 3100 INPUT B$ 3200 CLS 3300 PRINT "INPUT LOCATION"; 3400 INPUT C$ 3500 CLS 3600 PRINT "INPUT COUNTY AND STATE"; 3700 INPUT D$ 3800 CLS 3900 PRINT "INPUT FORMATION NAME"; 4000 INPUT F$ 4100 CLS 4200 COLOR 10:PRINT"INPUT BHT, TD, & SURFACE TEMPERATURE"; 4300 INPUT BHT,TD,TS 4400 CLS 4500 PRINT"INPUT Rmf TEMPERATURE"; 4600 INPUT RMFT 4700 CLS 4800 COLOR 11:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++" 4900 COLOR 15:PRINT" BHT TD Ts Rmf TEMP" 5000 COLOR 14:PRINT BHT,TD,TS,RMFT 5100 COLOR 11:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++" 5200 PRINT 5300 COLOR 12:PRINT"ARE THESE CORRECT"; 5400 INPUT T$ 5500 IF T$=N$ OR T$=NN$ THEN 4100 5600 CLS 5700 FOR X=1 TO N 5800 COLOR 11:PRINT "INPUT DEPTH and Rt"; 5900 INPUT DEPTH(X),RT(X) 6000 PRINT 6100 NEXT X 6200 CLS 6300 PRINT 6400 FOR X=1 TO N 6500 PRINT "DEPTH=";DEPTH(X) 6600 PRINT "INPUT Rhob and Rxo"; 6700 INPUT RHOB(X),RXO(X) 6800 PRINT 6900 NEXT X 7000 CLS 7100 COLOR 12:PRINT " - LOG DATA - 7200 COLOR 14:PRINT"##############################################################" 7300 COLOR 11:PRINT" DEPTH Rt Rxo Rhob" 7400 FOR X=1 TO N 7500 COLOR 15:PRINT DEPTH(X) ,RT(X) ,RXO(X) ,RHOB(X) 7600 NEXT X 7700 COLOR 14:PRINT"##############################################################" 7800 PRINT 7900 COLOR 10:PRINT "ARE THESE DATA CORRECT"; 8000 INPUT E$ 8100 IF E$=N$ OR E$=NN$ THEN 5600 8200 CLS 8300 COLOR 11:PRINT"IS YOUR SHALLOW RESISTIVITY LOG AN RXO DEVICE"; 8400 INPUT G$ 8500 IF G$=N$ OR G$=NN$ THEN 13200 8600 CLS 8700 PRINT"INPUT Rmf & Rw"; 8800 INPUT RMF,RW 8900 CLS 9000 COLOR 10:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" 9100 COLOR 11:PRINT" Rmf Rw" 9200 COLOR 15:PRINT RMF,RW 9300 COLOR 10:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" 9400 PRINT 9500 COLOR 14:PRINT"ARE THESE CORRECT"; 9600 INPUT T$ 9700 IF T$=N$ OR T$=NN$ THEN 8600 9800 LET GRAD=(BHT-TS)/TD 9900 LET TF=(GRAD*DEPTH(1))+TS 10000 LET RMF=(RMF*(RMFT+6.77))/(TF+6.77) 10100 CLS 10200 COLOR 11:PRINT"IN ORDER TO CALCULATE RESISTIVITY DERIVED POROSITY A VALUE FOR" 10300 PRINT"Sxo MUST BE SELECTED" 10400 PRINT"*****************************************************************" 10500 PRINT" API Sxo" 10600 PRINT"GAS 60-90" 10700 PRINT"HIGH GRAVITY OIL 40-50 90-95" 10800 PRINT"MEDIUM GRAVITY OIL 20-40 80-90" 10900 PRINT"LOW GRAVITY OIL 10-20 70-80" 11000 PRINT"NO Sxo CORRECTION 100" 11100 PRINT"FOR VARIABLE Sxo(Sxo=Sw^0.2) INPUT ZERO" 11200 PRINT"***************************************************************** 11300 COLOR 15:PRINT"SELECT A VALUE FOR Sxo"; 11400 INPUT SXO(1) 11500 LET SXO(1)=SXO(1)/100 11600 IF SXO(1)=O THEN 12200 11700 FOR X=1 TO N 11800 LET SXO(X)=SXO(1) 11900 NEXT X 12000 CLS 12100 GOTO 12800 12200 FOR X= 1 TO N 12300 LET PHIE(X)=(2.65-RHOB(X))/(2.65-1) 12400 LET SXO(X)=((1/(PHIE(X)^2))*(RW/RT(X)))^.5 12500 LET SXO(X)=(SXO(X))^.2 12600 NEXT X 12700 CLS 12800 FOR X=1 TO N 12900 LET PHIT(X)=(((RMF/RXO(X)))/(SXO(X)^2))^.5 13000 NEXT X 13100 GOTO 19100 13200 CLS 13300 PRINT"INPUT Rmf & Rw"; 13400 INPUT RMF,RW 13500 CLS 13600 COLOR 10:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" 13700 COLOR 11:PRINT" Rmf RW" 13800 COLOR 15:PRINT RMF,RW 13900 COLOR 10:PRINT"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" 14000 PRINT 14100 COLOR 14:PRINT"ARE THESE CORRECT"; 14200 INPUT T$ 14300 IF T$=N$ OR T$=NN$ THEN 13200 14400 CLS 14500 LET GRAD=(BHT-TS)/TD 14600 LET TF=(GRAD*DEPTH(1))+TS 14700 LET RMF=(RMF*(RMFT+6.77))/(TF+6.77) 14800 IF G$=Y$ OR G$=YY$ THEN 10100 14900 CLS 15000 COLOR 10:PRINT"IN ORDER TO CALCULATE MIXED WATER RESISTIVITY(Rz) A VALUE FOR z" 15100 PRINT"MUST BE SELECTED 15200 PRINT 15300 PRINT"*****************************************************************" 15400 PRINT"FINE GRAINED SS, LOW PERMEABILITY z=0.10" 15500 PRINT"AVERAGE SS z=0.075" 15600 PRINT"COARSE GRAINED SS, HIGH PERMEABILITY z=0.05" 15700 PRINT"*****************************************************************" 15800 COLOR 14:PRINT"SELECT A VALUE"; 15900 INPUT Z 16000 CLS 16100 LET RZ=1/((Z/RW)+((l-Z)/RMF)) 16200 COLOR 11:PRINT"IN ORDER TO CALCULATE RESISTIVITY DERIVED POROSITY A VALUE FOR" 16300 PRINT"SXo MUST BE SELECTED" 16400 PRINT 16500 PRINT"*****************************************************************" 16600 PRINT" API Sxo" 16700 PRINT"GAS 60-90" 16800 PRINT"HIGH GRAVITY OIL 40-50 90-95" 16900 PRINT"MEDIUM GRAVITY OIL 20-40 80-90" 17000 PRINT"LOW GRAVITY OIL 10-20 70-80" 17100 PRINT"NO Sxo CORRECTION 100" 17200 PRINT"FOR VARIABLE Sxo(Sxo=Sw^0.2) INPUT ZERO" 17300 PRINT"*****************************************************************" 17400 COLOR 15:PRINT"SELECT A VALUE FOR SXO"; 17500 INPUT SXO(1) 17600 LET SXO(1)=SXO(1)/100 17700 IF SXO(1)=O THEN 18200 17800 FOR X=1 TO N 17900 LET SXO(X)=SXO(1) 18000 NEXT X 18100 GOTO 18700 18200 FOR X=1 TO N 18300 LET PHIE(X)=(2.65-RHOB(X))/(2.65-1) 18400 LET SXO(X)=((1/(PHIE(X)^2))*(RW/RT(X)))^.5 18500 LET SXO(X)=(SXO(X))^.2 18600 NEXT X 18700 FOR X=1 TO N 18800 LET PHIT(X)=(((RZ/RXO(X)))/(SXO(X)^2))^.5 18900 NEXT X 19000 CLS 19100 FOR X=1 TO N 19200 LET PHIE(X)=(2.65-RHOB(X))/(2.65-1!) 19300 LET SW(X)=((l/PHIE(X)^2)*(RW/RT(X)))^.5 19400 LET BVWT(X)=SW(X)*PHIE(X) 19500 LET BVWE(X)=BVWT(X)-(PHIT(X)-PHIE(X)) 19600 LET SWE(X)=BVWE(X)/PHIE(X) 19700 NEXT X 19800 FOR X=l TO N 19900 LET PHIT(X)=INT(PHIT(X)*1000+.5)/1000 20000 LET PHIE(X)=INT(PHIE(X)*1000+.5)/1000 20100 LET SW(X)=INT(SW(X)*1000+.5)/1000 20200 LET BVWT(X)=INT(BVWT(X)*1000+.5)/1000 20300 LET BVWE(X)=INT(BVWE(X)*1000+.5)/1000 20400 LET SWE(X)=INT(SWE(X)*1000+.5)/1000 20500 LET SXO(X)=INT(SXO(X)*1000+.5)/1000 20600 NEXT X 20700 PRINT"WELL NAME: ";A$ 20800 PRINT"OPERATOR: ";B$ 20900 PRINT"LOCATION: ";C$ 21000 PRINT"COUNTY & STATE: ";D$ 21100 PRINT"FORMATION: ";F$ 21200 LOCATE 6,1:COLOR 10:PRINT USING"Rmf=#.###";RMF:LOCATE 6,12:COLOR 14:PRINT USING "Rz=#.###";RZ:LOCATE 6,23:COLOR 13:PRINT USING "Tf=###.#";TF 21300 COLOR 11:PRINT USING"Rw=#.###";RW 21400 COLOR 12:PRINT USING"z=#.###";Z 21500 PRINT 21600 COLOR 14:PRINT"##################################################################" 21700 COLOR 10:PRINT" DEPTH PHIt Swt PHIe Swe Sxo " 21800 FOR X=1 TO N 21900 COLOR 11:PRINT DEPTH(X);PHIT(X),SW(X),PHIE(X),SWE(X),SXO(X) 22000 NEXT X 22100 COLOR 14:PRINT"###########################X######################################" 22200 COLOR 12:PRINT"WARNING:IF Rw APPROXIMATELY EQUAL TO Rmf Swe UNRELIABLE" 22300 PRINT 22400 COLOR 15:PRINT"DO YOU WANT TO CALCULATE Swe WITH DIFFERENT Sxo OR z VALUES"; 22500 INPUT X$ 22600 IF X$=Y$ OR X$=YY$ THEN 14800 22700 CLS 22800 GOTO 1900 22900 END