Esempio n. 1
0
    def setUp(self):
        mesh = discretize.TensorMesh([30])
        sigma = np.ones(mesh.nC)
        model = np.log(sigma)

        # prob = DC.Simulation3DCellCentered(mesh, rhoMap=Maps.ExpMap(mesh))

        receivers = survey.BaseRx(20 * [[0.0]])
        source = survey.BaseSrc([receivers])
        sim = simulation.ExponentialSinusoidSimulation(
            mesh=mesh, survey=survey.BaseSurvey([source]), model_map=maps.ExpMap(mesh)
        )

        synthetic_data = sim.make_synthetic_data(model)
        dobs = synthetic_data.dobs

        self.relative = 0.01
        self.noise_floor = 1e-8

        synthetic_data.relative_error = self.relative
        synthetic_data.noise_floor = self.noise_floor

        dmis = data_misfit.L2DataMisfit(simulation=sim, data=synthetic_data)

        self.model = model
        self.mesh = mesh
        self.sim = sim
        self.survey = sim.survey
        # self.survey = survey
        # self.prob = prob
        self.data = synthetic_data
        self.dmis = dmis
Esempio n. 2
0
    def setUp(self):
        mesh = discretize.TensorMesh([100])
        self.sim = simulation.ExponentialSinusoidSimulation(
            mesh=mesh,
            model_map=maps.IdentityMap(mesh),
        )

        mtrue = np.zeros(mesh.nC)
        mtrue[mesh.vectorCCx > 0.3] = 1.0
        mtrue[mesh.vectorCCx > 0.45] = -0.5
        mtrue[mesh.vectorCCx > 0.6] = 0

        self.mtrue = mtrue
Esempio n. 3
0
    def setUp(self):
        mesh = discretize.TensorMesh([30])
        sigma = np.ones(mesh.nC)
        model = np.log(sigma)

        receivers = survey.BaseRx(20 * [[0.0]])
        source = survey.BaseSrc([receivers])

        self.sim = simulation.ExponentialSinusoidSimulation(
            mesh=mesh,
            survey=survey.BaseSurvey([source]),
            model_map=maps.ExpMap(mesh))

        self.dobs = self.sim.dpred(model)