• Title/Summary/Keyword: Hessian

Search Result 109, Processing Time 0.026 seconds

Construction the pseudo-Hessian matrix in Gauss-Newton Method and Seismic Waveform Inversion (Gauss-Newton 방법에서의 유사 Hessian 행렬의 구축과 이를 이용한 파형역산)

  • Ha, Tae-Young
    • Geophysics and Geophysical Exploration
    • /
    • v.7 no.3
    • /
    • pp.191-196
    • /
    • 2004
  • Seismic waveform inversion can be solved by using the classical Gauss-Newton method, which needs to construct the huge Hessian by the directly computed Jacobian. The property of Hessian mainly depends upon a source and receiver aperture, a velocity model, an illumination Bone and a frequency content of source wavelet. In this paper, we try to invert the Marmousi seismic data by controlling the huge Hessian appearing in the Gauss-Newton method. Wemake the two kinds of he approximate Hessian. One is the banded Hessian and the other is the approximate Hessian with automatic gain function. One is that the 1st updated velocity model from the banded Hessian is nearly the same of the result from the full approximate Hessian. The other is that the stability using the automatic gain function is more improved than that without automatic gain control.

Frequency domain elastic full waveform inversion using the new pseudo-Hessian matrix: elastic Marmousi-2 synthetic test (향상된 슈도-헤시안 행렬을 이용한 탄성파 완전 파형역산)

  • Choi, Yun-Seok;Shin, Chang-Soo;Min, Dong-Joo
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 2007.06a
    • /
    • pp.329-336
    • /
    • 2007
  • For scaling of the gradient of misfit function, we develop a new pseudo-Hessian matrix constructed by combining amplitude field and pseudo-Hessian matrix. Since pseudo- Hessian matrix neglects the calculation of the zero-lag auto-correlation of impulse responses in the approximate Hessian matrix, the pseudo-Hessian matrix has a limitation to scale the gradient of misfit function compared to the approximate Hessian matrix. To validate the new pseudo- Hessian matrix, we perform frequency-domain elastic full waveform inversion using this Hessian matrix. By synthetic experiments, we show that the new pseudo-Hessian matrix can give better convergence to the true model than the old one does. Furthermore, since the amplitude fields are intrinsically obtained in forward modeling procedure, we do not have to pay any extra cost to compute the new pseudo-Hessian. We think that the new pseudo-Hessian matrix can be used as an alternative of the approximate Hessian matrix of the Gauss-Newton method.

  • PDF

FROBENIUS ENDOMORPHISMS OF BINARY HESSIAN CURVES

  • Gyoyong Sohn
    • East Asian mathematical journal
    • /
    • v.39 no.5
    • /
    • pp.529-536
    • /
    • 2023
  • This paper introduces the Frobenius endomophisms on the binary Hessian curves. It provides an efficient and computable homomorphism for computing point multiplication on binary Hessian curves. As an application, it is possible to construct the GLV method combined with the Frobenius endomorphism to accelerate scalar multiplication over the curve.

Laplace-domain Waveform Inversion using the Pseudo-Hessian of the Logarithmic Objective Function and the Levenberg-Marquardt Algorithm (로그 목적함수의 유사 헤시안을 이용한 라플라스 영역 파형 역산과 레벤버그-마쿼트 알고리듬)

  • Ha, Wansoo
    • Geophysics and Geophysical Exploration
    • /
    • v.22 no.4
    • /
    • pp.195-201
    • /
    • 2019
  • The logarithmic objective function used in waveform inversion minimizes the logarithmic differences between the observed and modeled data. Laplace-domain waveform inversions usually adopt the logarithmic objective function and the diagonal elements of the pseudo-Hessian for optimization. In this case, we apply the Levenberg-Marquardt algorithm to prevent the diagonal elements of the pseudo-Hessian from being zero or near-zero values. In this study, we analyzed the diagonal elements of the pseudo-Hessian of the logarithmic objective function and showed that there is no zero or near-zero value in the diagonal elements of the pseudo-Hessian for acoustic waveform inversion in the Laplace domain. Accordingly, we do not need to apply the Levenberg-Marquardt algorithm when we regularize the gradient direction using the pseudo-Hessian of the logarithmic objective function. Numerical examples using synthetic and field datasets demonstrate that we can obtain inversion results without applying the Levenberg-Marquardt method.

