def plotDry(self,kbri,koil,pref,ksat_well,gdry,phit,SO,gr,facies,vsh,kdry_well): plt.rcParams.update({'font.size': 12,'legend.fontsize': 12}) kdry = Gassman.kdry(ksat_well,phit,SO,vsh,kbri,koil) n1 = len(ksat_well) dry_ratio = np.zeros(n1) #for i in range(n1): #print(vsh[i]) fig = figure(4, figsize=(13, 4)) ax = plt.subplot(1,2,1) #ax.set_xlabel('porosity $\phi$') ax.set_ylabel('K dry (Pa)') #plot(self._phi[:],self._kdry_l[1,:],'b') #plot(self._phi[:],self._kdry_l[2,:],'b') plot(self._phi[:],self._kdry_l[pref,:],'b') plot(self._phi[:],self._kdry_c[:],'b') #plot(self._phi[:],self._kdryConstant[3,:],'b') plot(self._phi[:],self._kdryConstant[5,:],'b') plot(self._phi[:],self._kdryConstant[7,:],'b') plot(self._phi[:],self._kdryConstant[9,:],'b') #ax.set_xticklabels([]) ax.set_xlabel('porosity $\phi$') #plot(self._phi[:],self._kdryConstant[18,:],'b') ylim(1.e9,1.5e10) xlim(0,0.5) p=ax.scatter(phit,kdry_well,c=vsh,edgecolor='none') #p=ax.scatter(phit,kdry,c='r') p.set_clim([0,.5]) cbar= plt.colorbar(p) cbar.set_label('Vsh') #plt.rcParams.update({'font.size': 12,'legend.fontsize': 12}) bx = plt.subplot(1,2,2) bx.set_xlabel('porosity $\phi$') bx.set_ylabel('G dry (Pa)') plot(self._phi[:],self._gdry_l[pref,:],'b') plot(self._phi[:],self._gdry_c[:],'b') plot(self._phi[:],self._gdryConstant[5,:],'b') plot(self._phi[:],self._gdryConstant[7,:],'b') plot(self._phi[:],self._gdryConstant[9,:],'b') ylim(1e9,1.5e10) xlim(0,0.5) #bx.set_yticklabels([]) p=bx.scatter(phit,gdry,c=vsh,edgecolor='none') cbar= fig.colorbar(p) cbar.set_label('Vsh') p.set_clim([0,.5]) fig.savefig('model_calibration.png', bbox_inches='tight') fig = figure(5, figsize=(6, 6)) ax = plt.subplot(1,1,1) ax.set_xlabel('porosity $\phi$') ax.set_ylabel('K dry (Pa)') #plot(self._phi[:],self._kdry_l[1,:],'b') #plot(self._phi[:],self._kdry_l[2,:],'b') plot(self._phi[:],self._kdry_l[pref,:],'b') plot(self._phi[:],self._kdry_c[:],'b') #plot(self._phi[:],self._kdryConstant[2,:],'b') plot(self._phi[:],self._kdryConstant[5,:],'b') plot(self._phi[:],self._kdryConstant[7,:],'b') plot(self._phi[:],self._kdryConstant[9,:],'b') ylim(1.e9,4e10) xlim(0,0.5) fig.savefig('teory_model.pdf')
def plotDryratio(self,title,kbri,koil,ksat_well,gdry,kdry_well,phit,SO,gr,facies,depth,dmin,dmax,vsh): #phid,phit,ksat_well,vsh_iw,g,SO,depth,well,facies,sid = np.loadtxt(fileName,unpack=True) kdry= Gassman.kdry(ksat_well,phit,SO,vsh,kbri,koil) n1 = len(ksat_well) dry_ratio = np.zeros(n1) for i in range(n1): #print(vsh[i],depth[i]) dry_ratio[i]=(kdry_well[i]+4.0/3.0*gdry[i])/gdry[i] dry_ratio[i]=np.power(dry_ratio[i],0.5) fig = figure(7, figsize=(8, 6)) fig.suptitle(title) bx = plt.subplot(1,3,1) bx.set_xlabel('Gamma Ray') bx.set_ylabel('Depth (feet)') bx.xaxis.set_major_locator(MaxNLocator(4)) p=bx.plot(gr,depth,c='k') bx.set_ylim(dmax,dmin) #bx.set_title('') #bx.set_yticklabels([]) colors=[1.0,2.,3.,4.,5.,6.,7.,8.,9.,10.0,20.0,30.0,40.0,50.0,60.0,70.0,80.0,90.0] n = len(depth) for i in range(n): for j in range(len(colors)): if facies[i] == colors[j]: if colors[j] >=5 and colors[j]<=9: j = 4 c1='k' p=bx.scatter(gr[i],depth[i],c=c1) elif colors[j] >=50 and colors[j]<=90: j = 13 c1='k' p=bx.scatter(gr[i],depth[i],c=c1) else: c1=cm.hsv((j+1)/float(len(colors)),1) p=bx.scatter(gr[i],depth[i],c=c1,edgecolor='none') bx.set_xlim(0,160) bx = plt.subplot(1,3,2) bx.set_xlabel('Vp/Vs dry ratio') #bx.set_ylabel('Depth (ft)') #p=bx.scatter(dry_ratio,depth,c=SO,edgecolor='none') p1=bx.plot(dry_ratio,depth,c='k') bx.xaxis.set_major_locator(MaxNLocator(4)) #p.set_clim([0,0.5]) #cbar= plt.colorbar(p) #cbar.set_label('Oil saturation') ylim(dmax,dmin) bx.set_yticklabels([]) xlim(1.4,1.65) ''' bx = plt.subplot(1,5,3) bx.set_xlabel('Oil saturation') #bx.set_ylabel('Depth (ft)') #p=bx.scatter(dry_ratio,depth,c=SO,edgecolor='none') p1=bx.plot(SO,depth,c='k') bx.xaxis.set_major_locator(MaxNLocator(4)) #p.set_clim([0,0.5]) #cbar= plt.colorbar(p) #cbar.set_label('Oil saturation') bx.set_yticklabels([]) ylim(dmax,dmin) xlim(0,0.5) ''' kbri=2.8e9 koil=0.829e9 rhobri=1030 rhooil=0.66e3 #kgassman(kdry_well,phit,SO,vsh,kbri,koil,rhobri,rhooil) ksatR,rhosatR,soR,kfluid=Gassman.kgassman(kdry_well,phit,SO,vsh,kbri,koil,rhobri,rhooil) bx = plt.subplot(1,3,3) bx.set_xlabel('Kdry(GPa)') #bx.set_ylabel('Depth (ft)') #p=bx.scatter(dry_ratio,depth,c=SO,edgecolor='none') bx.plot(kdry_well/1e9,depth,c='k') bx.plot(kdry/1e9,depth,c='r',label='kdry obs') bx.set_xlim(0.0,17.0) #bx2 = bx.twiny() #bx2.plot(kfluid/1e9,depth,c='b',label='kfluid') #plt.legend() #bx2.set_xlabel('Kfluid(GPa)',color='b') bx.set_yticklabels([]) bx.xaxis.set_major_locator(MaxNLocator(4)) #p.set_clim([0,0.5]) #cbar= plt.colorbar(p) #cbar.set_label('Oil saturation') bx.set_ylim(dmax,dmin) #bx2.set_ylim(dmax,dmin) fig.savefig(title+'QC2.pdf') #bx2.set_xlim(0.0 ,2.5) ''' bx = plt.subplot(1,5,3) bx.set_xlabel('Oil saturation') #bx.set_ylabel('Depth (ft)') #p=bx.scatter(dry_ratio,depth,c=SO,edgecolor='none') #p1=bx.plot(SO,depth,c='k') p2=bx.plot(soR,depth,c='k') bx.xaxis.set_major_locator(MaxNLocator(4)) #p.set_clim([0,0.5]) #cbar= plt.colorbar(p) #cbar.set_label('Oil saturation') bx.set_yticklabels([]) ylim(dmax,dmin) xlim(0,0.5) ''' '''