The PageRank method for automatic detection of microseismic events Huiyu Zhu*, Jie Zhang, University of Science and Technology of China (USTC)

Transcription

1 ownloaded /5/ to Redistribution subject to SEG license or copyright; see Terms of Use at The PageRank method for automatic detection of microseismic events Huiyu Zhu*, Jie Zhang, University of Science and Technology of hina (UST) Summary In microseismic monitoring, the fundamental task is to detect events automatically. The majority of microseismic events may suffer the limits from low signal-to-noise ratio for detection. Using a single master event to correlate the entire dataset may not work because events do not necessarily all exhibit sufficient similarity to the same event. uring hydraulic fracturing, however, it is likely that event may be similar to event, and event may be similar to event, but event may not be similar to event on the same standard. The situation of this event connectivity is similar to web pages, in which page content similarity is detected by PageRank. PageRank is one of the initial web page searching algorithms implemented by Google to estimate the links between documents. pplying that technology to microseismic problems, we rank windows of recorded data for similarity measurement, thus detect which windows likely contain potential events. We generate synthetic microseismic events with different source mechanisms, and demonstrate that the approach can handle noise and variations in waveform data reasonably well. Introduction Low-permeability oil reservoir and gas shales are difficult to produce. Hydraulic fracturing technology is often applied to connect the pathway and force hydrocarbons flow out. The process of hydrofracturing may induce microseismic events, which may be monitored by a surface or downhole seismic array. uring data processing, we build a velocity model and locate the microseismic events either by traveltimes or waveforms. These event locations may suggest where fracturing takes place. mong these efforts, however, the very first step is to detect hundreds to thousands of microseismic events from a large and noisy dataset. Most microseismic events are small and signal-to-noise ratio is low, thus, automatic detection is difficult. One may use master events to correlate data traces and detect weak events. Previous studies have also shown that the correlation detection can be effective as long as the separation between the master event and the target event is less than the dominant wavelength, and the master event and target events are similar (Gibbons and Ringdal, 6). In real problems, sometimes using a single master event to correlate the entire dataset may not work. We often face with a problem that event may be similar to event, and event may be similar to event, but event may not be similar to event on the same standard. This is due to source focal mechanisms associated with fracturing. The problem is similar to mining important pages among massive web pages, and it is solved by applying PageRank method to solve a connectivity matrix problem (Page et al., 999). Motivated by guiar and eroza () applying PageRank to detect weak tremors during a large earthquake, we explore the method for detecting weak microseismic events. The PageRank method can be applied to rank windows of recorded data according to mutual links. page will have a high rank if the sum of the ranks of its backlinks is high. Figure shows a possible relationship among four events. 6 x x x x x Figure : Event relationship: window is linked (similar) to window, window is linked to window, window is linked to window, but window is indirectly linked to window (since -> -> -> ). This relationship could be due to source focal mechanism associated with consistent fracturing pattern and variations of the pattern. The PageRank method The PageRank is the probability of surfing a certain web page under random surfing model. In the initial condition, the probability of each page is the same value. Like Markov hain, PageRank is computed from the previous stage and it calculates the value iteratively. SEG SEG enver nnual Meeting OI Page 6

2 ownloaded /5/ to Redistribution subject to SEG license or copyright; see Terms of Use at Let n be the number of pages. To describe their relationship on links, we generate a connective matrix G. If there is a hyperlink from page i to page j, then gij=, otherwise gij=. Then G matrix can be very huge, but very sparse. Let cj be the column sum of G: c j g () Then, let p be the probability that the random walk follows a link. typical value for p is.85 and -p is the probability that some arbitrary page is chosen. So we can use δto represent a particular page that is randomly chosen, see the Equation. ( p) n () We generate matrix, whose size is n by n. The elements in are p gij cj : cj aij () n: cj The matrix is the transition probability matrix of the Markov hain. The elements in are all strictly between zero and one, and the column sums are all equal to one. We can apply Perron-Frobenius theorem, an important result in matrix theory, to such matrix. It concludes that a nonzero solution of the Equation exists. x x () nd it is unique to within a scaling factor. If this scaling factor is chosen so that xi (5) i Then, x is the state vector of the Markov hain and is the PageRank value. We could solve the same matrix problem by replacing page with microseismic events in the above method. pplication and examples We apply PageRank method in microseismic event detection. In synthetic test, we design a six-layer model (Table ). Thickness(km) Vs Vp/Vs ensity(g/cm ) (km/s) Table : Model layer parameters. i ij We apply elastic wave modeling of a point earthquake source in a multi-layered half space using the Thompson- Haskell propagator matrix technique (Zhu and Rivera, ) to calculate synthetic waveforms with different source focal mechanisms. We assume receivers on the surface. We calculate eight events in total with magnitudes varying from - to. Table shows the event information. epth magnitude Strike ip Rake (km) Table : There are eight synthetic events. The strike of all the events is assumed exactly same but with dip and rake varying. We divide the data into the same length of windows. Each window is lagged by one sample point. Then we correlate each window with all other windows and calculate the correlation coefficient () values of paired windows. uring the calculation, we find the population of the correlation coefficient () values following a normal distribution approximately as shown in Figure. We can establish a threshold of detection according to the population. normalized histogram distribution of values theoretical values for all window pairs Figure : The blue bars are the distribution of orrelation oefficient () values for window pairs. The red line is the theoretical normal distribution plot. The distribution of values allows norm distribution. SEG SEG enver nnual Meeting OI Page 6