Computer Aided Optimal Circuit Design (전자계산기에 의한 최적회로설계 방식 연구)

  • 김덕진;김선영
    • Journal of the Korean Institute of Telematics and Electronics
    • /
    • v.14 no.4
    • /
    • pp.22-31
    • /
    • 1977
  • A general equation by which the Hessian matrix of an error function can be determined directly, has been derived. It was verified to be useful in optimization processes that include the Hessian matrix. A few design examples had shown that this method had accelerated the processes of finding the minimums. The advantage of this technique is the possibility of optimizing functions that composed of both the phases and magnitudes.

  • PDF

Pseudo-multiscale Waveform Inversion for Velocity Modeling

  • Yang Dongwoo;Shin Changsoo;Yoon Kwangjin;Yang Seungjin;Suh Junghee;Hong Soonduk
    • Proceedings of the KSEEG Conference
    • /
    • 2002.04a
    • /
    • pp.159-162
    • /
    • 2002
  • We tried to obtain an initial velocity model for prestack depth migration via waveform inversion. For application of any field data we chose a smooth background layered velocity model (v=v0 + k x z) as an initial velocity model. Newton type waveform inversion needs to invert huge Hessian matrix. In order to compute full Hessian matrix arising from full aperture data and full illumination zone, we meet insurmountable difficulties of paying astronomical computing cost. For the layered media, approximate Hessian emerging from single shot aperture data can be used repeatedly for split spread source configuration. In our work of using this Hessian characteristic of layered media we attempted to obtain the approximate velocity model as close as possible to the true velocity model in first iteration.

  • PDF

F-Hessian SIFT-Based Railroad Level-Crossing Vision System (F-Hessian SIFT기반의 철도건널목 영상 감시 시스템)

  • Lim, Hyung-Sup;Yoon, Hak-Sun;Kim, Chel-Huan;Ryu, Deung-Ryeol;Cho, Hwang;Lee, Key-Seo
    • The Journal of the Korea institute of electronic communication sciences
    • /
    • v.5 no.2
    • /
    • pp.138-144
    • /
    • 2010
  • This paper presents the experimental analysis of a F-Hessian SIFT-Based Railroad Level-Crossing Safety Vision System. Region of surveillance, region of interests, data matching based on extracting feature points has been examined under the laboratory condition by the model rig on a small scale. Real-time system were observed by using SIFT based on F-Hessian feature tracking method and other common algorithm.

Rock Fracture Centerline Extraction based on Hessian Matrix and Steger algorithm

  • Wang, Weixing;Liang, Yanjie
    • KSII Transactions on Internet and Information Systems (TIIS)
    • /
    • v.9 no.12
    • /
    • pp.5073-5086
    • /
    • 2015
  • The rock fracture detection by image analysis is significant for fracture measurement and assessment engineering. The paper proposes a novel image segmentation algorithm for the centerline tracing of a rock fracture based on Hessian Matrix at Multi-scales and Steger algorithm. A traditional fracture detection method, which does edge detection first, then makes image binarization, and finally performs noise removal and fracture gap linking, is difficult for images of rough rock surfaces. To overcome the problem, the new algorithm extracts the centerlines directly from a gray level image. It includes three steps: (1) Hessian Matrix and Frangi filter are adopted to enhance the curvilinear structures, then after image binarization, the spurious-fractures and noise are removed by synthesizing the area, circularity and rectangularity; (2) On the binary image, Steger algorithm is used to detect fracture centerline points, then the centerline points or segments are linked according to the gap distance and the angle differences; and (3) Based on the above centerline detection roughly, the centerline points are searched in the original image in a local window along the direction perpendicular to the normal of the centerline, then these points are linked. A number of rock fracture images have been tested, and the testing results show that compared to other traditional algorithms, the proposed algorithm can extract rock fracture centerlines accurately.