Changes between Version 7 and Version 8 of glass


Ignore:
Timestamp:
05/29/12 16:21:56 (8 years ago)
Author:
branden
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • glass

    v7 v8  
    199199 
    200200Most of the components of the affinity statistic generated for an origin-pick association are not tunable. The composite affinity statistic is generated via the following equation: 
     201{{{ 
    201202Affinity = Origin.Gap_Affinity * Origin.Number_of_Arrivals_Affinity * 
    202203           Pick.Residual_Affinity * P.Distance_Affinity * Pick.PPD_Affinity 
     204}}} 
     205 
    203206Certain components of the affinity statistic can be turned on and off via the OPCALC_fAffinityStatistics variable. (This variable is currently hardcoded and cannot be changed in a config file.) 
    204207 
    205208The component affinity statistics are generated in the following manner 
    206  
     209{{{ 
    207210Origin.Gap_Affinity 
    208211Value:          Origin.Gap_Affinity = 4 * Bellcurve_Value(Origin.Gap / 360) 
     
    274277               globe, which have no energy from the quake associated with their 
    275278               pick, and may corrupt the origin parameters. 
    276  
    277  
    278  
    279 Composite Affinity 
    280  
    281  
    282  
    283  
     279}}} 
     280 
     281''Composite Affinity'' 
     282 
     283{{{ 
    284284Because the composite affinity statistic is a product of all the components, 
    285285any affinity with a value of 0 automatically nixes the association.  The 
     
    332332pick would not be eligible for assignment to the origin. 
    333333 
    334  
    335  
    336  
    337  
    338 Relevant Origin Time Range 
    339  
    340  
    341  
    342  
     334}}} 
     335 
     336''Relevant Origin Time Range'' 
     337 
     338{{{ 
    343339In it's primary role of trying to assign a new pick to its most likely origin, 
    344340the Associator grabs only the origins that are time relevant to the pick.  The 
     
    346342starts 2400 seconds pick (hardcoded as OPCALC_dSecondaryAssociationPrePickTime) 
    347343prior to the pick and ends at the time of the pick. 
    348  
    349  
    350  
    351  
    352  
    353 Relevant Pick Time Ranges For An Origin 
    354  
    355  
    356  
    357  
     344}}} 
     345 
     346''Relevant Pick Time Ranges For An Origin'' 
     347 
     348{{{ 
    358349When the Associator runs in it's alternate scope of reevaluating how all of 
    359350the time relevant picks associate with an origin after that origin has been 
    360351modified, it uses various time windows.  (this is a bug/feature). 
    361352 
    362  
    363  
    364 Picks not assigned to any origin 
     353* Picks not assigned to any origin 
    365354 
    366355The Associator uses a pick time range of (Origin_Time - 2400 seconds to 
     
    368357OPCALC_dSecondaryAssociationPrePickTime) 
    369358 
    370 Picks assigned to a different origin 
     359* Picks assigned to a different origin 
    371360 
    372361The Associator uses a pick time range of (Origin_Time - 2000 seconds to 
    373362Origin_Time)  (2000 seconds is hardcoded) 
    374  
    375  
    376  
    377  
    378  
    379  
    380 The Locator 
    381  
    382  
    383  
    384  
    385  
    386  
    387 Overview 
    388  
    389  
    390  
    391  
     363}}} 
     364 
     365 
     366 
     367==== The Locator ==== 
     368 
     369===== Overview ===== 
     370 
     371{{{ 
    392372The Locator refines Origin parameters based on a least-squares inversion of 
    393373the residual vectors for associated phases.  (or something like that)  The 
     
    395375pick along with the dtdx(horizontal slowness) and dtdz (vertical slowness) 
    396376for the phase association. 
    397  
    398  
    399  
    400  
    401  
    402 Parameters Tunable 
    403  
    404  
    405  
    406  
     377}}} 
     378 
     379===== Parameters Tunable ===== 
     380 
     381{{{ 
    407382The Locator has two tunable parameters: the location weight of each phase 
    408383(set in the traveltime tables), and the number of iterations it will run 
    409384each time it is invoked. 
    410  
    411  
    412  
    413  
    414  
    415 Location Weight 
    416  
    417  
    418  
    419  
     385}}} 
     386 
     387''Location Weight'' 
     388 
     389{{{ 
    420390The Location Weight is set in the traveltime tables.  See Tunable Affinity 
    421391Statistics Traveltime Tables under the Associator section for a description 
     
    423393Locator to apply that weight to picks of that phase.  The default weight 
    424394is 1.0. 
    425  
    426  
    427  
    428  
    429  
    430 Number of Iterations 
    431  
    432  
    433  
    434  
     395}}} 
     396 
     397''Number of Iterations'' 
     398 
     399{{{ 
    435400A line similar to the following 
    436401 
     
    440405locator to be run for three iterations each time the Locator is invoked. 
    441406The default is 1. 
    442  
    443  
    444  
    445  
    446  
    447 Parameters Untunable 
    448  
    449  
    450  
     407}}} 
     408 
     409===== Parameters Untunable ===== 
    451410 
    452411The Locator has no untunable parameters. 
    453412 
    454  
    455  
    456  
    457  
    458 The Filter 
    459  
    460  
    461  
    462  
    463  
    464  
    465 Overview 
    466  
    467  
    468  
    469  
     413==== The Filter ==== 
     414 
     415===== Overview ===== 
     416{{{ 
    470417The Filter is the portion of Glass that eliminates origins or pick-origin 
    471418associations because they don't pass certain validity tests.  The Filter 
     
    474421pick-associations from the Glass output.  The filters can be thought of as 
    475422falling into 3 groups:  core, network specific, and output. 
    476  
    477  
    478  
    479  
    480  
    481 Core Filters 
    482  
    483  
    484  
    485  
     423}}} 
     424 
     425''Core Filters'' 
     426{{{ 
    486427The Core filters are very simple.  They ensure that existing origins and 
    487428origin/pick associations don't violate the original association requirements. 
     
    491432has dropped due to changes in the parameters of the origin. 
    492433 
    493  
    494  
    495  
    496434Validate Origin 
    497  
    498  
    499435 
    500436A filter will delete an origin if it discovers any of the following 3 
    501437conditions: 
    502  
    503  
    504  
    505438   
    506439Number of phases associated with the origin dropped below the number 
    507440  of points required by the Nucleator to create a new origin. 
    508441 
    509    
    510442Number of phases used by the locator to compute the latest hypocenter 
    511443  for the origin dropped below the number of points required by the Nucleator 
    512444  to create a new origin. 
    513  
    514445   
    515446Number of P phases associated with the origin drops below (number of 
     
    517448  against excessive origin deletion. 
    518449 
    519  
    520  
    521  
    522  
    523450Validate Associations 
    524  
    525  
    526  
    527451 
    528452A filter will delete an origin / pick association (unassign the pick from the 
     
    535459the locator swings the origin in some direction that is not favorable to the 
    536460pick. 
    537  
    538  
    539  
    540  
    541  
    542 Network Specific Filters 
    543  
    544  
    545  
    546  
     461}}} 
     462 
     463''Network Specific Filters'' 
     464{{{ 
    547465Glass contains several origin oriented network specific filters that are 
    548466capable of invalidating a given origin.  The key reason for these filters is 
     
    554472coincident/noise events in North America to be produced. 
    555473 
    556  
    557  
    558474The current set of filters defines 5 regions: 
    559  
    560  
    561  
    562    
     475  
    563476The Inter-mountain West region (46,-104   -  37.5,-113) 
    564477 
    565    
    566478Western hemisphere (85,-140  -  45,-40) 
    567  
    568479   
    569480Northwest quadrant (85,-170  -  0,-40) 
    570  
    571481   
    572482North polar region (Lat > 70) 
    573  
    574483   
    575484Rest of the world 
    576  
    577  
    578  
    579485 
    580486For each of these regions the filter checks for a minimum number of quality 
     
    583489UU, IW, and MB are not allowed to count towards the minimum requirement. 
    584490 
    585  
    586  
    587491If the epicenter of the origin falls in regions 1 - 4, then the origin will 
    588492be invalidated unless one of the following criteria is met:  
    589493 
    590  
    591  
    592    
    59349410 quality P phases are associated with the origin. (10 is hardcoded 
    594495  as MIN_QLTY_ARR) 
    595496 
    596    
    5974974 quality short-range P phases are associated with the origin, where 
    598498  short-range is defined as distance less than 10 degrees. (4 is hardcoded 
    599499  as MIN_QLTY_CLOSEIN_ARR) 
    600  
    601  
    602  
    603500 
    604501If the epicenter of the origin falls outside regions 1-4 (inside region 5), 
     
    607504This is an integral equation, so if number of points is 8 then you need 5, 
    608505if 9 then 5, if 10 then 6, in order to avoid invalidation. 
    609  
    610  
    611  
    612  
    613  
    614 Output Filters 
    615  
    616  
    617  
    618  
     506}}} 
     507 
     508''Output Filters'' 
     509 
     510{{{ 
    619511Glass contains one final set of filters that do not have a dramatic effect 
    620512on the world of Glass, but affect what information that Glass exports to the 
     
    622514that periodically exports event messages to the rest of Hydra. 
    623515 
    624  
    625  
    626  
    627516Minimum Number of Phases 
    628  
    629  
    630  
    631517 
    632518A filter prevents the publication of origins that have fewer than the minimum 
     
    636522world. 
    637523 
    638  
    639  
    640  
    641524Ancillary Phases 
    642  
    643  
    644  
    645525 
    646526A filter tries to prevent certain coincidental associations with an origin, 
     
    652532for that origin, but all 4 of the PKPdf phases will be published. 
    653533 
    654  
    655  
    656  
    657534Pdif Phases 
    658  
    659  
    660535 
    661536A filter prevents Pdif phases > 110 degrees from being published.  I 
     
    663538degrees for the current picker are almost always coincidental associations. 
    664539 
    665  
    666  
    667  
    668540Old Events 
    669  
    670  
    671541 
    672542A filter prevents quakes older than a threshold number of days from being 
     
    679549beyond the 5.0 day threshold. 
    680550 
    681  
    682  
    683  
    684551Dead Events 
    685  
    686  
    687552 
    688553A filter prevents an origin from being published once it is dead.  The 
     
    694559automatic processing and is being handled further down the processing stream. 
    695560 
    696  
    697  
    698  
    699561Publication Rate 
    700  
    701  
    702562 
    703563The publication module controls how often glass updates the outside world on 
     
    707567the glass_pub_params.d config file, and the controlling parameters are 
    708568documented within that file. 
    709  
    710  
    711  
    712  
    713  
    714 Other Stuff 
    715  
    716  
    717  
    718  
     569}}} 
     570 
     571==== Other Stuff ==== 
     572{{{ 
    719573There are several other mechanisms within Glass that affect processing, but do 
    720574not affect output to the degree of the major four.  The Depth Refocus module 
    721575attempts to discover alternate hypocenter depths that elude the Locator.  The 
    722576State Processing orchestrates how origins and picks are processed by Glass. 
    723  
    724  
    725  
    726  
    727 Depth Refocus 
    728  
    729  
    730  
     577}}} 
     578 
     579===== Depth Refocus ===== 
     580{{{ 
    731581Depth Refocus is a module that searches for an improved origin depth and 
    732582time that cannot be discovered by the Locator.  Refocus assumes a correct 
     
    737587time/depth grid.  It was created to get the Locator out of local residual 
    738588minima.  Currently there are no tunable parameters in the Focus module. 
    739  
    740  
    741  
    742  
    743 State Processing 
    744  
    745  
    746  
     589}}} 
     590 
     591===== State Processing ===== 
     592{{{ 
    747593Glass contains a state processing engine that controls how picks and origins 
    748594are processed.  The engine is entity oriented where an entity is either a 
     
    751597considered to be an engineering issue and is not meant to be modified by a 
    752598seismologist/network operator. 
    753  
    754  
    755  
    756  
    757  
    758 Hydra Administrator 
    759  
    760  
    761  
    762  
     599}}} 
     600 
     601=== Hydra Administrator === 
     602{{{ 
    763603This section to be filled in later with information regarding how to configure 
    764604Glass to run in a Hydra/Earthworm near real-time environment. 
    765  
    766  
    767  
    768  
    769  
    770 Configuring Glass 
    771  
    772  
    773  
    774  
    775  
    776  
    777 glass.d 
    778  
    779  
    780  
    781  
    782  
    783  
    784 Error Reporting 
    785  
    786  
    787  
     605}}} 
     606 
     607==== Configuring Glass ==== 
     608 
     609===== glass.d ===== 
     610 
     611''Error Reporting'' 
    788612 
    789613Glass, by default writes 3 levels of status messages: 
    790  
    791  
    792  
    793    
     614{{{ 
    794615MAJOR ERROR 
    795616  There are 39 Major Errors.  They include things like: 
    796  
    797    
    798  
    799617     
    800618Can't open config file. 
    801  
    802619     
    803620Can't load DLL (Module) 
    804  
    805621     
    806622Can't allocate memory 
    807623 
    808      
    809624Can't load traveltime tables 
    810  
    811      
     625   
    812626Can't get EW message types 
    813  
    814      
     627    
    815628Can't  access required module 
    816  
    817      
     629    
    818630Can't attach to shared memory 
    819631 
    820      
    821632Can't create internal messages 
    822  
    823633     
    824634Can't parse config file command 
    825  
    826635     
    827636Can't write to temporary files. 
    828  
    829    
    830  
    831637   
    832638MAJOR WARNING 
     
    836642  morning. 
    837643 
    838    
    839644MINOR ERROR 
    840645  There are somewhere around 90 Minor Errors.  They cover everything under the 
    841646  sun, including Missed Messages, and so forth so on, and are a general catch 
    842647  all for stuff that the operator or programmer may want to know about. 
    843  
    844  
    845  
    846648 
    847649To make matters more complicated, while the internal error levels are 
     
    850652recommend it) configure glass via the .d file, to issue status messages for 
    851653each DEBUG/INFO message and not issue them for any of the Major Errors. 
    852  
    853  
    854  
     654}}} 
    855655From the default glass.d file: 
    856  
    857  
     656{{{ 
    858657#  DebugLevel commands 
    859658#  CatalogDebugLevel | EarthwormDebugLevel | GlassDebugLevel | GlintDebugLevel | LocatorDebugLevel | PublisherDebugLevel 
     
    888687#  Example: 
    889688#  EarthwormDebugLevel  I:Level=2 I:OTF=1 I:OTD=0 I:OTE=0 I:OTS=0 I:OSM=0 
    890  
    891  
    892  
    893  
    894  
    895 User 
    896  
    897  
    898  
    899  
     689}}} 
     690 
     691=== User === 
     692{{{ 
    900693This section to be filled in later with information regarding how to operate 
    901694the Glass displays including the ManQuake display for manually testing a 
    902695hypothetical hypocenter. 
    903  
    904  
    905  
    906  
    907  
    908 Programmer 
    909  
    910  
    911  
    912  
     696}}} 
     697 
     698=== Programmer === 
     699{{{ 
    913700This section to be filled in later with information regarding the 
    914701architectural layout of Glass, along with notes on how to maintain/modify it. 
     702}}} 
    915703 
    916704