A method for calculating subterranean mass distribution etc, from the gravity value at the surface is discussed for the three dimensional case. The principle of the method was already described in the preceding papers.
In the three dimensional case problems newly arise which are relate to the upper boundary of the spectrum of the gravity data. That is, different unit function or unit gravity must be used according to the assumed spectrum boundary.
In the present paper three kinds of unit functions are proposed, and these are summarized as followings.
1) When the upper limit of spatial wave number of the gravity data is independently given for x and y direction the unit function represented by
1)
is used.
2) When the upper limit of wave number is given in such a way as √ω
12+ω
22<Ω, the unit function represented by
is used,
2) where ω
1, ω
2 respectively represents wave number of x and y direction.
3) When the limit of wave number is limited in the hexagonal region on the wave number plane, the unit function represented by
is used.
The unit function δ_??_ is product form of trle unit function aiscussea in me two aimensional case. The gravity value given at lattice points is replaced by the smoothed function whose value concides with given value at each lattice points. The operations are weighted sum of that unit function. The mass distribution etc, are given by the weighted sum of the response for that unit gravity. A week point of this unit function is that the assumed limit of wave number varies according to the choice of the co-ordinate.
As it is natural that the wave number in the two dimensional spectrum is given by √ω
12 the unit function δ_??_ is reasonably used. This unit function is regorous but contains inconyenience, if we only consider simplicity of calculations. That is, the smoothed function which takes the place of the given data is given by the weighted integral of continuous data. And the mass distribution etc. are also given by the weighted integral of the response for that unit function. As a complomise of 1) and 2), it is desirable that both the upper limit of wave number is given as the case 2) and the results are obtained by the weighted sum of the data. In order to satisfy these conditions as much as possible, a hexagonal unit function δ_??_ is proposed. If this unit function is adopted, and if the data are given at the mesh points which are represented by the intersection points of the groups of straight lines having equal distance and each group are mutually inclined by 60°, the smoothed function and the necessary results are obtained as weighted sum of the unit function or of response for that unit function.
抄録全体を表示