import m8r # top of reservoir top = m8r.math(n1=301,d1=1,o1=0,output=0.06,label1='Trace')[0] # bottom of reservoir bot = m8r.math(output='input*(2-x1/300)')[top] # wedge model wedge=m8r.cat(axis=2).unif2(n1=181,o1=0,d1=0.001,label1='Time',unit1='s',v00=[10,15,20])[top,bot] # seismic convolution modeling seis = m8r.ai2refl.ricker1(frequency=25)[wedge] # plotting seis.grey(title='Wedge Model',color='G').show() m8r.window(j2=10)[seis].wiggle(title="Wedge Model",poly=True,yreverse=True,transp=True).show()
import m8r
import matplotlib.pyplot as plt
model = m8r.math(n1=201,o1=-6.283185,d1=.06283185,
n2=201,o2=-6.283185,d2=.06283185,
n3=101,o3=-3.141592,d3=.06283185,
output='sin(x1*x3) * cos(x2)*sin(x3) + sin(0.5*x3)')[0]
model.grey(color='j',gainpanel='a',title='RSF').show()
plt.imshow(model[0,:,:])
plt.show()
import m8r
import matplotlib.pyplot as plt
model = m8r.math(n1=201,
o1=-6.283185,
d1=.06283185,
n2=201,
o2=-6.283185,
d2=.06283185,
n3=101,
o3=-3.141592,
d3=.06283185,
output='sin(x1*x3) * cos(x2)*sin(x3) + sin(0.5*x3)')[0]
model.grey(color='j', gainpanel='a', title='RSF').show()
plt.imshow(model[0, :, :])
plt.show()
import m8r
# top of reservoir
top = m8r.math(n1=301, d1=1, o1=0, output=0.06, label1='Trace')[0]
# bottom of reservoir
bot = m8r.math(output='input*(2-x1/300)')[top]
# wedge model
wedge = m8r.cat(axis=2).unif2(n1=181,
o1=0,
d1=0.001,
label1='Time',
unit1='s',
v00=[10, 15, 20])[top, bot]
# seismic convolution modeling
seis = m8r.ai2refl.ricker1(frequency=25)[wedge]
# plotting
seis.grey(title='Wedge Model', color='G').show()
m8r.window(j2=10)[seis].wiggle(title="Wedge Model",
poly=True,
yreverse=True,
transp=True).show()