3 ownloaded /5/ to Redistribution subject to SEG license or copyright; see Terms of Use at We establish our detection threshold on the basis of Gaussian distribution with zero mean (where the mean is -., very close to zero, so we regard it as zero). We focus on the large positive values to declare a positive detection. For a normal distribution, we can use σ or σ as a measure to establish a threshold of detection. The positive threshold of σ corresponds to a two-sided significance level of 65%, while σ corresponds to a twosided significance level of 99.7%. ut there is a tradeoff between threshold and positive detection. s if we choose a higher threshold, it will provide more confident matches, but sacrifice some positive correlation results for the low- SNR data. If we give a lower threshold, it will give more positive correlation results, but there should be less confidence on the matches. Figure shows a comparison between two different levels of threshold. x x - Page Rank threshold σ 5 x - Page Rank threshold σ The Signal-to-noise ratio is Figure : : synthetic data without noise. : synthetic data with noise, and its SNR is and are corresponding to the PageRank values for σ and σ as a threshold. Figure gives us ideas on choosing the threshold value. In our test, we prefer to use σ to be the threshold, which should provide more potential matches. 5 x x - Page Rank The Signal-to-noise ratio is Figure : The result for low signal-to-noise data. The signal-tonoise ratio is The power of noise is bigger than signal. The red circles correspond to events. Figure shows the PageRank for synthetic data in Z direction. When we have calculated the PageRank for these window pairs, we got to know which windows have high probabilities of being linked, that is, the microseismic events we wanted to detect. In Figure we find the fourth circle and the fifth one are not divided clearly. ecause of the influence by noise, in this direction, we cannot recognize the events from noise. Thus, when the level of noise improves, the PageRank method works but not so robust. In the real problem, this situation often occurs: from one direction or on a single trace, the signal-to-noise ratio is too low to detect. To solve the problem, we can combine the other two directions (X direction, Y direction), or utilize multiple traces for calculation. The theory is the same. Figure 5 presents a case of multiple traces and the detection results. SEG SEG enver nnual Meeting OI Page 65

4 ownloaded /5/ to Redistribution subject to SEG license or copyright; see Terms of Use at - x - Page Rank x - Page Rank x - Page Rank x - Page Rank x - Page Rank Figure 5: PageRank values of multiple traces. Events are marked by red lines. onclusions: We applied the PageRank method for automatic event detection. It helps ranking microseismic events from noisy data. In our test, we divide the data into small windows and correlate these windows and calculate the orrelation oefficient () values. These windows are ranked following their links, and a matrix problem is solved. We use the statistical theory to generate a threshold to distinguish noise and signal. We make a comparison on different thresholds. We tested a single trace and multitrace data, and found that multi-trace data offers more information for detection. The approach presents a great potential for event detection. cknowledgments We thank our research group for helpful advice during this project. SEG SEG enver nnual Meeting OI Page 66

5 ownloaded /5/ to Redistribution subject to SEG license or copyright; see Terms of Use at EITE REFERENES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the SEG Technical Program Expanded bstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. REFERENES guiar,.., and G.. eroza,, pagerank for earthquakes: Seismological Research Letters, 85, no., 5, Gibbons, S. J., and F. Ringdal, 6, The detection of low magnitude seismic events using array-based waveform correlation: Geophysical Journal International, 65, no., 9 66, Kummerow, J.,, Using the value of the crosscorrelation coefficient to locate microseismic events: Geophysics, 75, no., M7 M5, Munro, K.,, utomatic event detection and picking of P-wave arrivals: REWES Research Report, 6,. Maxwell, S.,, Microseismic: Growth born from success: The Leading Edge, 9, 8, Page, L., S. rin, R. Motwani, and T. Winograd, 999, The pagerank citation ranking: ringing order to the web: Technical Report Report, Stanford InfoLab. Song, F., H. S. Kuleli, M. N. Toksöz, E. y, and H. Zhang,, n improved method for hydrofractureinduced microseismic event detection and phase picking: Geophysics, 75, no. 6, 7 5, Warpinski, N., 9, Microseismic monitoring: Inside and out: Journal of Petroleum Technology, 6, no., 8 85, Zhu, L., and L.. Rivera,, note on the dynamic and static displacements from a point source in multilayered media : Geophysical Journal International, 8, no., 69 67, SEG SEG enver nnual Meeting OI Page 67