Пример #1
0
    def test_akirichards_alt(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.

        reflect = avo.akirichards_alt(vp1, vs1, rho1, vp2, vs2, rho2, theta)

        self.assertAlmostEquals(reflect, -0.112236, places=3)
Пример #2
0
    def test_akirichards_alt(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = 40.


        reflect = avo.akirichards_alt( vp1,vs1,rho1,vp2,
                                        vs2,rho2,theta )

        self.assertAlmostEquals( reflect, -0.112236, places = 3 )
Пример #3
0
    def test_akirichards_alt(self):

        vp1 = 12250.
        vp2 = 11600.

        vs1 = 6200.
        vs2 = 6650.

        rho1 = 2.66
        rho2 = 2.34

        theta = np.arange(45)

        reflect = avo.akirichards_alt(vp1, vs1, rho1, vp2, vs2, rho2, theta)
        reflect_zoep = avo.zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)

        # See if it is within .1 of zoep for < 45 deg
        test = np.allclose(reflect, reflect_zoep, rtol=self.tolerance)
        self.assertTrue(test)
Пример #4
0
def akirichards_alt(Rp0, Rp1, theta1):
    reflection.akirichards_alt(Rp0.vp, Rp0.vs, Rp0.rho, Rp1.vp, Rp1.vs, Rp1.rho, theta1)